Status of single-cell RNA sequencing for reproductive toxicology in zebrafish and the transcriptomic trade-off.

The utilization of transcriptomic studies identifying profiles of gene expression, especially in toxicogenomics, has catapulted next-generation sequencing to the forefront of reproductive toxicology. An innovative yet underutilized RNA sequencing technique emerging into this field is single-cell RNA sequencing (scRNA-seq), which provides sequencing at the individual cellular level of gonad tissue. ScRNA-seq provides a novel and unique perspective for identifying distinct cellular profiles, including identification of rare cell subtypes. The specificity of scRNA-seq is a powerful tool for reproductive toxicity research, especially for translational animal models including zebrafish. Studies to date not only have focused on 'tissue atlassing' or characterizing what cell types make up different tissues but have also begun to include toxicant exposure as a factor that this review aims to explore. Future scRNA-seq studies will contribute to understanding exposure-induced outcomes; however, the trade-offs with traditional methods need to be considered.

Adult-Onset Transcriptomic Effects of Developmental Exposure to Benzene in Zebrafish (Danio rerio): Evaluating a Volatile Organic Compound of Concern.

Urban environments are afflicted by mixtures of anthropogenic volatile organic compounds (VOCs). VOC sources that drive human exposure include vehicle exhaust, industrial emissions, and oil spillage. The highly volatile VOC benzene has been linked to adverse health outcomes. However, few studies have focused on the later-in-life effects of low-level benzene exposure during the susceptible window of early development. Transcriptomic responses during embryogenesis have potential long-term consequences at levels equal to or lower than 1 ppm, therefore justifying the analysis of adult zebrafish that were exposed during early development. Previously, we identified transcriptomic alteration following controlled VOC exposures to 0.1 or 1 ppm benzene during the first five days of embryogenesis using a zebrafish model. In this study, we evaluated the adult-onset transcriptomic responses to this low-level benzene embryogenesis exposure (n = 20/treatment). We identified key genes, including col1a2 and evi5b, that were differentially expressed in adult zebrafish in both concentrations. Some DEGs overlapped at the larval and adult stages, specifically nfkbiaa, mecr, and reep1. The observed transcriptomic results suggest dose- and sex-dependent changes, with the highest impact of benzene exposure to be on cancer outcomes, endocrine system disorders, reproductive success, neurodevelopment, neurological disease, and associated pathways. Due to molecular pathways being highly conserved between zebrafish and mammals, developmentally exposed adult zebrafish transcriptomics is an important endpoint for providing insight into the long term-effects of VOCs on human health and disease.

Evaluating Phenotypic and Transcriptomic Responses Induced by Low-Level VOCs in Zebrafish: Benzene as an Example.

Urban environments are plagued by complex mixtures of anthropogenic volatile organic compounds (VOCs), such as mixtures of benzene, toluene, ethylene, and xylene (BTEX). Sources of BTEX that drive human exposure include vehicle exhaust, industrial emissions, off-gassing of building material, as well as oil spillage and leakage. Among the BTEX mixture, benzene is the most volatile compound and has been linked to numerous adverse health outcomes. However, few studies have focused on the effects of low-level benzene on exposure during early development, which is a susceptible window when hematological, immune, metabolic, and detoxification systems are immature. In this study, we used zebrafish to conduct a VOC exposure model and evaluated phenotypic and transcriptomic responses following 0.1 and 1 ppm benzene exposure during the first five days of embryogenesis (n = 740 per treatment). The benzene body burden was 2 mg/kg in 1 ppm-exposed larval zebrafish pools and under the detection limit in 0.1 ppm-exposed fish. No observable phenotypic changes were found in both larvae except for significant skeletal deformities in 0.1 ppm-exposed fish (p = 0.01) compared with unexposed fish. Based on transcriptomic responses, 1 ppm benzene dysregulated genes that were implicated with the development of hematological system, and the regulation of oxidative stress response, fatty acid metabolism, immune system, and inflammatory response, including apob, nfkbiaa, serpinf1, foxa1, cyp2k6, and cyp2n13 from the cytochrome P450 gene family. Key genes including pik3c2b, pltp, and chia.2 were differentially expressed in both 1 and 0.1 ppm exposures. However, fewer transcriptomic changes were induced by 0.1 ppm compared with 1 ppm. Future studies are needed to determine if these transcriptomic responses during embryogenesis have long-term consequences at levels equal to or lower than 1 ppm.

Insight into 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced disruption of zebrafish spermatogenesis via single cell RNA-seq.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a potent and environmentally persistent endocrine disrupting chemical. Our previous work demonstrated the latent reproductive maladies of early-life TCDD exposure in zebrafish. Zebrafish acutely exposed to low, environmentally relevant levels of TCDD (50 pg/mL) during two windows of sexual differentiation in development (1 hour of exposure at 3 and 7 weeks postfertilization) were later infertile, showed a reduction in sperm, and exhibited gene expression consistent with an altered microenvironment, even months after exposure. Due to the highly heterogeneous cell- type and -stage landscape of the testes, we hypothesized various cell types contribute markedly different profiles toward the pathology of TCDD exposure. To investigate the contributions of the diverse cell types in the adult zebrafish testes to TCDD-induced pathology, we utilized single-cell RNA-seq and the 10x Genomics platform. The method successfully captured every stage of testicular germ cell development. Testes of adult fish exposed during sexual differentiation to TCDD contained sharply decreased populations of late spermatocytes, spermatids, and spermatozoa. Spermatogonia and early spermatocyte populations were, in contrast, enriched following exposure. Pathway analysis of differentially expressed genes supported previous findings that TCDD exposure resulted in male infertility, and suggested this outcome is due to apoptosis of spermatids and spermatozoa, even years after exposure cessation. Increased germ cell apoptosis was confirmed histologically. These results provide support for an environmental exposure explanation of idiopathic male infertility.

Multi- and Transgenerational Effects of Developmental Exposure to Environmental Levels of PFAS and PFAS Mixture in Zebrafish (Danio rerio).

Per- and polyfluoroalkyl substances (PFASs) are ubiquitous in the environment and are tied to myriad health effects. Despite the phasing out of the manufacturing of two types of PFASs (perfluorosulfonic acid (PFOS) and perfluorooctanoic acid (PFOA)), chemical composition renders them effectively indestructible by ambient environmental processes, where they thus remain in water. Exposure via water can affect both human and aquatic wildlife. PFASs easily cross the placenta, exposing the fetus at critical windows of development. Little is known about the effects of low-level exposure during this period; even less is known about the potential for multi- and transgenerational effects. We examined the effects of ultra-low, very low, and low-level PFAS exposure (7, 70, and 700 ng/L PFOA; 24, 240, 2400 ng/L PFOS; and stepwise mixtures) from 0-5 days post-fertilization (dpf) on larval zebrafish (Danio rerio) mortality, morphology, behavior and gene expression and fecundity in adult F0 and F1 fish. As expected, environmentally relevant PFAS levels did not affect survival. Morphological abnormalities were not observed until the F1 and F2 generations. Behavior was affected differentially by each chemical and generation. Gene expression was increasingly perturbed in each generation but consistently showed lipid pathway disruption across all generations. Dysregulation of behavior and gene expression is heritable, even in larvae with no direct or indirect exposure. This is the first report of the transgenerational effects of PFOA, PFOS, and their mixture in terms of zebrafish behavior and untargeted gene expression.

Comparative Toxicotranscriptomics of Single Cell RNA-Seq and Conventional RNA-Seq in TCDD-Exposed Testicular Tissue.

In this report, we compare the outcomes and limitations of two methods of transcriptomic inquiry on adult zebrafish testes exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) during sexual differentiation: conventional or bulk RNA-seq (bulk-seq) and single cell RNA sequencing (scRNA-seq) data. scRNA-seq has emerged as a valuable tool for uncovering cell type-specific transcriptome dynamics which exist in heterogeneous tissue. Our lab previously showed the toxicological value of the scRNA-seq pipeline to characterize the sequelae of TCDD exposure in testes, demonstrating that loss of spermatids and spermatozoa, but not other cell types, contributed to the pathology of infertility in adult male zebrafish exposed during sexual differentiation. To investigate the potential for technical artifacts in scRNA-seq such as cell dissociation effects and reduced transcriptome coverage, we compared bulk-sequenced and scRNA-seq-paired samples from control and TCDD-exposed samples to understand what is gained and lost in scRNA-seq vs bulk-seq, both transcriptomically and toxicologically. We hypothesized that the testes may be sensitive to tissue disruption as they contain multiple cell types under constant division and/or maturation, and that TCDD exposure may mediate the extent of sensitivity. Thus, we sought to understand the extent to which this dissociation impacts the toxicological value of data returned from scRNA-seq. We confirm that the required dissociation of individual cells from intact tissue has a significant impact on gene expression, affecting gene pathways with the potential to confound toxicogenomics studies on exposures if findings are not well-controlled and well-situated in context. Additionally, a common scRNA-seq method using cDNA amplified from the 3' end of mRNA under-detects low-expressing transcripts including transcription factors. We confirm this, and show TCDD-related genes may be overlooked by scRNA-seq, however, this under-detection effect is not mediated by TCDD exposure. Even so, scRNA-seq generally extracted toxicologically relevant information better than the bulk-seq method in the present study. This report aims to inform future experimental design for transcriptomic investigation in the growing field of toxicogenomics by demonstrating the differential information extracted from sequencing cells-despite being from the same tissue and exposure scheme-is influenced by the specific protocol used, with implications for the interpretation of exposure-induced risk.

Developmental Phenotypic and Transcriptomic Effects of Exposure to Nanomolar Levels of 4-Nonylphenol, Triclosan, and Triclocarban in Zebrafish (Danio rerio).

Triclosan, triclocarban and 4-nonylphenol are all chemicals of emerging concern found in a wide variety of consumer products that have exhibited a wide range of endocrine-disrupting effects and are present in increasing amounts in groundwater worldwide. Results of the present study indicate that exposure to these chemicals at critical developmental periods, whether long-term or short-term in duration, leads to significant mortality, morphologic, behavioral and transcriptomic effects in zebrafish (Danio rerio). These effects range from total mortality with either long- or short-term exposure at 100 and 1000 nM of triclosan, to abnormalities in uninflated swim bladder seen with long-term exposure to triclocarban and short-term exposure to 4-nonylphenol, and cardiac edema seen with short-term 4-nonylphenol exposure. Additionally, a significant number of genes involved in neurological and cardiovascular development were differentially expressed after the exposures, as well as lipid metabolism genes and metabolic pathways after exposure to each chemical. Such changes in behavior, gene expression, and pathway abnormalities caused by these three known endocrine disruptors have the potential to impact not only the local ecosystem, but human health as well.

Adipocyte-driven unfolded protein response is a shared transcriptomic signature of metastatic prostate carcinoma cells.

A critical unknown in the field of skeletal metastases is how cancer cells find a way to thrive under harsh conditions, as exemplified by metastatic colonization of adipocyte-rich bone marrow by prostate carcinoma cells. To begin understanding molecular processes that enable tumor cells to survive and progress in difficult microenvironments such as bone, we performed unbiased examination of the transcriptome of two different prostate cancer cell lines in the absence or presence of bone marrow adipocytes. Our RNAseq analyses and subsequent quantitative PCR and protein-based assays reveal that upregulation of endoplasmic reticulum (ER) stress and unfolded protein response (UPR) genes is a shared signature between metastatic prostate carcinoma cell lines of different origin. Pathway analyses and pharmacological examinations highlight the ER chaperone BIP as an upstream coordinator of this transcriptomic signature. Additional patient-based data support our overall conclusion that ER stress and UPR induction are shared, important factors in the response and adaptation of metastatic tumor cells to their micro-environment. Our studies pave the way for additional mechanistic investigations and offer new clues towards effective therapeutic interventions in metastatic disease.