Folate Intake Alters Mutation Frequency and Profiles in a Tissue- and Dose-Specific Manner in MutaMouse Male Mice.

BACKGROUND: While cancer is common, its incidence varies widely by tissue. These differences are attributable to variable risk factors, such as environmental exposure, genetic inheritance, and lifetime number of stem cell divisions in a tissue. Folate deficiency is generally associated with increased risk for colorectal cancer (CRC) and acute lymphocytic leukemia (ALL). Conversely, high folic acid (FA) intake has also been associated with higher CRC risk. OBJECTIVE: Our objective was to compare the effect of folate intake on mutant frequency (MF) and types of mutations in the colon and bone marrow of mice. METHODS: Five-week-old MutaMouse male mice were fed a deficient (0 mg FA/kg), control (2 mg FA/kg), or supplemented (8 mg FA/kg) diet for 20 wk. Tissue MF was assessed using the lacZ mutant assay and comparisons made by 2-factor ANOVA. LacZ mutant plaques were sequenced using next-generation sequencing, and diet-specific mutation profiles within each tissue were compared by Fisher's exact test. RESULTS: In the colon, the MF was 1.5-fold and 1.3-fold higher in mice fed the supplemented diet compared with mice fed the control (P = 0.001) and deficient (P = 0.008) diets, respectively. This contrasted with the bone marrow MF in the same mice where the MF was 1.7-fold and 1.6-fold higher in mice fed the deficient diet compared with mice fed the control (P = 0.02) and supplemented (P = 0.03) diets, respectively. Mutation profiles and signatures (mutation context) were tissue-specific. CONCLUSIONS: Our data indicate that dietary folate intake affects mutagenesis in a tissue- and dose-specific manner in mice. Mutation profiles were generally tissue- but not dose-specific, suggesting that altered cellular folate status appears to interact with endogenous mutagenic mechanisms in each tissue to create a permissive context in which specific mutation types accumulate. These data illuminate potential mechanisms underpinning differences in observed associations between folate intake/status and cancer.

Transgenerational hypocortisolism and behavioral disruption are induced by the antidepressant fluoxetine in male zebrafish Danio rerio.

The global prevalence of depression is high during childbearing. Due to the associated risks to the mother and baby, the selective serotonin reuptake inhibitor fluoxetine (FLX) is often the first line of treatment. Given that FLX readily crosses the placenta, a fetus may be susceptible to the disruptive effects of FLX during this highly plastic stage of development. Here, we demonstrate that a 6-day FLX exposure to a fetus-relevant concentration at a critical developmental stage suppresses cortisol levels in the adult zebrafish (F0). This effect persists for three consecutive generations in the unexposed descendants (F1 to F3) without diminution and is more pronounced in males. We also show that the in vivo cortisol response of the interrenal (fish "adrenal") to an i.p. injection of adrenocorticotropic hormone was also reduced in the males from the F0 and F3 FLX lineages. Transcriptomic profiling of the whole kidney containing the interrenal cells revealed that early FLX exposure significantly modified numerous pathways closely associated with cortisol synthesis in the male adults from the F0 and F3 generations. We also show that the low cortisol levels are linked to significantly reduced exploratory behaviors in adult males from the F0 to F2 FLX lineages. This may be a cause for concern given the high prescription rates of FLX to pregnant women and the potential long-term negative impacts on humans exposed to these therapeutic drugs.

Transcriptional profiling of male CD-1 mouse lungs and Harderian glands supports the involvement of calcium signaling in acrylamide-induced tumors.

Acrylamide (AA) exposure causes increased incidence of forestomach, lung, and Harderian gland tumors in male mice. One hypothesized mode of action (MOA) for AA-carcinogenicity includes genotoxicity/mutagenicity as a key event, possibly resulting from AA metabolism to the direct genotoxic metabolite glycidamide. Alternatively, altered calcium signaling (CS) has been proposed as a central key event in the MOA. To examine the plausibility of these proposed MOAs, RNA-sequencing was performed on tumor target tissues: Harderian glands (the most sensitive tumor target tissue in the rodent 2-year cancer bioassay) and lungs of AA-exposed male CD-1 mice. Animals were exposed to 0.0, 1.5, 3.0, 6.0, 12.0, or 24.0 mg AA/kg bw-day in drinking water for 5, 15, or 31 days. We observed a pronounced effect on genes involved in CS and cytoskeletal processes in both tissues, but no evidence supporting a genotoxic MOA. Benchmark dose modeling suggests transcriptional points of departure (PODs) of 0.54 and 2.21 mg/kg bw-day for the Harderian glands and lungs, respectively. These are concordant with PODs of 0.17 and 1.27 mg/kg bw-day derived from the cancer bioassay data for these tissues in male mice, respectively. Overall, this study supports the involvement of CS in AA-induced mouse carcinogenicity, which is consistent with a recently proposed CS-based MOA in rat thyroid, and with other published reports of aberrant CS in malignant tumors in rodents and humans.

Genetic Toxicity of Complex Mixtures of Polycyclic Aromatic Hydrocarbons: Evaluating Dose-Additivity in a Transgenic Mouse Model.

This study evaluates the risk assessment approach currently employed for polycyclic aromatic hydrocarbon (PAH)-contaminated media, wherein carcinogenic hazards are evaluated using a dose-addition model that employs potency equivalency factors (PEFs) for targeted carcinogenic PAHs. Here, MutaMouse mice were subchronically exposed to PAH mixtures (p.o.), and mutagenic potency (MP) values were determined for five tissues. Predicted dose-additive mixture MPs were generated by summing the products of the concentrations and MPs of the individual targeted PAHs; values were compared to the experimental MPs of the mixtures to evaluate dose-additivity. Additionally, the PEF-determined BaP-equivalent concentrations were compared to those determined using a bioassay-derived method (BDM) (i.e., an additivity-independent approach). In bone marrow, mixture mutagenicity was less than dose-additive and the PEF-method provided higher estimates of BaP-equivalents than the BDM. Conversely, mixture mutagenicity in site-of-contact tissues (e.g., small intestine) was generally more than dose-additive and the PEF-method provided lower estimates of BaP-equivalents than the BDM. Overall, this study demonstrates that dose-additive predictions of mixture mutagenic potency based on the concentrations and potencies of a small number of targeted PAHs results in values that are surprisingly close to those determined experimentally, providing support for the dose-additive assumption employed for human health risk assessment of PAH mixtures.

Characterizing Benzo[a]pyrene-induced lacZ mutation spectrum in transgenic mice using next-generation sequencing.

BACKGROUND: The transgenic rodent mutation reporter assay provides an efficient approach to identify mutagenic agents in vivo. A major advantage of this assay is that mutant reporter transgenes can be sequenced to provide information on the mode of action of a mutagen and to identify clonally expanded mutations. However, conventional DNA sequence analysis is laborious and expensive for long transgenes, such as lacZ (3096 bp), and is not normally implemented in routine screening. METHODS: We developed a high-throughput next-generation sequencing (NGS) approach to simultaneously sequence large numbers of barcoded mutant lacZ transgenes from different animals. We collected 3872 mutants derived from the bone marrow DNA of six Muta™Mouse males exposed to the well-established mutagen benzo[a]pyrene (BaP) and six solvent-exposed controls. Mutants within animal samples were pooled, barcoded, and then sequenced using NGS. RESULTS: We identified 1652 mutant sequences from 1006 independent mutations that underwent clonal expansion. This deep sequencing analysis of mutation spectrum demonstrated that BaP causes primarily guanine transversions (e.g. G:C → T:A), which is highly consistent with previous studies employing Sanger sequencing. Furthermore, we identified novel mutational hotspots in the lacZ transgene that were previously uncharacterized by Sanger sequencing. Deep sequencing also allowed for an unprecedented ability to correct for clonal expansion events, improving the sensitivity of the mutation reporter assay by 50 %. CONCLUSION: These results demonstrate that the high-throughput nature and reduced costs offered by NGS provide a sensitive and fast approach for elucidating and comparing mutagenic mechanisms of various agents among tissues and enabling improved evaluation of genotoxins.

Identification of thyroid hormone receptor binding sites in developing mouse cerebellum.

BACKGROUND: Thyroid hormones play an essential role in early vertebrate development as well as other key processes. One of its modes of action is to bind to the thyroid hormone receptor (TR) which, in turn, binds to thyroid response elements (TREs) in promoter regions of target genes. The sequence motif for TREs remains largely undefined as does the precise chromosomal location of the TR binding sites. A chromatin immunoprecipitation on microarray (ChIP-chip) experiment was conducted using mouse cerebellum post natal day (PND) 4 and PND15 for the thyroid hormone receptor (TR) beta 1 to map its binding sites on over 5000 gene promoter regions. We have performed a detailed computational analysis of these data. RESULTS: By analysing a recent spike-in study, the optimal normalization and peak identification approaches were determined for our dataset. Application of these techniques led to the identification of 211 ChIP-chip peaks enriched for TR binding in cerebellum samples. ChIP-PCR validation of 25 peaks led to the identification of 16 true positive TREs. Following a detailed literature review to identify all known mouse TREs, a position weight matrix (PWM) was created representing the classic TRE sequence motif. Various classes of promoter regions were investigated for the presence of this PWM, including permuted sequences, randomly selected promoter sequences, and genes known to be regulated by TH. We found that while the occurrence of the TRE motif is strongly correlated with gene regulation by TH for some genes, other TH-regulated genes do not exhibit an increased density of TRE half-site motifs. Furthermore, we demonstrate that an increase in the rate of occurrence of the half-site motifs does not always indicate the specific location of the TRE within the promoter region. To account for the fact that TR often operates as a dimer, we introduce a novel dual-threshold PWM scanning approach for identifying TREs with a true positive rate of 0.73 and a false positive rate of 0.2. Application of this approach to ChIP-chip peak regions revealed the presence of 85 putative TREs suitable for further in vitro validation. CONCLUSIONS: This study further elucidates TRß gene regulation in mouse cerebellum, with 211 promoter regions identified to bind to TR. While we have identified 85 putative TREs within these regions, future work will study other mechanisms of action that may mediate the remaining observed TR-binding activity.

Brief Developmental Exposure to Fluoxetine Causes Life-Long Alteration of the Brain Transcriptome in Zebrafish.

Fluoxetine (FLX) and other selective serotonin reuptake inhibitors are widely used to treat depressive disorders during pregnancy. Early-life exposure to FLX is known to disrupt the normal function of the stress axis in humans, rodents, and teleosts. We used a zebrafish line with a cortisol-inducible fluorescent transgene to study the effects of developmental daily exposure to FLX (54 µg/L) on the transcriptomic profile of brain tissues in exposed larvae and later as 6-month-old adults. High throughput RNA sequencing was conducted on brain tissues in unstressed and stressed conditions. Long-lasting effects of FLX were observed in telencephalon (Tel) and hypothalamus (Hyp) of adult zebrafish with 1927 and 5055 genes significantly (≥1.2 fold-change, false-discovery p-value < 0.05) dysregulated in unstressed condition, respectively. Similar findings were observed in Hyp with 1245 and 723 genes being significantly dysregulated in stressed adults, respectively. Differentially expressed genes converted to Homo sapiens orthologues were used for Ingenuity Pathway Analysis. The results showed alteration of pathways involved in neuroendocrine signaling, cholesterol metabolism and synaptogenesis. Enriched networks included lipid metabolism, molecular transport, and nervous system development. Analysis of putative upstream transcription regulators showed potential dysregulation of clocka and nr3c1 which control circadian rhythm, stress response, cholesterol metabolism and histone modifications. Several genes involved in epigenetic regulation were also affected by FLX, including dnmt3a, adarb1, adarb2, hdac4, hdac5, hdac8, and atf2. We report life-long disruptive effects of FLX on pathways associated with neuroendocrine signaling, stress response and the circadian rhythm, and all of which are implicated in the development of depressive disorders in humans. Our results raise concern for the persistent endocrine-disrupting potential of brief antidepressant exposure during embryonic development.

Integrated Genotoxicity Testing of three anti-infective drugs using the TGx-DDI transcriptomic biomarker and high-throughput CometChip® assay in TK6 cells.

Genotoxicity testing relies on the detection of gene mutations and chromosome damage and has been used in the genetic safety assessment of drugs and chemicals for decades. However, the results of standard genotoxicity tests are often difficult to interpret due to lack of mode of action information. The TGx-DDI transcriptomic biomarker provides mechanistic information on the DNA damage-inducing (DDI) capability of chemicals to aid in the interpretation of positive in vitro genotoxicity data. The CometChip® assay was developed to assess DNA strand breaks in a higher-throughput format. We paired the TGx-DDI biomarker with the CometChip® assay in TK6 cells to evaluate three model agents: nitrofurantoin (NIT), metronidazole (MTZ), and novobiocin (NOV). TGx-DDI was analyzed by two independent labs and technologies (nCounter® and TempO-Seq®). Although these anti-infective drugs are, or have been, used in human and/or veterinary medicine, the standard genotoxicity testing battery showed significant genetic safety findings. Specifically, NIT is a mutagen and causes chromosome damage, and MTZ and NOV cause chromosome damage in conventional in vitro tests. Herein, the TGx-DDI biomarker classified NIT and MTZ as non-DDI at all concentrations tested, suggesting that NIT's mutagenic activity is bacterial specific and that the observed chromosome damage by MTZ might be a consequence of in vitro test conditions. In contrast, NOV was classified as DDI at the second highest concentration tested, which is in line with the fact that NOV is a bacterial DNA-gyrase inhibitor that also affects topoisomerase II at high concentrations. The lack of DNA damage for NIT and MTZ was confirmed by the CometChip® results, which were negative for all three drugs except at overtly cytotoxic concentrations. This case study demonstrates the utility of combining the TGx-DDI biomarker and CometChip® to resolve conflicting genotoxicity data and provides further validation to support the reproducibility of the biomarker.

Thyroid hormone-regulated gene expression in juvenile mouse liver: identification of thyroid response elements using microarray profiling and in silico analyses.

BACKGROUND: Disruption of thyroid hormone signalling can alter growth, development and energy metabolism. Thyroid hormones exert their effects through interactions with thyroid receptors that directly bind thyroid response elements and can alter transcriptional activity of target genes. The effects of short-term thyroid hormone perturbation on hepatic mRNA transcription in juvenile mice were evaluated, with the goal of identifying genes containing active thyroid response elements. Thyroid hormone disruption was induced from postnatal day 12 to 15 by adding goitrogens to dams' drinking water (hypothyroid). A subgroup of thyroid hormone-disrupted pups received intraperitoneal injections of replacement thyroid hormones four hours prior to sacrifice (replacement). An additional group received only thyroid hormones four hours prior to sacrifice (hyperthyroid). Hepatic mRNA was extracted and hybridized to Agilent mouse microarrays. RESULTS: Transcriptional profiling enabled the identification of 28 genes that appeared to be under direct thyroid hormone-regulation. The regulatory regions of the genome adjacent to these genes were examined for half-site sequences that resemble known thyroid response elements. A bioinformatics search identified 33 thyroid response elements in the promoter regions of 13 different genes thought to be directly regulated by thyroid hormones. Thyroid response elements found in the promoter regions of Tor1a, 2310003H01Rik, Hect3d and Slc25a45 were further validated by confirming that the thyroid receptor is associated with these sequences in vivo and that it can bind directly to these sequences in vitro. Three different arrangements of thyroid response elements were identified. Some of these thyroid response elements were located far up-stream (> 7 kb) of the transcription start site of the regulated gene. CONCLUSIONS: Transcriptional profiling of thyroid hormone disrupted animals coupled with a novel bioinformatics search revealed new thyroid response elements associated with genes previously unknown to be responsive to thyroid hormone. The work provides insight into thyroid response element sequence motif characteristics.

Transcriptomic pathway and benchmark dose analysis of Bisphenol A, Bisphenol S, Bisphenol F, and 3,3',5,5'-Tetrabromobisphenol A in H9 human embryonic stem cells.

Bisphenol A (BPA) is a chemical used in the manufacturing of plastics to which human exposure is ubiquitous. Numerous studies have linked BPA exposure to many adverse health outcomes prompting the replacement of BPA with various analogues including bisphenol-F (BPF) and bisphenol S (BPS). Other bisphenols are used in various consumer applications, such as 3,3',5,5'-Tetrabromobisphenol A (TBBPA), which is used as a flame retardant. Few studies to date have examined the effects of BPA and its analogues in stem cells to explore potential developmental impacts. Here we used transcriptomics to investigate similarities and differences of BPA and three of its analogues in the estrogen receptor negative, human embryonic stem cell line H9 (WA09). H9 cells were exposed to increasing concentrations of the bisphenols and analyzed using RNA-sequencing. Our data indicate that BPA, BPF, and BPS have similar potencies in inducing transcriptional changes and perturb many of the same pathways. TBBPA, the least structurally similar bisphenol of the group, exhibited much lower potency. All bisphenols robustly impacted gene expression in these cells, albeit at concentrations well above those observed in estrogen-positive cells. Overall, we provide a foundational data set against which to explore the transcriptional similarities of other bisphenols in embryonic stem cells, which may be used to assess the suitability of chemical grouping for read-across and for preliminary potency evaluation.

Potency Ranking of Per- and Polyfluoroalkyl Substances Using High-Throughput Transcriptomic Analysis of Human Liver Spheroids.

Per- and polyfluoroalkyl substances (PFAS) are some of the most prominent organic contaminants in human blood. Although the toxicological implications of human exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) are well established, data on lesser-understood PFAS are limited. New approach methodologies (NAMs) that apply bioinformatic tools to high-throughput data are being increasingly considered to inform risk assessment for data-poor chemicals. The aim of this study was to compare the potencies (ie, benchmark concentrations: BMCs) of PFAS in primary human liver microtissues (3D spheroids) using high-throughput transcriptional profiling. Gene expression changes were measured using TempO-seq, a templated, multiplexed RNA-sequencing platform. Spheroids were exposed for 1 or 10 days to increasing concentrations of 23 PFAS in 3 subgroups: carboxylates (PFCAs), sulfonates (PFSAs), and fluorotelomers and sulfonamides. PFCAs and PFSAs exhibited trends toward increased transcriptional potency with carbon chain-length. Specifically, longer-chain compounds (7-10 carbons) were more likely to induce changes in gene expression and have lower transcriptional BMCs. The combined high-throughput transcriptomic and bioinformatic analyses support the capability of NAMs to efficiently assess the effects of PFAS in liver microtissues. The data enable potency ranking of PFAS for human liver cell spheroid cytotoxicity and transcriptional changes, and assessment of in vitro transcriptomic points of departure. These data improve our understanding of the possible health effects of PFAS and will be used to inform read-across for human health risk assessment.

A Modern Genotoxicity Testing Paradigm: Integration of the High-Throughput CometChip® and the TGx-DDI Transcriptomic Biomarker in Human HepaRG™ Cell Cultures.

Higher-throughput, mode-of-action-based assays provide a valuable approach to expedite chemical evaluation for human health risk assessment. In this study, we combined the high-throughput alkaline DNA damage-sensing CometChip® assay with the TGx-DDI transcriptomic biomarker (DDI = DNA damage-inducing) using high-throughput TempO-Seq®, as an integrated genotoxicity testing approach. We used metabolically competent differentiated human HepaRG™ cell cultures to enable the identification of chemicals that require bioactivation to cause genotoxicity. We studied 12 chemicals (nine DDI, three non-DDI) in increasing concentrations to measure and classify chemicals based on their ability to damage DNA. The CometChip® classified 10/12 test chemicals correctly, missing a positive DDI call for aflatoxin B1 and propyl gallate. The poor detection of aflatoxin B1 adducts is consistent with the insensitivity of the standard alkaline comet assay to bulky lesions (a shortcoming that can be overcome by trapping repair intermediates). The TGx-DDI biomarker accurately classified 10/12 agents. TGx-DDI correctly identified aflatoxin B1 as DDI, demonstrating efficacy for combined used of these complementary methodologies. Zidovudine, a known DDI chemical, was misclassified as it inhibits transcription, which prevents measurable changes in gene expression. Eugenol, a non-DDI chemical known to render misleading positive results at high concentrations, was classified as DDI at the highest concentration tested. When combined, the CometChip® assay and the TGx-DDI biomarker were 100% accurate in identifying chemicals that induce DNA damage. Quantitative benchmark concentration (BMC) modeling was applied to evaluate chemical potencies for both assays. The BMCs for the CometChip® assay and the TGx-DDI biomarker were highly concordant (within 4-fold) and resulted in identical potency rankings. These results demonstrate that these two assays can be integrated for efficient identification and potency ranking of DNA damaging agents in HepaRG™ cell cultures.

High-Throughput Transcriptomic Analysis of Human Primary Hepatocyte Spheroids Exposed to Per- and Polyfluoroalkyl Substances as a Platform for Relative Potency Characterization.

Per- and poly-fluoroalkyl substances (PFAS) are widely found in the environment because of their extensive use and persistence. Although several PFAS are well studied, most lack toxicity data to inform human health hazard and risk assessment. This study focused on 4 model PFAS: perfluorooctanoic acid (PFOA; 8 carbon), perfluorobutane sulfonate (PFBS; 4 carbon), perfluorooctane sulfonate (PFOS; 8 carbon), and perfluorodecane sulfonate (PFDS; 10 carbon). Human primary liver cell spheroids (pooled from 10 donors) were exposed to 10 concentrations of each PFAS and analyzed at 4 time points. The approach aimed to: (1) identify gene expression changes mediated by the PFAS, (2) identify similarities in biological responses, (3) compare PFAS potency through benchmark concentration analysis, and (4) derive bioactivity exposure ratios (ratio of the concentration at which biological responses occur, relative to daily human exposure). All PFAS induced transcriptional changes in cholesterol biosynthesis and lipid metabolism pathways, and predicted PPARα activation. PFOS exhibited the most transcriptional activity and had a highly similar gene expression profile to PFDS. PFBS induced the least transcriptional changes and the highest benchmark concentration (ie, was the least potent). The data indicate that these PFAS may have common molecular targets and toxicities, but that PFOS and PFDS are the most similar. The transcriptomic bioactivity exposure ratios derived here for PFOA and PFOS were comparable to those derived using rodent apical endpoints in risk assessments. These data provide a baseline level of toxicity for comparison with other known PFAS using this testing strategy.

Prolonged Low-Dose Dioxin Exposure Impairs Metabolic Adaptability to High-Fat Diet Feeding in Female but Not Male Mice.

CONTEXT: Human studies consistently show an association between exposure to persistent organic pollutants, including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD, aka "dioxin"), and increased diabetes risk. We previously showed that a single high-dose TCDD exposure (20 µg/kg) decreased plasma insulin levels in male and female mice in vivo, but effects on glucose homeostasis were sex-dependent. OBJECTIVE: The current study assessed whether prolonged exposure to a physiologically relevant low-dose of TCDD impacts glucose homeostasis and/or the islet phenotype in a sex-dependent manner in chow-fed or high-fat diet (HFD)-fed mice. METHODS: Male and female mice were exposed to 20 ng/kg/d TCDD 2×/week for 12 weeks and simultaneously fed standard chow or a 45% HFD. Glucose homeostasis was assessed by glucose and insulin tolerance tests, and glucose-induced plasma insulin levels were measured in vivo. Histological analysis was performed on pancreas from male and female mice, and islets were isolated from females for TempO-Seq transcriptomic analysis. RESULTS: Low-dose TCDD exposure did not lead to adverse metabolic consequences in chow-fed male or female mice, or in HFD-fed males. However, TCDD accelerated the onset of HFD-induced hyperglycemia and impaired glucose-induced plasma insulin levels in females. TCDD caused a modest increase in islet area in males but reduced the percent beta cell area within islets in females. TempO-Seq analysis suggested abnormal changes to endocrine and metabolic pathways in female TCDDHFD islets. CONCLUSION: Our data suggest that prolonged low-dose TCDD exposure has minimal effects on glucose homeostasis and islet morphology in chow-fed male and female mice but promotes maladaptive metabolic responses in HFD-fed females.

Depot-Specific Analysis of Human Adipose Cells and Their Responses to Bisphenol S.

Exposure to endocrine-disrupting chemicals (EDCs) is associated with adverse health outcomes including obesity and diabetes. Obesity, and more specifically visceral obesity, is correlated with metabolic disease. The adipose tissue is an endocrine organ and a potential target for many environmental pollutants including bisphenols. The subcutaneous (Sc) and the omental (Om, visceral) depots are composed of mature adipocytes and residing progenitors, which may be different between the depots and may be EDCs targets. Bisphenol A (BPA) is a suspected metabolic disruptor, and is being replaced with structurally similar compounds such as bisphenol S (BPS). Like BPA, BPS induces adipogenesis in murine and primary human Sc preadipocytes. However, the effect of BPS on Om preadipocytes is not known. In this study, we show that human primary progenitors from Om depots have a distinct transcriptomic signature as compared to progenitors derived from donor-matched Sc depots. Furthermore, we show that BPS increases adipogenesis both of Om and Sc preadipocytes and can mimic the action of glucocorticoids or peroxisome proliferator-activated receptor γ (PPARγ) agonists. We also show that BPS treatment, at 0.1 µM and 25 µM, modifies the adipokine profiles both of Om- and Sc-derived adipocytes in a depot-specific manner. Taken together our data show distinct gene expression profiles in the Om vs Sc progenitors and similar responses to the BPA analogue, BPS.

Flow cytometric micronucleus assay and TGx-DDI transcriptomic biomarker analysis of ten genotoxic and non-genotoxic chemicals in human HepaRG™ cells.

BACKGROUND: Modern testing paradigms seek to apply human-relevant cell culture models and integrate data from multiple test systems to accurately inform potential hazards and modes of action for chemical toxicology. In genetic toxicology, the use of metabolically competent human hepatocyte cell culture models provides clear advantages over other more commonly used cell lines that require the use of external metabolic activation systems, such as rat liver S9. HepaRG™ cells are metabolically competent cells that express Phase I and II metabolic enzymes and differentiate into mature hepatocyte-like cells, making them ideal for toxicity testing. We assessed the performance of the flow cytometry in vitro micronucleus (MN) test and the TGx-DDI transcriptomic biomarker to detect DNA damage-inducing (DDI) chemicals in human HepaRG™ cells after a 3-day repeat exposure. The biomarker, developed for use in human TK6 cells, is a panel of 64 genes that accurately classifies chemicals as DDI or non-DDI. Herein, the TGx-DDI biomarker was analyzed by Ion AmpliSeq whole transcriptome sequencing to assess its classification accuracy using this more modern gene expression technology as a secondary objective. METHODS: HepaRG™ cells were exposed to increasing concentrations of 10 test chemicals (six genotoxic chemicals, including one aneugen, and four non-genotoxic chemicals). Cytotoxicity and genotoxicity were measured using the In Vitro MicroFlow® kit, which was run in parallel with the TGx-DDI biomarker. RESULTS: A concentration-related decrease in relative survival and a concomitant increase in MN frequency were observed for genotoxic chemicals in HepaRG™ cells. All five DDI and five non-DDI agents were correctly classified (as genotoxic/non-genotoxic and DDI/non-DDI) by pairing the test methods. The aneugenic agent (colchicine) yielded the expected positive result in the MN test and negative (non-DDI) result by TGx-DDI. CONCLUSIONS: This next generation genotoxicity testing strategy is aligned with the paradigm shift occurring in the field of genetic toxicology. It provides mechanistic insight in a human-relevant cell-model, paired with measurement of a conventional endpoint, to inform the potential for adverse health effects. This work provides support for combining these assays in an integrated test strategy for accurate, higher throughput genetic toxicology testing in this metabolically competent human progenitor cell line.

The use of in vitro transcriptional data to identify thresholds of effects in a human lens epithelial cell-line exposed to ionizing radiation.

PURPOSE: The International Commission on Radiological Protection (ICRP) recently recommended reducing the occupational equivalent dose limit for the lens of the eye. Based primarily on a review of epidemiological data, the absorbed dose threshold is now considered to be 0.5 Gy independent of dose-rate and severity of opacification, reduced from the previous threshold of 2 Gy. However, direct mechanistic evidence to support an understanding of the underlying molecular mechanisms of damage is still lacking. To this end, we explored the effects of a broad dose-range of ionizing radiation exposure on gene expression changes in a human lens epithelial (HLE) cell-line in order to better understand the shape of the dose-response relationship and identify transcriptional thresholds of effects. METHODS: HLE cells were exposed to doses of 0, 0.01, 0.05, 0.25, 0.5, 2, and 5 Gy of X-ray radiation at two dose rates (1.62 cGy/min and 38.2 cGy/min). Cell culture lysates were collected 20 h post-exposure and analyzed using whole-genome RNA-sequencing. Pathways and dose-thresholds of biological effects were identified using benchmark dose (BMD) modeling. RESULTS: Transcriptional responses were minimal at doses less than 2 Gy. At higher doses, there were a significant number of differentially expressed genes (DEGs) (p≤.05, fold change≥|1.5|) at both dose rates, with 1308 DEGs for the low dose rate (LDR) and 840 DEGs for the high dose rate (HDR) exposure. Dose-response modeling showed that a number of genes exhibited non-linear bi-phasic responses, which was verified by digital droplet PCR. BMD analysis showed the majority of the pathways responded at BMD median values in the dose range of 1.5-2.5 Gy, with the lowest BMD median value being 0.6 Gy for the HDR exposure. The minimum pathway BMD median value for LDR exposure, however, was 2.5 Gy. Although the LDR and HDR exposures shared pathways involved in extracellular matrix reorganization and collagen production with BMD median value of 2.9 Gy, HDR exposures were more effective in activating pathways associated with DNA damage response, apoptosis, and cell cycling relative to LDR exposure. CONCLUSIONS: Overall, the results suggest that radiation induces complex non-linear transcriptional dose-response relationships that are dose-rate dependent. Pathways shared between the two dose rates may be important contributors to radiation-induced cataractogenesis. BMD analysis suggests that the majority of pathways are activated above 0.6 Gy, which supports current ICRP identified dose thresholds for deterministic effects to the lens of the eye of 0.5 Gy.

Hepatic transcriptional dose-response analysis of male and female Fischer rats exposed to hexabromocyclododecane.

Hexabromocyclododecane (HBCD) is a brominated flame retardant found in the environment and human tissues. The toxicological effects of HBCD exposure are not clearly understood. We employed whole-genome RNA-sequencing on liver samples from male and female Fischer rats exposed to 0, 250, 1250, and 5000 mg technical mixture of HBCD/kg diet for 28 days to gain further insight into HBCD toxicity. HBCD altered 428 and 250 gene transcripts in males and females, respectively, which were involved in metabolism of xenobiotics, oxidative stress, immune response, metabolism of glucose and lipids, circadian regulation, cell cycle, fibrotic activity, and hormonal balance. Signature analysis supported that HBCD operates through the constitutive androstane and pregnane X receptors. The median transcriptomic benchmark dose (BMD) for the lowest statistically significant pathway was within 1.5-fold of the BMD for increased liver weight, while the BMD for the lowest pathway with at least three modeled genes (minimum 5% of pathway) was similar to the lowest apical endpoint BMD. The results show how transcriptional analyses can inform mechanisms underlying chemical toxicity and the doses at which potentially adverse effects occur. This experiment is part of a larger study exploring the use of toxicogenomics and high-throughput screening for human health risk assessment.

Paternal exposure to benzo(a)pyrene induces genome-wide mutations in mouse offspring.

Understanding the effects of environmental exposures on germline mutation rates has been a decades-long pursuit in genetics. We used next-generation sequencing and comparative genomic hybridization arrays to investigate genome-wide mutations in the offspring of male mice exposed to benzo(a)pyrene (BaP), a common environmental pollutant. We demonstrate that offspring developing from sperm exposed during the mitotic or post-mitotic phases of spermatogenesis have significantly more de novo single nucleotide variants (1.8-fold; P < 0.01) than controls. Both phases of spermatogenesis are susceptible to the induction of heritable mutations, although mutations arising from post-fertilization events are more common after post-mitotic exposure. In addition, the mutation spectra in sperm and offspring of BaP-exposed males are consistent. Finally, we report a significant increase in transmitted copy number duplications (P = 0.001) in BaP-exposed sires. Our study demonstrates that germ cell mutagen exposures induce genome-wide mutations in the offspring that may be associated with adverse health outcomes.

Demodifying RNA for Transcriptomic Analyses of Archival Formalin-Fixed Paraffin-Embedded Samples.

Archival formalin-fixed paraffin-embedded (FFPE) tissue samples offer a vast but largely untapped resource for genomic research. The primary technical issues limiting use of FFPE samples are RNA yield and quality. In this study, we evaluated methods to demodify RNA highly fragmented and crosslinked by formalin fixation. Primary endpoints were RNA recovery, RNA-sequencing quality metrics, and transcriptional responses to a reference chemical (phenobarbital, PB). Frozen mouse liver samples from control and PB groups (n = 6/group) were divided and preserved for 3 months as follows: frozen (FR); 70% ethanol (OH); 10% buffered formalin for 18 h followed by ethanol (18F); or 10% buffered formalin (3F). Samples from OH, 18F, and 3F groups were processed to FFPE blocks and sectioned for RNA isolation. Additional sections from 3F received the following demodification protocols to mitigate RNA damage: short heated incubation with Tris-Acetate-EDTA buffer; overnight heated incubation with an organocatalyst using 2 different isolation kits; or overnight heated incubation without organocatalyst. Ribo-depleted, stranded, total RNA libraries were built and sequenced using the Illumina HiSeq 2500 platform. Overnight incubation (± organocatalyst) increased RNA yield >3-fold and RNA integrity numbers and fragment analysis values by > 1.5- and >3.0-fold, respectively, versus 3F. Postsequencing metrics also showed reduced bias in gene coverage and deletion rates for overnight incubation groups. All demodification groups had increased overlap for differentially expressed genes (77%-84%) and enriched pathways (91%-97%) with FR, with the highest overlap in the organocatalyst groups. These results demonstrate simple changes in RNA isolation methods that can enhance genomic analyses of FFPE samples.