Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes.

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

The power of TOPMed imputation for the discovery of Latino-enriched rare variants associated with type 2 diabetes.

AIMS/HYPOTHESIS: The Latino population has been systematically underrepresented in large-scale genetic analyses, and previous studies have relied on the imputation of ungenotyped variants based on the 1000 Genomes (1000G) imputation panel, which results in suboptimal capture of low-frequency or Latino-enriched variants. The National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) released the largest multi-ancestry genotype reference panel representing a unique opportunity to analyse rare genetic variations in the Latino population. We hypothesise that a more comprehensive analysis of low/rare variation using the TOPMed panel would improve our knowledge of the genetics of type 2 diabetes in the Latino population. METHODS: We evaluated the TOPMed imputation performance using genotyping array and whole-exome sequence data in six Latino cohorts. To evaluate the ability of TOPMed imputation to increase the number of identified loci, we performed a Latino type 2 diabetes genome-wide association study (GWAS) meta-analysis in 8150 individuals with type 2 diabetes and 10,735 control individuals and replicated the results in six additional cohorts including whole-genome sequence data from the All of Us cohort. RESULTS: Compared with imputation with 1000G, the TOPMed panel improved the identification of rare and low-frequency variants. We identified 26 genome-wide significant signals including a novel variant (minor allele frequency 1.7%; OR 1.37, p=3.4 × 10-9). A Latino-tailored polygenic score constructed from our data and GWAS data from East Asian and European populations improved the prediction accuracy in a Latino target dataset, explaining up to 7.6% of the type 2 diabetes risk variance. CONCLUSIONS/INTERPRETATION: Our results demonstrate the utility of TOPMed imputation for identifying low-frequency variants in understudied populations, leading to the discovery of novel disease associations and the improvement of polygenic scores. DATA AVAILABILITY: Full summary statistics are available through the Common Metabolic Diseases Knowledge Portal ( https://t2d.hugeamp.org/downloads.html ) and through the GWAS catalog ( https://www.ebi.ac.uk/gwas/ , accession ID: GCST90255648). Polygenic score (PS) weights for each ancestry are available via the PGS catalog ( https://www.pgscatalog.org , publication ID: PGP000445, scores IDs: PGS003443, PGS003444 and PGS003445).

Characterizing Serpinb2 as a Modulator of TCDD-Induced Suppression of the B Cell.

2,3,7,8-Tetrachlordibenzo- p-dioxin (TCDD) is an environmental pollutant that can cause various toxic effects, including chloracne, metabolic syndrome, and immune suppression. Most of the toxicity associated with TCDD is mediated through activation of the aryl hydrocarbon receptor (AHR). Recent research has suggested the presence of a wide-range of interindividual variability in TCDD-mediated suppression of the Immunoglobulin-M (IgM) response across the human population. In an attempt to identify putative modifiers of AHR-mediated immunosuppression beyond the AHR, B cells were isolated from a panel of genetically diverse mouse strain to scan for modulators that drive interstrain differences in TCDD-mediated suppression of the IgM response. Results implicated a region of mouse Chromosome 1 near a gene encoding serine peptidase inhibitor, clade B, member 2 ( Serpinb2) whose human ortholog is plasminogen activator inhibitor 2 (PAI2). Further downstream analyses indicated that Serpinb2 is dysregulated by TCDD and, furthermore, that B cells from Serpinb2 -/- mice are significantly more sensitive to TCDD-mediated suppression as compared to littermate controls. This study suggests a protective role of Serpinb2 within TCDD-mediated immunosuppression and, furthermore, a novel function of Serpinb2-related activity in the IgM response.

Mitochondrial-targeted aryl hydrocarbon receptor and the impact of 2,3,7,8-tetrachlorodibenzo-p-dioxin on cellular respiration and the mitochondrial proteome.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor within the Per-Arnt-Sim (PAS) domain superfamily. Exposure to the most potent AHR ligand, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), is associated with various pathological effects including metabolic syndrome. While research over the last several years has demonstrated a role for oxidative stress and metabolic dysfunction in AHR-dependent TCDD-induced toxicity, the role of the mitochondria in this process has not been fully explored. Our previous research suggested that a portion of the cellular pool of AHR could be found in the mitochondria (mitoAHR). Using a protease protection assay with digitonin extraction, we have now shown that this mitoAHR is localized to the inter-membrane space (IMS) of the organelle. TCDD exposure induced a degradation of mitoAHR similar to that of cytosolic AHR. Furthermore, siRNA-mediated knockdown revealed that translocase of outer-mitochondrial membrane 20 (TOMM20) was involved in the import of AHR into the mitochondria. In addition, TCDD altered cellular respiration in an AHR-dependent manner to maintain respiratory efficiency as measured by oxygen consumption rate (OCR). Stable isotope labeling by amino acids in cell culture (SILAC) identified a battery of proteins within the mitochondrial proteome influenced by TCDD in an AHR-dependent manner. Among these, 17 proteins with fold changes≥2 are associated with various metabolic pathways, suggesting a role of mitochondrial retrograde signaling in TCDD-mediated pathologies. Collectively, these studies suggest that mitoAHR is localized to the IMS and AHR-dependent TCDD-induced toxicity, including metabolic dysfunction, wasting syndrome, and hepatic steatosis, involves mitochondrial dysfunction.

Hepatic polybrominated diphenyl ether (PBDE) levels in Wisconsin river otters (Lontra canadensis) and Michigan bald eagles (Haliaeetus leucocephalus).

Polybrominated diphenyl ethers (PBDEs) are persistent and toxic flame-retardant chemicals widespread in the Great Lakes ecosystem. These chemicals are now being regulated and phased-out of the region; therefore it remains important to understand the extent of contamination in order to track the efficacy of recent actions. Here, Σ4PBDE congeners (PBDE-47, 99, 100, 153;wetweight basis unless indicated)were determined in liver tissues from Wisconsin river otters (Lontra canadensis; n = 35; 2009-2010) and Michigan bald eagles (Haliaeetus leucocephalus; n = 33; 2009-2011). In otters, Σ4PBDE ranged from0.5 to 72.9 ng/g, with a mean (±SD) and median (25th-75th percentile inter-quartile range) of 16.3 ± 16.4 ng/g and 11.3 (5.6-18.9) ng/g, respectively. The mean lipid-adjusted Σ4PBDE was 1377 ± 1485 ng/g. In eagles, Σ4PBDE ranged from 0 to 1,538.8 ng/g, with a mean and median of 74.3 ± 266.7 ng/g and 21.2 (5.7-28.9) ng/g, respectively. The mean lipid-adjusted Σ4PBDE was 5274.5 ± 19,896.1 ng/g. In both species, PBDE-47 accounted for >50% of the Σ4PBDE, followed by PBDE-99 and PBDE-100 (each ~17-19% of the total). The PBDE levels reported here in otters are similar to mammalian wildlife elsewhere, though the levels in eagles are among the highest worldwide across studied birds. The findings indicate that apex Great Lakes wildlife remain exposed to appreciable levels of PBDEs and more work is needed to understand whether such exposures are associated with adverse health outcomes.

Mercury exposure and neurochemical biomarkers in multiple brain regions of Wisconsin river otters (Lontra canadensis).

River otters are fish-eating wildlife that bioaccumulate high levels of mercury (Hg). Mercury is a proven neurotoxicant to mammalian wildlife, but little is known about the underlying, sub-clinical effects. Here, the overall goal was to increase understanding of Hg's neurological risk to otters. First, Hg values across several brain regions and tissues were characterized. Second, in three brain regions with known sensitivity to Hg (brainstem, cerebellum, and occipital cortex), potential associations among Hg levels and neurochemical biomarkers [N-methyl-D-aspartic acid (NMDA) and gamma-aminobutyric acid (GABAA) receptor] were explored. There were no significant differences in Hg levels across eight brain regions (rank order, highest to lowest: frontal cortex, cerebellum, temporal cortex, occipital cortex, parietal cortex, basal ganglia, brainstem, and thalamus), with mean values ranging from 0.7 to 1.3 ug/g dry weight. These brain levels were significantly lower than mean values in the muscle (2.1 ± 1.4 ug/g), liver (4.7 ± 4.3 ug/g), and fur (8.8 ± 4.8 ug/g). While a significant association was found between Hg and NMDA receptor levels in the brain stem (P = 0.028, rp = -0.293), no relationships were found in the cerebellum and occipital cortex. For the GABA receptor, no relationships were found. The lack of consistent Hg-associated neurochemical changes is likely due to low brain Hg levels in these river otters, which are amongst the lowest reported.

The Type 2 Diabetes Knowledge Portal: An open access genetic resource dedicated to type 2 diabetes and related traits.

Associations between human genetic variation and clinical phenotypes have become a foundation of biomedical research. Most repositories of these data seek to be disease-agnostic and therefore lack disease-focused views. The Type 2 Diabetes Knowledge Portal (T2DKP) is a public resource of genetic datasets and genomic annotations dedicated to type 2 diabetes (T2D) and related traits. Here, we seek to make the T2DKP more accessible to prospective users and more useful to existing users. First, we evaluate the T2DKP's comprehensiveness by comparing its datasets with those of other repositories. Second, we describe how researchers unfamiliar with human genetic data can begin using and correctly interpreting them via the T2DKP. Third, we describe how existing users can extend their current workflows to use the full suite of tools offered by the T2DKP. We finally discuss the lessons offered by the T2DKP toward the goal of democratizing access to complex disease genetic results.

Evaluating human genetic support for hypothesized metabolic disease genes.

We investigate the extent to which human genetic data are incorporated into studies that hypothesize novel links between genes and metabolic disease. To lower the barriers to using genetic data, we present an approach to enable researchers to evaluate human genetic support for experimentally determined hypotheses.

A combined polygenic score of 21,293 rare and 22 common variants improves diabetes diagnosis based on hemoglobin A1C levels.

Polygenic scores (PGSs) combine the effects of common genetic variants1,2 to predict risk or treatment strategies for complex diseases3-7. Adding rare variation to PGSs has largely unknown benefits and is methodically challenging. Here, we developed a method for constructing rare variant PGSs and applied it to calculate genetically modified hemoglobin A1C thresholds for type 2 diabetes (T2D) diagnosis7-10. The resultant rare variant PGS is highly polygenic (21,293 variants across 154 genes), depends on ultra-rare variants (72.7% observed in fewer than three people) and identifies significantly more undiagnosed T2D cases than expected by chance (odds ratio = 2.71; P = 1.51 × 10-6). A PGS combining common and rare variants is expected to identify 4.9 million misdiagnosed T2D cases in the United States-nearly 1.5-fold more than the common variant PGS alone. These results provide a method for constructing complex trait PGSs from rare variants and suggest that rare variants will augment common variants in precision medicine approaches for common disease.

Rare loss of function variants in the hepatokine gene INHBE protect from abdominal obesity.

Identifying genetic variants associated with lower waist-to-hip ratio can reveal new therapeutic targets for abdominal obesity. We use exome sequences from 362,679 individuals to identify genes associated with waist-to-hip ratio adjusted for BMI (WHRadjBMI), a surrogate for abdominal fat that is causally linked to type 2 diabetes and coronary heart disease. Predicted loss of function (pLOF) variants in INHBE associate with lower WHRadjBMI and this association replicates in data from AMP-T2D-GENES. INHBE encodes a secreted protein, the hepatokine activin E. In vitro characterization of the most common INHBE pLOF variant in our study, indicates an in-frame deletion resulting in a 90% reduction in secreted protein levels. We detect associations with lower WHRadjBMI for variants in ACVR1C, encoding an activin receptor, further highlighting the involvement of activins in regulating fat distribution. These findings highlight activin E as a potential therapeutic target for abdominal obesity, a phenotype linked to cardiometabolic disease.

Peanut butter as an alternative dose delivery method to prevent strain-dependent orogastric gavage-induced stress in mouse teratogenicity studies.

INTRODUCTION: Animal-based studies are essential for assessing toxicity to environmental pollutants, especially when the potential targets are specific developmental time points, teratogenic, or multi-organ systems that cannot be modeled in vitro. Orogastric gavage is a widely used technique for exposure because of its increased accuracy of dose administration over free feeding. However, repeated use of this method has been reported to cause physiological stress on the exposed animals that could interfere with interpretation of results. Previous studies have shown that genetic background also contributes to the level of stress and can affect individual response. METHODS: To evaluate the impact of stress on repeated orogastric gavage, we exposed C67BL/6J and 129S1/SvImJ inbred mouse strains to 2,3,7,8 tetrachlorodibenzo-p-dioxin (TCDD), a potent xenobiotic that has been extensively studied in vivo. Pregnant females were dosed for ten days after mating using orogastric gavage with olive oil as vehicle or through diet using peanut butter as vehicle. Serum corticosterone levels, body weight, and reproduction endpoints were measured to evaluate levels of stress induced by the dosing technique. RESULTS: The levels of stress caused by orogastric gavage was strongly dependent on strain background and on the phenotypic endpoint. Orogastric gavage-induced stress was more detrimental in 129S1/SvlmJ pregnant female mice than in C57BL/6J. CONCLUSION: These results show that administration of xenobiotics via controlled diet can improve the reproducibility and rigor of exposure studies requiring orogastric delivery.

Genetics-Based Approach to Identify Novel Genes Regulated by the Aryl Hydrocarbon Receptor in Mouse Liver.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor in the Per-Arnt-Sim superfamily of environmental sensors that is linked to several metabolic diseases, including nonalcoholic fatty liver disease. Much remains unknown regarding the impact of genetic variation in AHR-driven disease, as past studies have focused on a small number of inbred strains. Recently, the presence of a wide range of interindividual variability amongst humans was reported in response to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), the prototypical ligand of the AHR. In this study, a panel of 14 diverse mouse strains was exposed to TCDD for 10 days to characterize the AHR-mediated response across genetic backgrounds. Responses to TCDD are heavily dependent on genetic background. Although mice carry 1 of 4 Ahr alleles known to impact the affinity to AHR-ligands, we observed significant intra-allelic variability suggesting the presence of novel genetic modifiers of AHR signaling. A regression-based approach was used to scan for genes regulated by the AHR and/or associated with TCDD-induced phenotypes. The approach identified 7 genes, 2 of which are novel, that are likely regulated by the AHR based on association with hepatic TCDD burden (p ≤ .05). Finally, we identified 1 gene, Dio1, which was associated with change in percent body fat across the diverse set of strains (p ≤ .05). Overall, the results in this study exemplify the power of genetics-based approaches in identifying novel genes that are putatively regulated by the AHR.

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes.

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

An Automated Method To Predict Mouse Gene and Protein Sequences Using Variant Data.

With recent advances in sequencing technologies, the scientific community has begun to probe the potential genetic bases behind complex phenotypes in humans and model organisms. In many cases, the genomes of genetically distinct strains of model organisms, such as the mouse (Mus musculus), have not been fully sequenced. Here, we report on a tool designed to use single-nucleotide polymorphism (SNP) and insertion-deletion (indel) data to predict gene, mRNA, and protein sequences for up to 36 genetically distinct mouse strains. By automated querying of freely accessible databases through a graphical interface, the software requires no data and little computational experience. As a proof of concept, we predicted the gene and amino acid sequence of the aryl hydrocarbon receptor (Ahr) for all inbred mouse strains of which variant data were currently available through Mouse Genome Project. Predicted sequences were compared with fully sequenced genomes to show that the tool is effective in predicting gene and protein sequences.

Characterizing the Role of HMG-CoA Reductase in Aryl Hydrocarbon Receptor-Mediated Liver Injury in C57BL/6 Mice.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor. The prototypical ligand of the AHR is an environmental contaminant called 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). TCDD exposure is associated with many adverse health outcomes in humans including non-alcoholic fatty liver disease (NAFLD). Previous studies suggest that AHR ligands alter cholesterol homeostasis in mice through repression of genes involved in cholesterol biosynthesis, such as Hmgcr, which encodes the rate-limiting enzyme of cholesterol biosynthesis called 3-hydroxy-3-methyl-glutaryl coenzyme A reductase (HMGCR). In this study, we sought to characterize the impact of HMGCR repression in TCDD-induced liver injury. C57BL/6 mice were exposed to TCDD in the presence or absence of simvastatin, a competitive inhibitor of HMGCR. Simvastatin exposure decreased TCDD-induced hepatic lipid accumulation in both sexes, but was most prominent in females. Simvastatin and TCDD (S + T) co-treatment increased hepatic AHR-battery gene expression and liver injury in male, but not female, mice. In addition, the S + T co-treatment led to an increase in hepatic glycogen content that coincides with heavier liver in female mice. Results from this study suggest that statins, which are amongst the most prescribed pharmaceuticals, may protect from AHR-mediated steatosis, but alter glycogen metabolism and increase the risk of TCDD-elicited liver damage in a sex-specific manner.

Incorporating population-level genetic variability within laboratory models in toxicology: From the individual to the population.

Humans respond to chemical exposures differently due to many factors, such as previous and concurrent stressors, age, sex, and genetic background. The vast majority of laboratory-based toxicology studies, however, have not considered the impact of population-level variability within dose-response relationships. The lack of data dealing with the influence of genetic diversity on the response to chemical exposure provides a difficult challenge for risk assessment as individuals within the population will display a wide-range of responses following toxicant challenge. Notably, the genetic background of individuals plays a major role in the variability seen in a population-level response to a drug or chemical and, thus, there is growing interest in including genetic diversity into laboratory-models. Here we outline several laboratory-based models that can be used to assay the influence of genetic variability on an individual's response to chemicals: 1) genetically-diverse cell lines, 2) human primary cells, 3) and genetically-diverse mouse panels. We also provide a succinct review for several seminal studies to highlight the capability, feasibility, and power of each of these models. This article is intended to highlight the need to include population-level genetic diversity into toxicological study designs via laboratory-based models with the goal to provide and supplement evidence in assessing the risk posed by chemicals to the human population. As such, incorporation of genetic variability will positively impact human-based risk assessment and provide empirical data to aid and influence decision-making processes in relation to chemical exposures.

A cell-free testing platform to screen chemicals of potential neurotoxic concern across twenty vertebrate species.

There is global demand for new in vitro testing tools for ecological risk assessment. The objective of the present study was to apply a set of cell-free neurochemical assays to screen many chemicals across many species in a relatively high-throughput manner. The platform assessed 7 receptors and enzymes that mediate neurotransmission of γ-aminobutyric acid, dopamine, glutamate, and acetylcholine. Each assay was optimized to work across 20 vertebrate species (5 fish, 5 birds, 7 mammalian wildlife, 3 biomedical species including humans). We tested the screening assay platform against 80 chemicals (23 pharmaceuticals and personal care products, 20 metal[loid]s, 22 polycyclic aromatic hydrocarbons and halogenated organic compounds, 15 pesticides). In total, 10 800 species-chemical-assay combinations were tested, and significant differences were found in 4041 cases. All 7 assays were significantly affected by at least one chemical in each species tested. Among the 80 chemicals tested, nearly all resulted in a significant impact on at least one species and one assay. The 5 most active chemicals were prochloraz, HgCl2 , Sn, benzo[a]pyrene, and vinclozolin. Clustering analyses revealed groupings according to chemicals, species, and chemical-assay combinations. The results show that cell-free assays can screen a large number of samples in a short period of time in a cost-effective manner in a range of animals not easily studied using traditional approaches. Strengths and limitations of this approach are discussed, as well as next steps. Environ Toxicol Chem 2017;36:3081-3090. © 2017 SETAC.

Data on AHR-dependent changes in the mitochondrial proteome in response to ,3,7,8-tetrachlorodibenzo-p-dioxin.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that is the principal regulator of a cell׳s response to many polyaromatic hydrocarbons, such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). To gain a better understanding of the impact of TCDD on the mitochondrial proteome, a stable isotope labeling by amino acids in cell culture (SILAC)-based proteomic analysis was performed. We used two mouse hepatoma cell lines that differ in AHR expression levels, hepa1c1c7 (AHR-expressing) and hepac12 (AHR-deficient). The cell lines were exposed to TCDD (10 nM) for 72 h; each treatment was assayed in triplicate and were analyzed as separate runs on the mass-spectrometer. Mitochondria were then isolated and mitochondrial proteins were separated by SDS-PAGE and subject to mass spectrometry. The data presented were collected from four independent SILAC experiments. Within each experiment, three isotopes were employed to compare protein ratios via mass-spectrometry: (1) light l-arginine/l-lysine HCl (Arg0, Lys0), (2) medium (15)N4-l-arginin/(13)C6l-lysine HCl (Arg4, Lys6), and (3) heavy (13)C6 (15)N4l-arginine/(13)C6 (15)N2l-lysine HCl (Arg10, Lys8). The raw data includes approximately 2500 annotated proteins. The datasets provided by this study can be a reference to other toxicologists investigating TCDD-induced mitochondrial dysfunction. The data presented here are associated with the research article, "Mitochondrial-targeted Aryl Hydrocarbon Receptor and the Impact of 2,3,7,8-Tetrachlorodibenzo-p-Dioxin on Cellular Respiration and the Mitochondrial Proteome" (Hwang et al. (2016) [1]).

The Influence of Human Interindividual Variability on the Low-Dose Region of Dose-Response Curve Induced by 2,3,7,8-Tetrachlorodibenzo-p-Dioxin in Primary B Cells.

The influence of interindividual variability is not typically assessed in traditional toxicological studies. Given that chemical exposures occur in heterogeneous populations, this knowledge gap has the potential to cause undue harm within the realms of public health and industrial and municipal finances. A recent report from the National Research Council (NRC) suggests that when accounting for interindividual variation in responses, traditionally assumed nonlinear dose-response relationships (DRRs) for noncancer-causing endpoints would better be explained with a linear relationship within the low-dose region. To address this knowledge gap and directly test the NRC's assumption, this study focused on assessing the DRR between 2,3,7,8-tetracholorodibenzo-p-dioxin (TCDD) exposure and immune suppression in a cohort of unique human donors. Human B cells were isolated from 51 individual donors and treated with logarithmically increasing concentrations of TCDD (0-30 nM TCDD). Two endpoints sensitive to TCDD were assessed: (1) number of IgM-secreting B cells and (2) quantity of IgM secreted. The results show that TCDD significantly suppressed both the number of IgM-secreting B cells and the quantity of IgM secreted (P < .05). Statistical model comparisons indicate that the low-dose region of the two DRRs is best explained with a nonlinear relationship. Rather than assuming low-dose linearity for all noncancer-causing DRRs, our study suggests the need to consider the specific mode of action of toxicants and pharmaceuticals during risk-management decision making.

Dose-Dependent Metabolic Reprogramming and Differential Gene Expression in TCDD-Elicited Hepatic Fibrosis.

We have previously shown that in response to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-elicited NAFLD progression, central carbon, glutaminolysis, and serine/folate metabolism are reprogrammed to support NADPH production and ROS defenses. To further investigate underlying dose-dependent responses associated with TCDD-induced fibrosis, female C57BL/6 mice were gavaged with TCDD every 4 days (d) for 28 d or 92 d. RNA-Seq, ChIP-Seq (2 h), and 28 d metabolomic (urine, serum, and hepatic extract) analyses were conducted with complementary serum marker assessments at 92 d. Additional vehicle and 30 µg/kg treatment groups were allowed to recover for 36 d following the 92-d treatment regimen to examine recovery from TCDD-elicited fibrosis. Histopathology revealed dose-dependent increases in hepatic fat accumulation, inflammation, and periportal collagen deposition at 92 days, with increased fibrotic severity in the recovery group. Serum proinflammatory and profibrotic interleukins-1ß, -2, -4, -6, and -10, as well as TNF-α and IFN-γ, exhibited dose-dependent induction. An increase in glucose tolerance was observed with a concomitant 3.0-fold decrease in hepatic glycogen linked to increased ascorbic acid biosynthesis and proline metabolism, consistent with increased fibrosis. RNA-Seq identified differential expression of numerous matrisome genes including an 8.8-fold increase in Tgfb2 indicating myofibroblast activation. Further analysis suggests reprogramming of glycogen, ascorbic acid, and amino acid metabolism in support of collagen deposition and the use of proline as a substrate for ATP production via the proline cycle. In summary, we demonstrate that glycogen, ascorbic acid, and amino acid metabolism are also reorganized to support remodeling of the extracellular matrix, progressing to hepatic fibrosis in response to chronic injury from TCDD.