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.

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.

Point-of-use carbon-block drinking water filters change gut microbiome of larval zebrafish.

Activated carbon block (ACB) point-of-use (PoU) drinking water filters can change the bacterial composition in drinking water. Consuming ACB PoU filtered water may also influence gut microbiomes. This study uses the zebrafish model to evaluate how the ACB PoU filter affects the gut microbiomes and phenotypic responses in larvae and adulthood. An ACB PoU filter manifold system was constructed to feed larval and adult zebrafish tap and filtered water at the early and late stages of the filter operation period. Adult zebrafish gut microbiomes were not affected by exposure to water types and filter stages. Unlike the adult, gut microbiomes of the larvae exposed to filtered water at the late stage of filter operation were dominated by more filter-relevant bacterial taxa, including Comamonadaceae and Brevundimonas, than the early stage-filtered-water- and tap water-exposed larvae. We also found some fish that were either exposed to filtered water at early and late stages or tap water supplied to the filter toward the end of the experiment showed hyperactive locomotion behaviour, and had significantly lower relative abundances of a Pseudomonas spp. (OTU3) than the normally behaved fish. Our findings indicate that ACB PoU filtered water can alter gut microbiomes and affect the behaviour patterns in larval zebrafish.

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.

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.

Developmental phenotypic and transcriptomic effects of exposure to nanomolar levels of metformin in zebrafish.

Metformin is found in the majority of lakes and streams in the United States, leading to widespread environmental exposure. Results of the present study indicate that extended duration metformin exposure at critical developmental periods leads to decreased survival rates in zebrafish (danio rerio), an NIH approved human model. Significant abnormalities are seen with extended duration metformin exposure from 4 h post fertilization up to 5 days post fertilization, although short term metformin exposure for 24 h at 4-5 days post fertilization did not lead to any significant abnormalities. Both extended and short term duration did however have an impact on locomotor activity of zebrafish, and several genes involved in neurological and cardiovascular development were differentially expressed after exposure to metformin. The changes seen in behavior, gene expression and morphological abnormalities caused by metformin exposure should be examined further in future studies in order to assess their potential human health implications as metformin prescriptions continue to increase worldwide.

The phenotypic and transcriptomic effects of developmental exposure to nanomolar levels of estrone and bisphenol A in zebrafish.

Estrone and BPA are two endocrine disrupting chemicals (EDCs) that are predicted to be less potent than estrogens such as 17ß-estradiol and 17α-ethinylestradiol. Human exposure concentrations to estrone and BPA can be as low as nanomolar levels. However, very few toxicological studies have focused on the nanomolar-dose effects. Low level of EDCs can potentially cause non-monotonic responses. In addition, exposures at different developmental stages can lead to different health outcomes. To identify the nanomolar-dose effects of estrone and BPA, we used zebrafish modeling to study the phenotypic and transcriptomic responses after extended duration exposure from 0 to 5 days post-fertilization (dpf) and short-term exposure at days 4-5 post fertilization. We found that non-monotonic transcriptomic responses occurred after extended duration exposures at 1 nM of estrone or BPA. At this level, estrone also caused hypoactivity locomotive behavior in zebrafish. After both extended duration and short-term exposures, BPA led to more apparent phenotypic responses, i.e. skeletal abnormalities and locomotion changes, and more significant transcriptomic responses than estrone exposure. After short-term exposure, BPA at concentrations equal or above 100 nM affected locomotive behavior and changed the expression of both estrogenic and non-estrogenic genes that are linked to neurological diseases. These data provide gaps of mechanisms between neurological genes expression and associated phenotypic response due to estrone or BPA exposures. This study also provides insights for assessing the acceptable concentration of BPA and estrone in aquatic environments.

Detection of endocrine disrupting chemicals in Danio rerio and Daphnia pulex: Step-one, behavioral screen.

Anthropogenic surface and ground water contamination by chemicals is a global problem, and there is an urgent need to develop tools to identify and elucidate biological effects. Contaminants of emerging concern (CECs) are not typically monitored or regulated and those with known or suspected endocrine disrupting potential have been termed endocrine disrupting chemicals (EDCs). Many CECs are known to be neurotoxic (e.g., insecticides) and many are incompletely characterized. Behavioral responses can identify chemicals with neuroactive properties, which can be relevant to EDC mechanisms (e.g., neuroendocrine disturbances). Two freshwater species, Daphnia pulex and Danio rerio, were evaluated for swimming behavior alterations resulting from 24-hr exposure to 9 CECs: triclosan, triclocarban, chlorpyrifos, dieldrin, 4-nonylphenol, bisphenol-A, atrazine, metformin, and estrone. This is the first step in the development of a bioassay for detecting estrogenic and/or anti-androgenic activity with the goal to evaluate complex mixtures of uncharacterized contaminants in water samples. The second step, described in a subsequent report, examines transcriptome alterations following chemical exposure. Significant differences in the swimming behavior response and sensitivity were found across chemicals within a species and across species for a given chemical in this unique optical bioassay system. In the concentration ranges studied, significant behavioral alterations were detected for 6 of 9 CECs for D. pulex and 4 of 9 CECs for D. rerio. These results underscore the utility of this bioassay to identify behavioral effects of sublethal concentrations of CECs before exploration of transcriptomic alterations for EDC detection.

TCDD-induced multi- and transgenerational changes in the methylome of male zebrafish gonads.

The legacy endocrine disrupting chemical and aryl hydrocarbon receptor agonist, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), is produced as a byproduct of industrial processes and causes adverse health effects ranging from skin irritation to cancer. TCDD endpoints are also observed in subsequent, unexposed generations; however, the mechanisms of these multi- and transgenerational effects are unknown. We hypothesized an epigenetic mechanism, specifically DNA methylation for the transgenerational, male-mediated reproductive effects of developmental TCDD exposure. Using whole genome bisulfite sequencing, we evaluated DNA methylation changes in three generations of zebrafish, the first of which was exposed to TCDD during sexual development at 50 ppt for 1 h at both 3- and 7-week post-fertilization. We discovered that TCDD induces multi- and transgenerational methylomic changes in testicular tissue from zebrafish with decreased reproductive capacity, but most significantly in the indirectly exposed F1 generation. In comparing differentially methylated genes to concurrent transcriptomic changes, we identified several genes and pathways through which transgenerational effects of low level TCDD exposure are likely inherited. These include significant differential methylation of genes involved in reproduction, endocrine function, xenobiotic metabolism, and epigenetic processing. Notably, a number of histone modification genes were both differentially methylated and expressed in all generations, and many differentially methylated genes overlapped between multiple generations. Collectively, our results suggest that DNA methylation is a promising mechanism to explain male-mediated transgenerational reproductive effects of TCDD exposure in zebrafish, and these effects are likely inherited through integration of multiple epigenetic pathways.

Developmental exposure to Pb2+ induces transgenerational changes to zebrafish brain transcriptome.

Lead (Pb2+) is a major public health hazard for urban children, with profound and well-characterized developmental and behavioral implications across the lifespan. The ability of early Pb2+ exposure to induce epigenetic changes is well-established, suggesting that Pb2+-induced neurobehavioral deficits may be heritable across generations. Understanding the long-term and multigenerational repercussions of lead exposure is crucial for clarifying both the genotypic alterations behind these behavioral outcomes and the potential mechanism of heritability. To study this, zebrafish (Danio rerio) embryos (<2 h post fertilization; EK strain) were exposed for 24 h to waterborne Pb2+ at a concentration of 10 µM. This exposed F0 generation was raised to adulthood and spawned to produce the F1 generation, which was subsequently spawned to produce the F2 generation. Previous avoidance conditioning studies determined that a 10 µM Pb2+ dose resulted in learning impairments persisting through the F2 generation. RNA was extracted from control- and 10 µM Pb2+-lineage F2 brains, (n = 10 for each group), sequenced, and transcript expression was quantified utilizing Quant-Seq. 648 genes were differentially expressed in the brains of F2 lead-lineage fish versus F2 control-lineage fish. Pathway analysis revealed altered genes in processes including synaptic function and plasticity, neurogenesis, endocrine homeostasis, and epigenetic modification, all of which are implicated in lead-induced neurobehavioral deficits and/or their inheritance. These data will inform future investigations to elucidate the mechanism of adult-onset and transgenerational health effects of developmental lead exposure.

Nanoplastics impact the zebrafish (Danio rerio) transcriptome: Associated developmental and neurobehavioral consequences.

Microplastics (MPs) are a ubiquitous pollutant detected not only in marine and freshwater bodies, but also in tap and bottled water worldwide. While MPs have been extensively studied, the toxicity of their smaller counterpart, nanoplastics (NPs), is not well documented. Despite likely large-scale human and animal exposure to NPs, the associated health risks remain unclear, especially during early developmental stages. To address this, we investigated the health impacts of exposures to both 50 and 200 nm polystyrene NPs in larval zebrafish. From 6 to 120 h post-fertilization (hpf), developing zebrafish were exposed to a range of fluorescent NPs (10-10,000 parts per billion). Dose-dependent increases in accumulation were identified in exposed larval fish, potentially coinciding with an altered behavioral response as evidenced through swimming hyperactivity. Notably, exposures did not impact mortality, hatching rate, or deformities; however, transcriptomic analysis suggests neurodegeneration and motor dysfunction at both high and low concentrations. Furthermore, results of this study suggest that NPs can accumulate in the tissues of larval zebrafish, alter their transcriptome, and affect behavior and physiology, potentially decreasing organismal fitness in contaminated ecosystems. The uniquely broad scale of this study during a critical window of development provides crucial multidimensional characterization of NP impacts on human and animal health.

Roles for a lipid phosphatase in the activation of its opposing lipid kinase.

Fig4 is a phosphoinositide phosphatase that converts PI3,5P2 to PI3P. Paradoxically, mutation of Fig4 results in lower PI3,5P2, indicating that Fig4 is also required for PI3,5P2 production. Fig4 promotes elevation of PI3,5P2, in part, through stabilization of a protein complex that includes its opposing lipid kinase, Fab1, and the scaffold protein Vac14. Here we show that multiple regions of Fig4 contribute to its roles in the elevation of PI3,5P2: its catalytic site, an N-terminal disease-related surface, and a C-terminal region. We show that mutation of the Fig4 catalytic site enhances the formation of the Fab1-Vac14-Fig4 complex, and reduces the ability to elevate PI3,5P2. This suggests that independent of its lipid phosphatase function, the active site plays a role in the Fab1-Vac14-Fig4 complex. We also show that the N-terminal disease-related surface contributes to the elevation of PI3,5P2 and promotes Fig4 association with Vac14 in a manner that requires the Fig4 C-terminus. We find that the Fig4 C-terminus alone interacts with Vac14 in vivo and retains some functions of full-length Fig4. Thus, a subset of Fig4 functions are independent of its phosphatase domain and at least three regions of Fig4 play roles in the function of the Fab1-Vac14-Fig4 complex.

Unanchored ubiquitin chains do not lead to marked alterations in gene expression in Drosophila melanogaster.

The small protein modifier ubiquitin regulates various aspects of cellular biology through its chemical conjugation onto proteins. Ubiquitination of proteins presents itself in numerous iterations, from a single mono-ubiquitination event to chains of poly-ubiquitin. Ubiquitin chains can be attached onto other proteins or can exist as unanchored species, i.e. free from another protein. Unanchored ubiquitin chains are thought to be deleterious to the cell and rapidly disassembled into mono-ubiquitin. We recently examined the toxicity and utilization of unanchored poly-ubiquitin in Drosophila melanogaster We found that free poly-ubiquitin species are largely innocuous to flies and that free poly-ubiquitin can be controlled by being degraded by the proteasome or by being conjugated onto another protein as a single unit. Here, to explore whether an organismal defense is mounted against unanchored chains, we conducted RNA-Seq analyses to examine the transcriptomic impact of free poly-ubiquitin in the fly. We found ∼90 transcripts whose expression is altered in the presence of different types of unanchored poly-ubiquitin. The set of genes identified was essentially devoid of ubiquitin-, proteasome-, or autophagy-related components. The seeming absence of a large and multipronged response to unanchored poly-ubiquitin supports the conclusion that these species need not be toxic in vivo and underscores the need to re-examine the role of free ubiquitin chains in the cell.

Developmental Dioxin Exposure Alters the Methylome of Adult Male Zebrafish Gonads.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental toxicant and endocrine disrupting compound with reproductive and developmental effects in humans and model organisms, including zebrafish. Our previous microarray and histological studies found defects in spermatogenesis and fertility of zebrafish in response to acute developmental TCDD exposure. These effects are apparent following exposure during reproductive development, modeling fetal basis of adult-onset disease. Some outcomes of these previous studies (reduced fertility, changes in sex ratio, transcriptomic alterations) are also transgenerational - persisting to unexposed generations - through the male germline. We hypothesized that DNA methylation could be a possible mechanism for these reproductive effects and performed whole genome bisulfite sequencing (WGBS), which identifies whole genome DNA methylation status at the base pair level, on testes of adult zebrafish exposed to TCDD (two separate hour-long exposures to 50 pg/mL TCDD at 3 and 7 weeks post fertilization). In response to TCDD exposure, multiple genes were differentially methylated; many of which are involved in reproductive processes or epigenetic modifications, suggesting a role of DNA methylation in later-life health outcomes. Additionally, several differentially methylated genes corresponded with gene expression changes identified in TCDD-exposed zebrafish testes, indicating a potential link between DNA methylation and gene expression. Ingenuity pathway analysis of WGBS and microarray data revealed genes involved in reproductive processes and development, RNA regulation, the cell cycle, and cellular morphology and development. We conclude that site-specific changes in DNA methylation of adult zebrafish testes occur in response to acute developmental TCDD exposure.