Gene transcription ontogeny of hypothalamic-pituitary-thyroid axis development in early-life stage fathead minnow and zebrafish.

The hypothalamic-pituitary-thyroid (HPT) axis is known to play a crucial role in the development of teleost fish. However, knowledge of endogenous transcription profiles of thyroid-related genes in developing teleosts remains fragmented. We selected two model teleost species, the fathead minnow (Pimephales promelas) and the zebrafish (Danio rerio), to compare the gene transcription ontogeny of the HPT axis. Control organisms were sampled at several time points during embryonic and larval development until 33 days post-fertilization. Total RNA was extracted from pooled, whole fish, and thyroid-related mRNA expression was evaluated using quantitative polymerase chain reaction. Gene transcripts examined included: thyrotropin-releasing hormone receptor (trhr), thyroid-stimulating hormone receptor (tshr), sodium-iodide symporter (nis), thyroid peroxidase (tpo), thyroglobulin (tg), transthyretin (ttr), deiodinases 1, 2, 3a, and 3b (dio1, dio2, dio3a and 3b), and thyroid hormone receptors alpha and beta (thrα and ß). A loess regression method was successful in identifying maxima and minima of transcriptional expression during early development of both species. Overall, we observed great similarities between the species, including maternal transfer, at least to some extent, of almost all transcripts (confirmed in unfertilized eggs), increasing expression of most transcripts during hatching and embryo-larval transition, and indications of a fully functional HPT axis in larvae. These data will aid in the development of hypotheses on the role of certain genes and pathways during development. Furthermore, this provides a background reference dataset for designing and interpreting targeted transcriptional expression studies both for fundamental research and for applications such as toxicology.

Re-evaluating the Significance of Estrone as an Environmental Estrogen.

Studies worldwide have demonstrated the occurrence of feminized male fish at sites impacted by human and animal wastes. A variety of chemicals could contribute to this phenomenon, but those receiving the greatest attention in terms of research and monitoring have been 17ß-estradiol (ß-E2) and 17α-ethinylestradiol, due both to their prevalence in the environment and strong estrogenic potency. A third steroid, estrone (E1), also can occur at high concentrations in surface waters but generally has been of lesser concern due to its relatively lower affinity for vertebrate estrogen receptors. In an initial experiment, male fathead minnow (Pimephales promelas) adults were exposed for 4-d to environmentally relevant levels of waterborne E1, which resulted in plasma ß-E2 concentrations similar to those found in reproductively active females. In a second exposure we used 13C-labeled E1, together with liquid chromatography-tandem mass spectrometry, to demonstrate that elevated ß-E2 measured in the plasma of the male fish was indeed derived from the external environment, most likely via a conversion catalyzed by one or more 17ß-hydroxysteroid dehydrogenases. The results of our studies suggest that the potential impact of E1 as an environmental estrogen currently is underestimated.

ETV6 Deficiency Unlocks ERG-Dependent Microsatellite Enhancers to Drive Aberrant Gene Activation in B-Lymphoblastic Leukemia.

Distal enhancers play critical roles in sustaining oncogenic gene-expression programs. We identify aberrant enhancer-like activation of GGAA tandem repeats as a characteristic feature of B-cell acute lymphoblastic leukemia (B-ALL) with genetic defects of the ETV6 transcriptional repressor, including ETV6-RUNX1+ and ETV6-null B-ALL. We show that GGAA repeat enhancers are direct activators of previously identified ETV6-RUNX1+/- like B-ALL "signature" genes, including the likely leukemogenic driver EPOR. When restored to ETV6-deficient B-ALL cells, ETV6 directly binds to GGAA repeat enhancers, represses their acetylation, downregulates adjacent genes, and inhibits B-ALL growth. In ETV6-deficient B-ALL cells, we find that the ETS transcription factor ERG directly binds to GGAA microsatellite enhancers and is required for sustained activation of repeat enhancer-activated genes. Together, our findings reveal an epigenetic gatekeeper function of the ETV6 tumor suppressor gene and establish microsatellite enhancers as a key mechanism underlying the unique gene-expression program of ETV6-RUNX1+/- like B-ALL. SIGNIFICANCE: We find a unifying mechanism underlying a leukemia subtype-defining gene-expression signature that relies on repetitive elements with poor conservation between humans and rodents. The ability of ETV6 to antagonize promiscuous, nonphysiologic ERG activity may shed light on other roles of these key regulators in hematolymphoid development and human disease. See related commentary by Mercher, p. 2. This article is highlighted in the In This Issue feature, p. 1.

Selective Enhancer Dependencies in MYC -Intact and MYC -Rearranged Germinal Center B-cell Diffuse Large B-cell Lymphoma.

High expression of MYC and its target genes define a subset of germinal center B-cell diffuse large B-cell lymphoma (GCB-DLBCL) associated with poor outcomes. Half of these high-grade cases show chromosomal rearrangements between the MYC locus and heterologous enhancer-bearing loci, while focal deletions of the adjacent non-coding gene PVT1 are enriched in MYC -intact cases. To identify genomic drivers of MYC activation, we used high-throughput CRISPR-interference (CRISPRi) profiling of candidate enhancers in the MYC locus and rearrangement partner loci in GCB-DLBCL cell lines and mantle cell lymphoma (MCL) comparators that lacked common rearrangements between MYC and immunoglobulin (Ig) loci. Rearrangements between MYC and non-Ig loci were associated with unique dependencies on specific enhancer subunits within those partner loci. Notably, fitness dependency on enhancer modules within the BCL6 super-enhancer ( BCL6 -SE) cluster regulated by a transcription factor complex of MEF2B, POU2F2, and POU2AF1 was higher in cell lines bearing a recurrent MYC::BCL6 -SE rearrangement. In contrast, GCB-DLBCL cell lines without MYC rearrangement were highly dependent on a previously uncharacterized 3' enhancer within the MYC locus itself (GCBME-1), that is regulated in part by the same triad of factors. GCBME-1 is evolutionarily conserved and active in normal germinal center B cells in humans and mice, suggesting a key role in normal germinal center B cell biology. Finally, we show that the PVT1 promoter limits MYC activation by either native or heterologous enhancers and demonstrate that this limitation is bypassed by 3' rearrangements that remove PVT1 from its position in cis with the rearranged MYC gene. Key points: CRISPR-interference screens identify a conserved germinal center B cell MYC enhancer that is essential for GCB-DLBCL lacking MYC rearrangements. Functional profiling of MYC partner loci reveals principles of MYC enhancer-hijacking activation by non-immunoglobulin rearrangements.

Case Study in 21st-Century Ecotoxicology: Using In Vitro Aromatase Inhibition Data to Predict Reproductive Outcomes in Fish In Vivo.

To reduce the use of intact animals for chemical safety testing, while ensuring protection of ecosystems and human health, there is a demand for new approach methodologies (NAMs) that provide relevant scientific information at a quality equivalent to or better than traditional approaches. The present case study examined whether bioactivity and associated potency measured in an in vitro screening assay for aromatase inhibition could be used together with an adverse outcome pathway (AOP) and mechanistically based computational models to predict previously uncharacterized in vivo effects. Model simulations were used to inform designs of 60-h and 10-21-day in vivo exposures of adult fathead minnows (Pimephales promelas) to three or four test concentrations of the in vitro aromatase inhibitor imazalil ranging from 0.12 to 260 µg/L water. Consistent with an AOP linking aromatase inhibition to reproductive impairment in fish, exposure to the fungicide resulted in significant reductions in ex vivo production of 17ß-estradiol (E2) by ovary tissue (≥165 µg imazalil/L), plasma E2 concentrations (≥74 µg imazalil/L), vitellogenin (Vtg) messenger RNA expression (≥165 µg imazalil/L), Vtg plasma concentrations (≥74 µg imazalil/L), uptake of Vtg into oocytes (≥260 µg imazalil/L), and overall reproductive output in terms of cumulative fecundity, number of spawning events, and eggs per spawning event (≥24 µg imazalil/L). Despite many potential sources of uncertainty in potency and efficacy estimates based on model simulations, observed magnitudes of apical effects were quite consistent with model predictions, and in vivo potency was within an order of magnitude of that predicted based on in vitro relative potency. Overall, our study suggests that NAMs and AOP-based approaches can support meaningful reduction and refinement of animal testing. Environ Toxicol Chem 2023;42:100-116. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Case Study in 21st Century Ecotoxicology: Using In Vitro Aromatase Inhibition Data to Predict Short-Term In Vivo Responses in Adult Female Fish.

The present study evaluated whether in vitro measures of aromatase inhibition as inputs into a quantitative adverse outcome pathway (qAOP) construct could effectively predict in vivo effects on 17ß-estradiol (E2) and vitellogenin (VTG) concentrations in female fathead minnows. Five chemicals identified as aromatase inhibitors in mammalian-based ToxCast assays were screened for their ability to inhibit fathead minnow aromatase in vitro. Female fathead minnows were then exposed to 3 of those chemicals: letrozole, epoxiconazole, and imazalil in concentration-response (5 concentrations plus control) for 24 h. Consistent with AOP-based expectations, all 3 chemicals caused significant reductions in plasma E2 and hepatic VTG transcription. Characteristic compensatory upregulation of aromatase and follicle-stimulating hormone receptor (fshr) transcripts in the ovary were observed for letrozole but not for the other 2 compounds. Considering the overall patterns of concentration-response and temporal concordance among endpoints, data from the in vivo experiments strengthen confidence in the qualitative relationships outlined by the AOP. Quantitatively, the qAOP model provided predictions that fell within the standard error of measured data for letrozole but not for imazalil and epoxiconazole. However, the inclusion of measured plasma concentrations of the test chemicals as inputs improved model predictions, with all predictions falling within the range of measured values. Results highlight both the utility and limitations of the qAOP and its potential use in 21st century ecotoxicology. Environ Toxicol Chem 2021;40:1155-1170. © 2020 SETAC. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

A method for CRISPR/Cas9 mutation of genes in fathead minnow (Pimephales promelas).

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 genome editing allows for the disruption or modification of genes in a multitude of model organisms. In the present study, we describe and employ the method for use in the fathead minnow (Pimephales promelas), in part, to assist in the development and validation of adverse outcome pathways (AOPs). The gene coding for an enzyme responsible for melanin production, tyrosinase (tyr), was the initial target chosen for development and assessment of the method since its disruption results in abnormal pigmentation, a phenotype obvious within 3-4 d after injection of fathead minnow embryos. Three tyrosinase-targeting guide strands were generated using the fathead minnow sequence in tandem with the CRISPOR guide strand selection tool. The strands targeted two areas: one stretch of sequence in a conserved region that demonstrated homology to EGF-like or laminin-like domains as determined by Protein Basic Local Alignment Search Tool in concert with the Conserved Domain Database, and a second area in the N-terminal region of the tyrosinase domain. To generate one cell embryos, in vitro fertilization was performed, allowing for microinjection of hundreds of developmentally-synchronized embryos with Cas9 proteins complexed to each of the three guide strands. Altered retinal pigmentation was observed in a portion of the tyr guide strand injected population within 3 d post fertilization (dpf). By 14 dpf, fish without skin and swim bladder pigmentation were observed. Among the three guide strands injected, the guide targeting the EGF/laminin-like domain was most effective in generating mutants. CRISPR greatly advances our ability to directly investigate gene function in fathead minnow, allowing for advanced approaches to AOP validation and development.

Adverse Outcome Pathway Network-Based Assessment of the Interactive Effects of an Androgen Receptor Agonist and an Aromatase Inhibitor on Fish Endocrine Function.

Predictive approaches to assessing the toxicity of contaminant mixtures have been largely limited to chemicals that exert effects through the same biological molecular initiating event. However, by understanding specific pathways through which chemicals exert effects, it may be possible to identify shared "downstream" nodes as the basis for forecasting interactive effects of chemicals with different molecular initiating events. Adverse outcome pathway (AOP) networks conceptually support this type of analysis. We assessed the utility of a simple AOP network for predicting the effects of mixtures of an aromatase inhibitor (fadrozole) and an androgen receptor agonist (17ß-trenbolone) on aspects of reproductive endocrine function in female fathead minnows. The fish were exposed to multiple concentrations of fadrozole and 17ß-trenbolone individually or in combination for 48 or 96 h. Effects on 2 shared nodes in the AOP network, plasma 17ß-estradiol (E2) concentration and vitellogenin (VTG) production (measured as hepatic vtg transcripts) responded as anticipated to fadrozole alone but were minimally impacted by 17ß-trenbolone alone. Overall, there were indications that 17ß-trenbolone enhanced decreases in E2 and vtg in fadrozole-exposed fish, as anticipated, but the results often were not statistically significant. Failure to consistently observe hypothesized interactions between fadrozole and 17ß-trenbolone could be due to several factors, including lack of impact of 17ß-trenbolone, inherent biological variability in the endpoints assessed, and/or an incomplete understanding of interactions (including feedback) between different pathways within the hypothalamic-pituitary-gonadal axis. Environ Toxicol Chem 2020;39:913-922. © 2020 SETAC.

First-generation annotations for the fathead minnow (Pimephales promelas) genome.

Ab initio gene prediction and evidence alignment were used to produce the first annotations for the fathead minnow (Pimephales promelas) genome. We also describe a genome browser, hosted by the Society of Environmental Toxicology and Chemistry, that provides simplified access to the annotation data in context with the genomic sequence. The present study extends the utility of the fathead minnow genome and supports the continued development of this species as a model organism for predictive toxicology. Environ Toxicol Chem 2017;36:3436-3442. Published 2017 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Editor's Highlight: Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS): A Web-Based Tool for Addressing the Challenges of Cross-Species Extrapolation of Chemical Toxicity.

Conservation of a molecular target across species can be used as a line-of-evidence to predict the likelihood of chemical susceptibility. The web-based Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS; https://seqapass.epa.gov/seqapass/) application was developed to simplify, streamline, and quantitatively assess protein sequence/structural similarity across taxonomic groups as a means to predict relative intrinsic susceptibility. The intent of the tool is to allow for evaluation of any potential protein target while remaining amenable to variable degrees of protein characterization, in the context of available information about the chemical/protein interaction and the molecular target itself. To accommodate this flexibility in the analysis, 3 levels of evaluation were developed. The first level of the SeqAPASS analysis compares primary amino acid sequences to a query sequence, calculating a metric for sequence similarity (including detection of orthologs); the second level evaluates sequence similarity within selected functional domains (eg, ligand-binding domain); and the third level of analysis compares individual amino acid residue positions of importance for protein conformation and/or interaction with the chemical upon binding. Each level of the SeqAPASS analysis provides additional evidence to apply toward rapid, screening-level assessments of probable cross species susceptibility. Such analyses can support prioritization of chemicals for further evaluation, selection of appropriate species for testing, extrapolation of empirical toxicity data, and/or assessment of the cross-species relevance of adverse outcome pathways. Three case studies are described herein to demonstrate application of the SeqAPASS tool: the first 2 focused on predictions of pollinator susceptibility to molt-accelerating compounds and neonicotinoid insecticides, and the third on evaluation of cross-species susceptibility to strobilurin fungicides. These analyses illustrate challenges in species extrapolation and demonstrate the broad utility of SeqAPASS for risk-based decision making and research.

Impaired anterior swim bladder inflation following exposure to the thyroid peroxidase inhibitor 2-mercaptobenzothiazole part I: Fathead minnow.

In the present study, a hypothesized adverse outcome pathway linking inhibition of thyroid peroxidase (TPO) activity to impaired swim bladder inflation was investigated in two experiments in which fathead minnows (Pimephales promelas) were exposed to 2-mercaptobenzothiazole (MBT). Continuous exposure to 1mg MBT/L for up to 22 days had no effect on inflation of the posterior chamber of the swim bladder, which typically inflates around 6 days post fertilization (dpf), a period during which maternally-derived thyroid hormone is presumed to be present. In contrast, inflation of the anterior swim bladder, which occurs around 14dpf, was impacted. Specifically, at 14dpf, approximately 50% of fish exposed to 1mg MBT/L did not have an inflated anterior swim bladder. In fish exposed to MBT through 21 or 22dpf, the anterior swim bladder was able to inflate, but the ratio of the anterior/posterior chamber length was significantly reduced compared to controls. Both abundance of thyroid peroxidase mRNA and thyroid follicle histology suggest that fathead minnows mounted a compensatory response to the presumed inhibition of TPO activity by MBT. Time-course characterization showed that fish exposed to MBT for at least 4 days prior to normal anterior swim bladder inflation had significant reductions in anterior swim bladder size, relative to the posterior chamber, compared to controls. These results, along with similar results observed in zebrafish (see part II, this issue) are consistent with the hypothesis that thyroid hormone signaling plays a significant role in mediating anterior swim bladder inflation and development in cyprinids, and that role can be disrupted by exposure to thyroid hormone synthesis inhibitors. Nonetheless, possible thyroid-independent actions of MBT on anterior swim bladder inflation cannot be ruled out based on the present results. Overall, although anterior swim bladder inflation has not been directly linked to survival as posterior swim bladder inflation has, potential links to adverse ecological outcomes are plausible given involvement of the anterior chamber in sound production and detection.