Whole exome and transcript profiling of liver following aflatoxin B1 exposure in rats.

We recently developed a rat whole exome sequencing (WES) panel and used it to evaluate early somatic mutations in archival liver tissues from F344/N rats exposed to the hepatocarcinogen, Aflatoxin B1 (AFB1), a widely studied, potent mutagen and hepatocarcinogen associated with hepatocellular carcinoma (HCC). Rats were exposed to 1-ppm AFB1 in feed for 14, 90, and 90 days plus a recovery 60-day, non-exposure period (150-day) timepoint. Isolated liver DNA was exome sequenced. We identified 172 sequence variants across all timepoints, of which 101 were non-synonymous variants. Well-annotated genes carried a diverse set of 29 non-synonymous mutations at 14 days, increasing to 39 mutations at 90 days and then decreasing to 33 mutations following the 60-day recovery. Gene Set Enrichment Analysis conducted on previously reported, available RNA expression data of the same exome sequenced archival samples identified altered transcripts in pathways associated with malignant transformation. These included HALLMARK gene sets associated with cell proliferation (MYC Targets Version 1 and Version 2, E2F targets), cell cycle (G2M checkpoint, mitotic spindle), cell death (apoptosis), and DNA damage (DNA repair, UV response Up, Reactive oxygen species) pathways. DriverNet Impact analysis integrated exome-seq and expression data to reveal somatic mutations in Mcm8, Bdp1, and Cct6a that may drive cancer formation. Connectivity with transcript expression changes identified these genes as the top-ranked candidate driver genes associated with hepatocellular transformation. In conclusion, exome sequencing revealed early somatic mutations that may play a role in cancer cell transformation that are translatable to aflatoxin-induced HCC.

Insights into Repeated Renal Injury Using RNA-Seq with Two New RPTEC Cell Lines.

Renal proximal tubule epithelial cells (RPTECs) are a primary site for kidney injury. We created two RPTEC lines from CD-1 mice immortalized with hTERT (human telomerase reverse transcriptase) or SV40 LgT antigen (Simian Virus 40 Large T antigen). Our hypothesis was that low-level, repeated exposure to subcytotoxic levels of 0.25-2.5 µM cisplatin (CisPt) or 12.5-100 µM aflatoxin B1 (AFB1) would activate distinctive genes and pathways in these two differently immortalized cell lines. RNA-seq showed only LgT cells responded to AFB1 with 1139 differentially expressed genes (DEGs) at 72 h. The data suggested that AFB1 had direct nephrotoxic properties on the LgT cells. However, both the cell lines responded to 2.5 µM CisPt from 3 to 96 h expressing 2000-5000 total DEGs. For CisPt, the findings indicated a coordinated transcriptional program of injury signals and repair from the expression of immune receptors with cytokine and chemokine secretion for leukocyte recruitment; robust expression of synaptic and substrate adhesion molecules (SAMs) facilitating the expression of neural and hormonal receptors, ion channels/transporters, and trophic factors; and the expression of nephrogenesis transcription factors. Pathway analysis supported the concept of a renal repair transcriptome. In summary, these cell lines provide in vitro models for the improved understanding of repeated renal injury and repair mechanisms. High-throughput screening against toxicant libraries should provide a wider perspective of their capabilities in nephrotoxicity.

Single nucleotide polymorphism patterns associated with a cancer resistant phenotype.

BACKGROUND AND AIMS: In this study ten mouse strains representing ~90% of genetic diversity in laboratory mice (B6C3F1/J, C57BL/6J, C3H/HeJ, A/J, NOD.B1oSnH2/J, NZO/HILtJ, 129S1/SvImJ, WSB/EiJ, PWK/PhJ, CAST/EiJ) were examined to identify the mouse strain with the lowest incidence of cancer. The unique single polymorphisms (SNPs) associated with this low cancer incidence are reported. METHODS: Evaluations of cancer incidence in the 10 mouse strains were based on gross and microscopic diagnosis of tumors. Single nucleotide polymorphisms (SNPs) in the coding regions of the genome were derived from the respective mouse strains located in the Sanger mouse sequencing database and the B6C3F1/N genome from the National Toxicology Program (NTP). RESULTS: The WSB strain had an overall lower incidence of both benign and malignant tumors compared to the other mouse strains. At 2 years, the incidence of total malignant tumors (Poly-3 incidence rate) ranged from 2% (WSB) to 92% (C3H) in males, and 14% (WSB) to 93% (NZO) in females, and the total incidence of benign and malignant tumor incidence ranged from 13% (WSB) to 99% (C3H) in males and 25% (WSB) to 96% (NOD) in females. Single nucleotide polymorphism (SNP) patterns were examined in the following strains: B6C3F1/N, C57BL/6J, C3H/HeJ, 129S1/SvImJ, A/J, NZO/HILtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ. We identified 7519 SNPs (involving 5751 Ensembl transcripts of 3453 Ensembl Genes) that resulted in a unique amino acid change in the coding region of the WSB strain. CONCLUSIONS: The inherited genetic patterns in the WSB cancer-resistant mouse strain occurred in genes involved in multiple cell functions including mitochondria, metabolic, immune, and membrane-related cell functions. The unique SNP patterns in a cancer resistant mouse strain provides insights for understanding and developing strategies for cancer prevention.

Leucine-rich repeats and calponin homology containing 4 (Lrch4) regulates the innate immune response.

Toll-like receptors (TLRs) are pathogen-recognition receptors that trigger the innate immune response. Recent reports have identified accessory proteins that provide essential support to TLR function through ligand delivery and receptor trafficking. Herein, we introduce leucine-rich repeats (LRRs) and calponin homology containing 4 (Lrch4) as a novel TLR accessory protein. Lrch4 is a membrane protein with nine LRRs in its predicted ectodomain. It is widely expressed across murine tissues and has two expression variants that are both regulated by lipopolysaccharide (LPS). Predictive modeling indicates that Lrch4 LRRs conform to the horseshoe-shaped structure typical of LRRs in pathogen-recognition receptors and that the best structural match in the protein database is to the variable lymphocyte receptor of the jawless vertebrate hagfish. Silencing Lrch4 attenuates cytokine induction by LPS and multiple other TLR ligands and dampens the in vivo innate immune response. Lrch4 promotes proper docking of LPS in lipid raft membrane microdomains. We provide evidence that this is through regulation of lipid rafts as Lrch4 silencing reduces cell surface gangliosides, a metric of raft abundance, as well as expression and surface display of CD14, a raft-resident LPS co-receptor. Taken together, we identify Lrch4 as a broad-spanning regulator of the innate immune response and a potential molecular target in inflammatory disease.

BMDExpress 2: enhanced transcriptomic dose-response analysis workflow.

SUMMARY: A new version (version 2) of the genomic dose-response analysis software, BMDExpress, has been created. The software addresses the increasing use of transcriptomic dose-response data in toxicology, drug design, risk assessment and translational research. In this new version, we have implemented additional statistical filtering options (e.g. Williams' trend test), curve fitting models, Linux and Macintosh compatibility and support for additional transcriptomic platforms with up-to-date gene annotations. Furthermore, we have implemented extensive data visualizations, on-the-fly data filtering, and a batch-wise analysis workflow. We have also significantly re-engineered the code base to reflect contemporary software engineering practices and streamline future development. The first version of BMDExpress was developed in 2007 to meet an unmet demand for easy-to-use transcriptomic dose-response analysis software. Since its original release, however, transcriptomic platforms, technologies, pathway annotations and quantitative methods for data analysis have undergone a large change necessitating a significant re-development of BMDExpress. To that end, as of 2016, the National Toxicology Program assumed stewardship of BMDExpress. The result is a modernized and updated BMDExpress 2 that addresses the needs of the growing toxicogenomics user community. AVAILABILITY AND IMPLEMENTATION: BMDExpress 2 is available at https://github.com/auerbachs/BMDExpress-2/releases. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

KRAS-retroviral fusion transcripts and gene amplification in arsenic-transformed, human prostate CAsE-PE cancer cells.

CAsE-PE cells are an arsenic-transformed, human prostate epithelial line containing oncogenic mutations in KRAS compared to immortalized, normal KRAS parent cells, RWPE-1. We previously reported increased copy number of mutated KRAS in CAsE-PE cells, suggesting gene amplification. Here, KRAS flanking genomic and transcriptomic regions were sequenced in CAsE-PE cells for insight into KRAS amplification. Comparison of DNA-Seq and RNA-Seq showed increased reads from background aligning to all KRAS exons in CAsE-PE cells, while a uniform DNA-Seq read distribution occurred in RWPE-1 cells with normal transcript expression. We searched for KRAS fusions in DNA and RNA sequencing data finding a portion of reads aligning to KRAS and viral sequence. After generation of cDNA from total RNA, short and long KRAS probes were generated to hybridize cDNA and KRAS enriched fragments were PacBio sequenced. More KRAS reads were captured from CAsE-PE cDNA versus RWPE-1 by each probe set. Only CAsE-PE cDNA showed KRAS viral fusion transcripts, primarily mapping to LTR and endogenous retrovirus sequences on either 5'- or 3'-ends of KRAS. Most KRAS viral fusion transcripts contained 4 to 6 exons but some PacBio sequences were in unusual orientations, suggesting viral insertions within the gene body. Additionally, conditioned media was extracted for potential retroviral particles. RNA-Seq of culture media isolates identified KRAS retroviral fusion transcripts in CAsE-PE media only. Truncated KRAS transcripts suggested multiple retroviral integration sites occurred within the KRAS gene producing KRAS retroviral fusions of various lengths. Findings suggest activation of endogenous retroviruses in arsenic carcinogenesis should be explored.

Identifying Compounds with Genotoxicity Potential Using Tox21 High-Throughput Screening Assays.

Genotoxicity is a critical component of a comprehensive toxicological profile. The Tox21 Program used five quantitative high-throughput screening (qHTS) assays measuring some aspect of DNA damage/repair to provide information on the genotoxic potential of over 10 000 compounds. Included were assays detecting activation of p53, increases in the DNA repair protein ATAD5, phosphorylation of H2AX, and enhanced cytotoxicity in DT40 cells deficient in DNA-repair proteins REV3 or KU70/RAD54. Each assay measures a distinct component of the DNA damage response signaling network; >70% of active compounds were detected in only one of the five assays. When qHTS results were compared with results from three standard genotoxicity assays (bacterial mutation, in vitro chromosomal aberration, and in vivo micronucleus), a maximum of 40% of known, direct-acting genotoxicants were active in one or more of the qHTS genotoxicity assays, indicating low sensitivity. This suggests that these qHTS assays cannot in their current form be used to replace traditional genotoxicity assays. However, despite the low sensitivity, ranking chemicals by potency of response in the qHTS assays revealed an enrichment for genotoxicants up to 12-fold compared with random selection, when allowing a 1% false positive rate. This finding indicates these qHTS assays can be used to prioritize chemicals for further investigation, allowing resources to focus on compounds most likely to induce genotoxic effects. To refine this prioritization process, models for predicting the genotoxicity potential of chemicals that were active in Tox21 genotoxicity assays were constructed using all Tox21 assay data, yielding a prediction accuracy up to 0.83. Data from qHTS assays related to stress-response pathway signaling (including genotoxicity) were the most informative for model construction. By using the results from qHTS genotoxicity assays, predictions from models based on qHTS data, and predictions from commercial bacterial mutagenicity QSAR models, we prioritized Tox21 chemicals for genotoxicity characterization.

Identification of Compounds That Inhibit Estrogen-Related Receptor Alpha Signaling Using High-Throughput Screening Assays.

The nuclear receptor, estrogen-related receptor alpha (ERRα; NR3B1), plays a pivotal role in energy homeostasis. Its expression fluctuates with the demands of energy production in various tissues. When paired with the peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α), the PGC/ERR pathway regulates a host of genes that participate in metabolic signaling networks and in mitochondrial oxidative respiration. Unregulated overexpression of ERRα is found in many cancer cells, implicating a role in cancer progression and other metabolism-related diseases. Using high throughput screening assays, we screened the Tox21 10K compound library in stably transfected HEK293 cells containing either the ERRα-reporter or the reporter plus PGC-1α expression plasmid. We identified two groups of antagonists that were potent inhibitors of ERRα activity and/or the PGC/ERR pathway: nine antineoplastic agents and thirteen pesticides. Results were confirmed using gene expression studies. These findings suggest a novel mechanism of action on bioenergetics for five of the nine antineoplastic drugs. Nine of the thirteen pesticides, which have not been investigated previously for ERRα disrupting activity, were classified as such. In conclusion, we demonstrated that high-throughput screening assays can be used to reveal new biological properties of therapeutic and environmental chemicals, broadening our understanding of their modes of action.

Whole exome sequencing in the rat.

BACKGROUND: The rat genome was sequenced in 2004 with the aim to improve human health altered by disease and environmental influences through gene discovery and animal model validation. Here, we report development and testing of a probe set for whole exome sequencing (WES) to detect sequence variants in exons and UTRs of the rat genome. Using an in-silico approach, we designed probes targeting the rat exome and compared captured mutations in cancer-related genes from four chemically induced rat tumor cell lines (C6, FAT7, DSL-6A/C1, NBTII) to validated cancer genes in the human database, Catalogue of Somatic Mutations in Cancer (COSMIC) as well as normal rat DNA. Paired, fresh frozen (FF) and formalin-fixed, paraffin-embedded (FFPE) liver tissue from naive rats were sequenced to confirm known dbSNP variants and identify any additional variants. RESULTS: Informatics analysis of available gene annotation from rat RGSC6.0/rn6 RefSeq and Ensembl transcripts provided 223,636 unique exons representing a total of 26,365 unique genes and untranslated regions. Using this annotation and the Rn6 reference genome, an in-silico probe design generated 826,878 probe sequences of which 94.2% were uniquely aligned to the rat genome without mismatches. Further informatics analysis revealed 25,249 genes (95.8%) covered by at least one probe and 23,603 genes (93.5%) had every exon covered by one or more probes. We report high performance metrics from exome sequencing of our probe set and Sanger validation of annotated, highly relevant, cancer gene mutations as cataloged in the human COSMIC database, in addition to several exonic variants in cancer-related genes. CONCLUSIONS: An in-silico probe set was designed to enrich the rat exome from isolated DNA. The platform was tested on rat tumor cell lines and normal FF and FFPE liver tissue. The method effectively captured target exome regions in the test DNA samples with exceptional sensitivity and specificity to obtain reliable sequencing data representing variants that are likely chemically induced somatic mutations. Genomic discovery conducted by means of high throughput WES queries should benefit investigators in discovering rat genomic variants in disease etiology and in furthering human translational research.

Exome Sequencing of Fresh-frozen or Formalin-fixed Paraffin-embedded B6C3F1/N Mouse Hepatocellular Carcinomas Arising Either Spontaneously or due to Chronic Chemical Exposure.

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide; however, the mutational properties of HCC-associated carcinogens remain largely uncharacterized. We hypothesized that mechanisms underlying chemical-induced HCC can be characterized by evaluating the mutational spectra of these tumors. To test this hypothesis, we performed exome sequencing of B6C3F1/N HCCs that arose either spontaneously in vehicle controls ( n = 3) or due to chronic exposure to gingko biloba extract (GBE; n = 4) or methyleugenol (MEG; n = 3). Most archived tumor samples are available as formalin-fixed paraffin-embedded (FFPE) blocks, rather than fresh-frozen (FF) samples; hence, exome sequencing from paired FF and FFPE samples was compared. FF and FFPE samples showed 63% to 70% mutation concordance. Multiple known (e.g., Ctnnb1T41A, BrafV637E) and novel (e.g., Erbb4C559S, Card10A700V, and Klf11P358L) mutations in cancer-related genes were identified. The overall mutational burden was greater for MEG than for GBE or spontaneous HCC samples. To characterize the mutagenic mechanisms, we analyzed the mutational spectra in the HCCs according to their trinucleotide motifs. The MEG tumors clustered closest to Catalogue of Somatic Mutations in Cancer signatures 4 and 24, which are, respectively, associated with benzo(a)pyrene- and aflatoxin-induced HCCs in humans. These results establish a novel approach for classifying liver carcinogens and understanding the mechanisms of hepatocellular carcinogenesis.

Molecular Changes in the Nasal Cavity after N, N-dimethyl-p-toluidine Exposure.

N, N-dimethyl-p-toluidine (DMPT; Cas No. 99-97-8), an accelerant for methyl methacrylate monomers in medical devices, is a nasal cavity carcinogen according to a 2-yr cancer study of male and female F344/N rats, with the nasal tumors arising from the transitional cell epithelium. In this study, we exposed male F344/N rats for 5 days to DMPT (0, 1, 6, 20, 60, or 120 mg/kg [oral gavage]) to explore the early changes in the nasal cavity after short-term exposure. Lesions occurred in the nasal cavity including hyperplasia of transitional cell epithelium (60 and 120 mg/kg). Nasal tissue was rapidly removed and preserved for subsequent laser capture microdissection and isolation of the transitional cell epithelium (0 and 120 mg/kg) for transcriptomic studies. DMPT transitional cell epithelium gene transcript patterns were characteristic of an antioxidative damage response (e.g., Akr7a3, Maff, and Mgst3), cell proliferation, and decrease in signals for apoptosis. The transcripts of amino acid transporters were upregulated (e.g., Slc7a11). The DMPT nasal transcript expression pattern was similar to that found in the rat nasal cavity after formaldehyde exposure, with over 1,000 transcripts in common. Molecular changes in the nasal cavity after DMPT exposure suggest that oxidative damage is a mechanism of the DMPT toxic and/or carcinogenic effects.

Differential DNA methylation profile of key genes in malignant prostate epithelial cells transformed by inorganic arsenic or cadmium.

Previous work shows altered methylation patterns in inorganic arsenic (iAs)- or cadmium (Cd)-transformed epithelial cells. Here, the methylation status near the transcriptional start site was assessed in the normal human prostate epithelial cell line (RWPE-1) that was malignantly transformed by 10µM Cd for 11weeks (CTPE) or 5µM iAs for 29weeks (CAsE-PE), at which time cells showed multiple markers of acquired cancer phenotype. Next generation sequencing of the transcriptome of CAsE-PE cells identified multiple dysregulated genes. Of the most highly dysregulated genes, five genes that can be relevant to the carcinogenic process (S100P, HYAL1, NTM, NES, ALDH1A1) were chosen for an in-depth analysis of the DNA methylation profile. DNA was isolated, bisulfite converted, and combined bisulfite restriction analysis was used to identify differentially methylated CpG sites, which was confirmed with bisulfite sequencing. Four of the five genes showed differential methylation in transformants relative to control cells that was inversely related to altered gene expression. Increased expression of HYAL1 (>25-fold) and S100P (>40-fold) in transformants was correlated with hypomethylation near the transcriptional start site. Decreased expression of NES (>15-fold) and NTM (>1000-fold) in transformants was correlated with hypermethylation near the transcriptional start site. ALDH1A1 expression was differentially expressed in transformed cells but was not differentially methylated relative to control. In conclusion, altered gene expression observed in Cd and iAs transformed cells may result from altered DNA methylation status.

RNA-Seq-based toxicogenomic assessment of fresh frozen and formalin-fixed tissues yields similar mechanistic insights.

Formalin-fixed, paraffin-embedded (FFPE) pathology specimens represent a potentially vast resource for transcriptomic-based biomarker discovery. We present here a comparison of results from a whole transcriptome RNA-Seq analysis of RNA extracted from fresh frozen and FFPE livers. The samples were derived from rats exposed to aflatoxin B1 (AFB1 ) and a corresponding set of control animals. Principal components analysis indicated that samples were separated in the two groups representing presence or absence of chemical exposure, both in fresh frozen and FFPE sample types. Sixty-five percent of the differentially expressed transcripts (AFB1 vs. controls) in fresh frozen samples were also differentially expressed in FFPE samples (overlap significance: P < 0.0001). Genomic signature and gene set analysis of AFB1 differentially expressed transcript lists indicated highly similar results between fresh frozen and FFPE at the level of chemogenomic signatures (i.e., single chemical/dose/duration elicited transcriptomic signatures), mechanistic and pathology signatures, biological processes, canonical pathways and transcription factor networks. Overall, our results suggest that similar hypotheses about the biological mechanism of toxicity would be formulated from fresh frozen and FFPE samples. These results indicate that phenotypically anchored archival specimens represent a potentially informative resource for signature-based biomarker discovery and mechanistic characterization of toxicity.

Development of a stable cell line with an intact PGC-1α/ERRα axis for screening environmental chemicals.

The estrogen-related receptor α (ERRα) and the peroxisome proliferator-activated receptor γ (PPARγ) coactivator 1α (PGC-1α) play critical roles in the control of several physiological functions, including the regulation of genes involved in energy homeostasis. However, little is known about the ability of environmental chemicals to disrupt or modulate this important bioenergetics pathway in humans. The goal of this study was to develop a cell-based assay system with an intact PGC-1α/ERRα axis that could be used as a screening assay for detecting such chemicals. To this end, we successfully generated several stable cell lines expressing PGC-1α and showed that the reporter driven by the native ERRα hormone response unit (AAB-Luc) is active in these cell lines and that the activation is PGC-1α-dependent. Furthermore, we show that this activation can be blocked by the ERRα selective inverse agonist, XCT790. In addition, we find that genistein and bisphenol A further stimulate the reporter activity, while kaempferol has minimal effect. These cell lines will be useful for identifying environmental chemicals that modulate this important pathway.

Gene-environment interactions within a precision environmental health framework.

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

The use of faecal microbiota transplant as treatment for recurrent or refractory Clostridioides difficile infection and other potential indications: second edition of joint British Society of Gastroenterology (BSG) and Healthcare Infection Society (HIS) guidelines.

The first British Society of Gastroenterology (BSG) and Healthcare Infection Society (HIS)-endorsed faecal microbiota transplant (FMT) guidelines were published in 2018. Over the past 5 years, there has been considerable growth in the evidence base (including publication of outcomes from large national FMT registries), necessitating an updated critical review of the literature and a second edition of the BSG/HIS FMT guidelines. These have been produced in accordance with National Institute for Health and Care Excellence-accredited methodology, thus have particular relevance for UK-based clinicians, but are intended to be of pertinence internationally. This second edition of the guidelines have been divided into recommendations, good practice points and recommendations against certain practices. With respect to FMT for Clostridioides difficile infection (CDI), key focus areas centred around timing of administration, increasing clinical experience of encapsulated FMT preparations and optimising donor screening. The latter topic is of particular relevance given the COVID-19 pandemic, and cases of patient morbidity and mortality resulting from FMT-related pathogen transmission. The guidelines also considered emergent literature on the use of FMT in non-CDI settings (including both gastrointestinal and non-gastrointestinal indications), reviewing relevant randomised controlled trials. Recommendations are provided regarding special areas (including compassionate FMT use), and considerations regarding the evolving landscape of FMT and microbiome therapeutics.

Proteomic profiling of S-acylated macrophage proteins identifies a role for palmitoylation in mitochondrial targeting of phospholipid scramblase 3.

S-Palmitoylation, the reversible post-translational acylation of specific cysteine residues with the fatty acid palmitate, promotes the membrane tethering and subcellular localization of proteins in several biological pathways. Although inhibiting palmitoylation holds promise as a means for manipulating protein targeting, advances in the field have been hampered by limited understanding of palmitoylation enzymology and consensus motifs. In order to define the complement of S-acylated proteins in the macrophage, we treated RAW 264.7 macrophage membranes with hydroxylamine to cleave acyl thioesters, followed by biotinylation of newly exposed sulfhydryls and streptavidin-agarose affinity chromatography. Among proteins identified by LC-MS/MS, S-acylation status was established by spectral counting to assess enrichment under hydroxylamine versus mock treatment conditions. Of 1183 proteins identified in four independent experiments, 80 proteins were significant for S-acylation at false discovery rate = 0.05, and 101 significant at false discovery rate = 0.10. Candidate S-acylproteins were identified from several functional categories, including membrane trafficking, signaling, transporters, and receptors. Among these were 29 proteins previously biochemically confirmed as palmitoylated, 45 previously reported as putative S-acylproteins in proteomic screens, 24 not previously associated with palmitoylation, and three presumed false-positives. Nearly half of the candidates were previously identified by us in macrophage detergent-resistant membranes, suggesting that palmitoylation promotes lipid raft-localization of proteins in the macrophage. Among the candidate novel S-acylproteins was phospholipid scramblase 3 (Plscr3), a protein that regulates apoptosis through remodeling the mitochondrial membrane. Palmitoylation of Plscr3 was confirmed through (3)H-palmitate labeling. Moreover, site-directed mutagenesis of a cluster of five cysteines (Cys159-161-163-164-166) abolished palmitoylation, caused Plscr3 mislocalization from mitochondrion to nucleus, and reduced macrophage apoptosis in response to etoposide, together suggesting a role for palmitoylation at this site for mitochondrial targeting and pro-apoptotic function of Plscr3. Taken together, we propose that manipulation of protein palmitoylation carries great potential for intervention in macrophage biology via reprogramming of protein localization.

Chronic PFOA exposure in vitro causes acquisition of multiple tumor cell characteristics in rat liver cells.

Perfluorooctanoic acid (PFOA) is tumorigenic in rats and mice and potentially tumorigenic in humans. Here, we studied long-term PFOA exposure with an in vitro transformation model using the rat liver epithelial cell, TRL 1215. Cells were cultured in 10 µM (T10), 50 µM (T50) and 100 µM (T100) PFOA for 38 weeks and compared to passage-matched control cells. T100 cells showed morphological changes, loss of cell contact inhibition, formation of multinucleated giant and spindle-shaped cells. T10, T50, and T100 cells showed increased LC50 values 20%, 29% to 35% above control with acute PFOA treatment, indicating a resistance to PFOA toxicity. PFOA-treated cells showed increases in Matrix metalloproteinase-9 secretion, cell migration, and developed more and larger colonies in soft agar. Microarray data showed Myc pathway activation at T50 and T100, associating Myc upregulation with PFOA-induced morphological transformation. Western blot confirmed that PFOA produced significant increases in c-MYC protein expression in a time- and concentration-related manner. Tumor invasion indicators MMP-2 and MMP-9, cell cycle regulator cyclin D1, and oxidative stress protein GST were all significantly overexpressed in T100 cells. Taken together, chronic in vitro PFOA exposure produced multiple cell characteristics of malignant progression and differential gene expression changes suggestive of rat liver cell transformation.

Testing an aflatoxin B1 gene signature in rat archival tissues.

Archival tissues from laboratory studies represent a unique opportunity to explore the relationship between genomic changes and agent-induced disease. In this study, we evaluated the applicability of qPCR for detecting genomic changes in formalin-fixed, paraffin-embedded (FFPE) tissues by determining if a subset of 14 genes from a 90-gene signature derived from microarray data and associated with eventual tumor development could be detected in archival liver, kidney, and lung of rats exposed to aflatoxin B1 (AFB1) for 90 days in feed at 1 ppm. These tissues originated from the same rats used in the microarray study. The 14 genes evaluated were Adam8, Cdh13, Ddit4l, Mybl2, Akr7a3, Akr7a2, Fhit, Wwox, Abcb1b, Abcc3, Cxcl1, Gsta5, Grin2c, and the C8orf46 homologue. The qPCR FFPE liver results were compared to the original liver microarray data and to qPCR results using RNA from fresh frozen liver. Archival liver paraffin blocks yielded 30 to 50 µg of degraded RNA that ranged in size from 0.1 to 4 kB. qPCR results from FFPE and fresh frozen liver samples were positively correlated (p ≤ 0.05) by regression analysis and showed good agreement in direction and proportion of change with microarray data for 11 of 14 genes. All 14 transcripts could be amplified from FFPE kidney RNA except the glutamate receptor gene Grin2c; however, only Abcb1b was significantly upregulated from control. Abundant constitutive transcripts, S18 and ß-actin, could be amplified from lung FFPE samples, but the narrow RNA size range (25-500 bp length) prevented consistent detection of target transcripts. Overall, a discrete gene signature derived from prior transcript profiling and representing cell cycle progression, DNA damage response, and xenosensor and detoxication pathways was successfully applied to archival liver and kidney by qPCR and indicated that gene expression changes in response to subchronic AFB1 exposure occurred predominantly in the liver, the primary target for AFB1-induced tumors. We conclude that an evaluation of gene signatures in archival tissues can be an important toxicological tool for evaluating critical molecular events associated with chemical exposures.

Deep Learning Image Analysis of High-Throughput Toxicology Assay Images.

High-throughput chemical screening approaches often employ microscopy to capture photomicrographs from multi-well cell culture plates, generating thousands of images that require time-consuming human analysis. To automate this subjective and time-consuming manual process, we have developed a method that uses deep learning to automatically classify digital assay images. We have trained a convolutional neural network (CNN) to perform binary and multi-class classification. The binary classifier binned assay images into healthy (comparable to untreated controls) and altered (not comparable to untreated-control) classes with >98% accuracy; the multi-class classifier assigned "Healthy," "Intermediate" and "Altered" labels to assay images with >95% accuracy. Our dataset comprised high-resolution assay images from primary human hepatocytes and undifferentiated (proliferating) and differentiated 2D cultures of HepaRG cells. In this study we have focused on testing and fine-tuning various CNN architectures, including ResNet 34, 50 and 101. To visualize regions in the images that the CNN model used for classification, we employed Class Activation Maps (CAM). This allowed us to better understand the inner workings of the neural network and led to additional optimizations of the algorithm. The results indicate a strong correspondence between dosage and classifier-predicted scores, suggesting that these scores might be useful in further characterizing benchmark dose. Together, these results clearly demonstrate that deep-learning based automated image classification of cell morphology changes upon chemical-induced stress can yield highly accurate and reproducible assessments of cytotoxicity across a variety of cell types.