Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

Ten years of using key characteristics of human carcinogens to organize and evaluate mechanistic evidence in IARC Monographs on the identification of carcinogenic hazards to humans: Patterns and associations.

Systematic review and evaluation of mechanistic evidence using the Key Characteristics approach was proposed by the International Agency for Research on Cancer (IARC) in 2012 and used by the IARC Monographs Working Groups since 2015. Key Characteristics are 10 features of agents known to cause cancer in humans. From 2015 to 2022, a total of 19 Monographs (73 agents combined) used Key Characteristics for cancer hazard classification. We hypothesized that a retrospective analysis of applications of the Key Characteristics approach to cancer hazard classification using heterogenous mechanistic data on diverse agents would be informative for systematic reviews in decision-making. We extracted information on the conclusions, data types, and the role mechanistic data played in the cancer hazard classification from each Monograph. Statistical analyses identified patterns in the use of Key Characteristics, as well as trends and correlations among Key Characteristics, data types, and ultimate decisions. Despite gaps in data for many agents and Key Characteristics, several significant results emerged. Mechanistic data from in vivo animal, in vitro animal, and in vitro human studies were most impactful in concluding that an agent could cause cancer via a Key Characteristic. To exclude the involvement of a Key Characteristic, data from large-scale systematic in vitro testing programs such as ToxCast, were most informative. Overall, increased availability of systemized data streams, such as human in vitro data, would provide the basis for more confident and informed conclusions about both positive and negative associations and inform expert judgments on cancer hazard.

Application of Ion Mobility Spectrometry-Mass Spectrometry for Compositional Characterization and Fingerprinting of a Library of Diverse Crude Oil Samples.

Exposure characterization of crude oils, especially in time-sensitive circumstances such as spills and disasters, is a well-known analytical chemistry challenge. Gas chromatography-mass spectrometry is commonly used for "fingerprinting" and origin tracing in oil spills; however, this method is both time-consuming and lacks the resolving power to separate co-eluting compounds. Recent advances in methodologies to analyze petroleum substances using high-resolution analytical techniques have demonstrated both improved resolving power and higher throughput. One such method, ion mobility spectrometry-mass spectrometry (IMS-MS), is especially promising because it is both rapid and high-throughput, with the ability to discern among highly homologous hydrocarbon molecules. Previous applications of IMS-MS to crude oil analyses included a limited number of samples and did not provide detailed characterization of chemical constituents. We analyzed a diverse library of 195 crude oil samples using IMS-MS and applied a computational workflow to assign molecular formulas to individual features. The oils were from 12 groups based on geographical and geological origins: non-US (1 group), US onshore (3), and US Gulf of Mexico offshore (8). We hypothesized that information acquired through IMS-MS data would provide a more confident grouping and yield additional fingerprint information. Chemical composition data from IMS-MS was used for unsupervised hierarchical clustering, as well as machine learning-based supervised analysis to predict geographic and source rock categories for each sample; the latter also yielded several novel prospective biomarkers for fingerprinting of crude oils. We found that IMS-MS data have complementary advantages for fingerprinting and characterization of diverse crude oils and that proposed polycyclic aromatic hydrocarbon biomarkers can be used for rapid exposure characterization. Environ Toxicol Chem 2023;42:2336-2349. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Ten Years of Using Key Characteristics of Human Carcinogens to Organize and Evaluate Mechanistic Evidence in IARC Monographs on the Identification of Carcinogenic Hazards to Humans: Patterns and Associations.

Systematic review and evaluation of the mechanistic evidence only recently been instituted in cancer hazard identification step of decision-making. One example of organizing and evaluating mechanistic evidence is the Key Characteristics approach of the International Agency for Research on Cancer (IARC) Monographs on the Identification of Carcinogenic Hazards to Humans. The Key Characteristics of Human Carcinogens were proposed almost 10 years ago and have been used in every IARC Monograph since 2015. We investigated the patterns and associations in the use of Key Characteristics by the independent expert Working Groups. We examined 19 Monographs (2015-2022) that evaluated 73 agents. We extracted information on the conclusions by each Working Group on the strength of evidence for agent-Key Characteristic combinations, data types that were available for decisions, and the role mechanistic data played in the final cancer hazard classification. We conducted both descriptive and association analyses within and across data types. We found that IARC Working Groups were cautious when evaluating mechanistic evidence: for only ∼13% of the agents was strong evidence assigned for any Key Characteristic. Genotoxicity and cell proliferation were most data-rich, while little evidence was available for DNA repair and immortalization Key Characteristics. Analysis of the associations among Key Characteristics revealed that only chemical's metabolic activation was significantly co-occurring with genotoxicity and cell proliferation/death. Evidence from exposed humans was limited, while mechanistic evidence from rodent studies in vivo was often available. Only genotoxicity and cell proliferation/death were strongly associated with decisions on whether mechanistic data was impactful on the final cancer hazard classification. The practice of using the Key Characteristics approach is now well-established at IARC Monographs and other government agencies and the analyses presented herein will inform the future use of mechanistic evidence in regulatory decision-making.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients.

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

Pleiotropic modifiers of age-related diabetes and neonatal intestinal obstruction in cystic fibrosis.

Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting ∼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF.

MagicalRsq: Machine-learning-based genotype imputation quality calibration.

Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub.

A resource for integrated genomic analysis of the human liver.

In this study, we generated whole-transcriptome RNA-Seq from n = 192 genotyped liver samples and used these data with existing data from the GTEx Project (RNA-Seq) and previous liver eQTL (microarray) studies to create an enhanced transcriptomic sequence resource in the human liver. Analyses of genotype-expression associations show pronounced enrichment of associations with genes of drug response. The associations are primarily consistent across the two RNA-Seq datasets, with some modest variation, indicating the importance of obtaining multiple datasets to produce a robust resource. We further used an empirical Bayesian model to compare eQTL patterns in liver and an additional 20 GTEx tissues, finding that MHC genes, and especially class II genes, are enriched for liver-specific eQTL patterns. To illustrate the utility of the resource to augment GWAS analysis with small sample sizes, we developed a novel meta-analysis technique to combine several liver eQTL data sources. We also illustrate its application using a transcriptome-enhanced re-analysis of a study of neutropenia in pancreatic cancer patients. The associations of genotype with liver expression, including splice variation and its genetic associations, are made available in a searchable genome browser.

Characterization of population variability of 1,3-butadiene derived protein adducts in humans and mice.

1,3-butadiene is a known human carcinogen and a chemical to which humans are exposed occupationally and through environmental pollution. Inhalation risk assessment of 1,3-butadiene was completed several decades ago before data on molecular biomarkers of exposure and effect have been reported from both human studies of workers and experimental studies in mice. To improve risk assessment of 1,3-butadiene, the quantitative characterization of uncertainty in estimations of inter-individual variability in cancer-related effects is needed. For this, we ought to take advantage of the availability of the data on 1,3-butadiene hemoglobin adducts, well established biomarkers of the internal dose of the reactive epoxides, from several large-scale human studies and from a study in a Collaborative Cross mouse population. We found that in humans, toxicokinetic uncertainty factor for 99th percentile of the population ranged from 3.27 to 7.9, depending on the hemoglobin adduct. For mice, these values ranged from less than 2 to 7.51, depending on the dose and the adduct. Quantitative estimated from this study can be used to reduce uncertainties in the parameter estimates used in the models to derive the inhalation unit risk, as well as to address possible differences in variability in 1,3-butadiene metabolism that may be dose-related.

Leveraging TOPMed imputation server and constructing a cohort-specific imputation reference panel to enhance genotype imputation among cystic fibrosis patients.

Cystic fibrosis (CF) is a severe genetic disorder that can cause multiple comorbidities affecting the lungs, the pancreas, the luminal digestive system and beyond. In our previous genome-wide association studies (GWAS), we genotyped approximately 8,000 CF samples using a mixture of different genotyping platforms. More recently, the Cystic Fibrosis Genome Project (CFGP) performed deep (approximately 30×) whole genome sequencing (WGS) of 5,095 samples to better understand the genetic mechanisms underlying clinical heterogeneity among patients with CF. For mixtures of GWAS array and WGS data, genotype imputation has proven effective in increasing effective sample size. Therefore, we first performed imputation for the approximately 8,000 CF samples with GWAS array genotype using the Trans-Omics for Precision Medicine (TOPMed) freeze 8 reference panel. Our results demonstrate that TOPMed can provide high-quality imputation for patients with CF, boosting genomic coverage from approximately 0.3-4.2 million genotyped markers to approximately 11-43 million well-imputed markers, and significantly improving polygenic risk score (PRS) prediction accuracy. Furthermore, we built a CF-specific CFGP reference panel based on WGS data of patients with CF. We demonstrate that despite having approximately 3% the sample size of TOPMed, our CFGP reference panel can still outperform TOPMed when imputing some CF disease-causing variants, likely owing to allele and haplotype differences between patients with CF and general populations. We anticipate our imputed data for 4,656 samples without WGS data will benefit our subsequent genetic association studies, and the CFGP reference panel built from CF WGS samples will benefit other investigators studying CF.

Intra- and Inter-Species Variability in Urinary N7-(1-Hydroxy-3-buten-2-yl)guanine Adducts Following Inhalation Exposure to 1,3-Butadiene.

1,3-Butadiene is a known carcinogen primarily targeting lymphoid tissues, lung, and liver. Cytochrome P450 activates butadiene to epoxides which form covalent DNA adducts that are thought to be a key mechanistic event in cancer. Previous studies suggested that inter-species, -tissue, and -individual susceptibility to adverse health effects of butadiene exposure may be due to differences in metabolism and other mechanisms. In this study, we aimed to examine the extent of inter-individual and inter-species variability in the urinary N7-(1-hydroxy-3-buten-2-yl)guanine (EB-GII) DNA adduct, a well-known biomarker of exposure to butadiene. For a population variability study in mice, we used the collaborative cross model. Female and male mice from five strains were exposed to filtered air or butadiene (590 ppm, 6 h/day, 5 days/week for 2 weeks) by inhalation. Urine samples were collected, and the metabolic activation of butadiene by DNA-reactive species was quantified as urinary EB-GII adducts. We quantified the degree of EB-GII variation across mouse strains and sexes; then, we compared this variation with the data from rats (exposed to 62.5 or 200 ppm butadiene) and humans (0.004-2.2 ppm butadiene). We show that sex and strain are significant contributors to the variability in urinary EB-GII levels in mice. In addition, we find that the degree of variability in urinary EB-GII in collaborative cross mice, when expressed as an uncertainty factor for the inter-individual variability (UFH), is relatively modest (≤threefold) possibly due to metabolic saturation. By contrast, the variability in urinary EB-GII (adjusted for exposure) observed in humans, while larger than the default value of 10-fold, is largely consistent with UFH estimates for other chemicals based on human data for non-cancer endpoints. Overall, these data demonstrate that urinary EB-GII levels, particularly from human studies, may be useful for quantitative characterization of human variability in cancer risks to butadiene.

Environmental impacts of Hurricane Florence flooding in eastern North Carolina: temporal analysis of contaminant distribution and potential human health risks.

BACKGROUND: Hurricane Florence made landfall in North Carolina in September 2018 causing extensive flooding. Several potential point sources of hazardous substances and Superfund sites sustained water damage and contaminants may have been released into the environment. OBJECTIVE: This study conducted temporal analysis of contaminant distribution and potential human health risks from Hurricane Florence-associated flooding. METHODS: Soil samples were collected from 12 sites across four counties in North Carolina in September 2018, January and May 2019. Chemical analyses were performed for organics by gas chromatography-mass spectrometry. Metals were analyzed using inductively coupled plasma mass spectrometry. Hazard index and cancer risk were calculated using EPA Regional Screening Level Soil Screening Levels for residential soils. RESULTS: PAH and metals detected downstream from the coal ash storage pond that leaked were detected and were indicative of a pyrogenic source of contamination. PAH at these sites were of human health concern because cancer risk values exceeded 1 × 10-6 threshold. Other contaminants measured across sampling sites, or corresponding hazard index and cancer risk, did not exhibit spatial or temporal differences or were of concern. SIGNIFICANCE: This work shows the importance of rapid exposure assessment following natural disasters. It also establishes baseline levels of contaminants for future comparisons.

Risk Characterization of Environmental Samples Using In Vitro Bioactivity and Polycyclic Aromatic Hydrocarbon Concentrations Data.

Methods to assess environmental exposure to hazardous chemicals have primarily focused on quantification of individual chemicals, although chemicals often occur in mixtures, presenting challenges to the traditional risk characterization framework. Sampling sites in a defined geographic region provide an opportunity to characterize chemical contaminants, with spatial interpolation as a tool to provide estimates for non-sampled sites. At the same time, the use of in vitro bioactivity measurements has been shown to be informative for rapid risk-based decisions. In this study, we measured in vitro bioactivity in 39 surface soil samples collected immediately after flooding associated with Hurricane Harvey in Texas in a residential area known to be inundated with polycyclic aromatic hydrocarbon (PAH) contaminants. Bioactivity data were from a number of functional and toxicity assays in 5 human cell types, such as induced pluripotent stem cell-derived hepatocytes, cardiomyocytes, neurons, and endothelial cells, as well as human umbilical vein endothelial cells. Data on concentrations of PAH in these samples were also available and the combination of data sources offered a unique opportunity to assess the joint spatial variation of PAH components and bioactivity. We found significant evidence of spatial correlation of a subset of PAH contaminants and of cell-based phenotypes. In addition, we show that the cell-based bioactivity data can be used to predict environmental concentrations for several PAH contaminants, as well as overall PAH summaries and cancer risk. This study's impact lies in its demonstration that cell-based profiling can be used for rapid hazard screening of environmental samples by anchoring the bioassays to concentrations of PAH. This work sets the stage for identification of the areas of concern and direct quantitative risk characterization based on bioactivity data, thereby providing an important supplement to traditional individual chemical analyses by shedding light on constituents that may be missed from targeted chemical monitoring.

Classification of estrogenic compounds by coupling high content analysis and machine learning algorithms.

Environmental toxicants affect human health in various ways. Of the thousands of chemicals present in the environment, those with adverse effects on the endocrine system are referred to as endocrine-disrupting chemicals (EDCs). Here, we focused on a subclass of EDCs that impacts the estrogen receptor (ER), a pivotal transcriptional regulator in health and disease. Estrogenic activity of compounds can be measured by many in vitro or cell-based high throughput assays that record various endpoints from large pools of cells, and increasingly at the single-cell level. To simultaneously capture multiple mechanistic ER endpoints in individual cells that are affected by EDCs, we previously developed a sensitive high throughput/high content imaging assay that is based upon a stable cell line harboring a visible multicopy ER responsive transcription unit and expressing a green fluorescent protein (GFP) fusion of ER. High content analysis generates voluminous multiplex data comprised of minable features that describe numerous mechanistic endpoints. In this study, we present a machine learning pipeline for rapid, accurate, and sensitive assessment of the endocrine-disrupting potential of benchmark chemicals based on data generated from high content analysis. The multidimensional imaging data was used to train a classification model to ultimately predict the impact of unknown compounds on the ER, either as agonists or antagonists. To this end, both linear logistic regression and nonlinear Random Forest classifiers were benchmarked and evaluated for predicting the estrogenic activity of unknown compounds. Furthermore, through feature selection, data visualization, and model discrimination, the most informative features were identified for the classification of ER agonists/antagonists. The results of this data-driven study showed that highly accurate and generalized classification models with a minimum number of features can be constructed without loss of generality, where these machine learning models serve as a means for rapid mechanistic/phenotypic evaluation of the estrogenic potential of many chemicals.

A New Liver Expression Quantitative Trait Locus Map From 1,183 Individuals Provides Evidence for Novel Expression Quantitative Trait Loci of Drug Response, Metabolic, and Sex-Biased Phenotypes.

Expression quantitative trait locus (eQTL) studies in human liver are crucial for elucidating how genetic variation influences variability in disease risk and therapeutic outcomes and may help guide strategies to obtain maximal efficacy and safety of clinical interventions. Associations between expression microarray and genome-wide genotype data from four human liver eQTL studies (n = 1,183) were analyzed. More than 2.3 million cis-eQTLs for 15,668 genes were identified. When eQTLs were filtered against a list of 1,496 drug response genes, 187,829 cis-eQTLs for 1,191 genes were identified. Additionally, 1,683 sex-biased cis-eQTLs were identified, as well as 49 and 73 cis-eQTLs that colocalized with genome-wide association study signals for blood metabolite or lipid levels, respectively. Translational relevance of these results is evidenced by linking DPYD eQTLs to differences in safety of chemotherapy, linking the sex-biased regulation of PCSK9 expression to anti-lipid therapy, and identifying the G-protein coupled receptor GPR180 as a novel drug target for hypertriglyceridemia.

Using Collaborative Cross Mouse Population to Fill Data Gaps in Risk Assessment: A Case Study of Population-Based Analysis of Toxicokinetics and Kidney Toxicodynamics of Tetrachloroethylene.

BACKGROUND: Interindividual variability in susceptibility remains poorly characterized for environmental chemicals such as tetrachloroethylene (PERC). Development of population-based experimental models provide a potential approach to fill this critical need in human health risk assessment. OBJECTIVES: In this study, we aimed to better characterize the contribution of glutathione (GSH) conjugation to kidney toxicity of PERC and the degree of associated interindividual toxicokinetic (TK) and toxicodynamic (TD) variability by using the Collaborative Cross (CC) mouse population. METHODS: Male mice from 45 strains were intragastrically dosed with PERC ([Formula: see text]) or vehicle (5% Alkamuls EL-620 in saline), and time-course samples were collected for up to 24 h. Population variability in TK of S-(1,2,2-trichlorovinyl)GSH (TCVG), S-(1,2,2-trichlorovinyl)-L-cysteine (TCVC), and N-acetyl-S-(1,2,2-trichlorovinyl)-L-cysteine (NAcTCVC) was quantified in serum, liver, and kidney, and analyzed using a toxicokinetic model. Effects of PERC on kidney weight, fatty acid metabolism-associated genes [ Acot1 (Acyl-CoA thioesterase 1), Fabp1 (fatty acid-binding protein 1), and Ehhadh (enoyl-coenzyme A, hydratase/3-hydroxyacyl coenzyme A dehydrogenase)], and a marker of proximal tubular injury [KIM-1 (kidney injury molecule-1)/Hepatitis A virus cellular receptor 1 ( Havcr1)] were evaluated. Finally, quantitative data on interstrain variability in both formation of GSH conjugation metabolites of PERC and its kidney effects was used to calculate adjustment factors for the interindividual variability in both TK and TD. RESULTS: Mice treated with PERC had significantly lower kidney weight, higher kidney-to-body weight (BW) ratio, and higher expression of fatty acid metabolism-associated genes ( Acot1, Fabp1, and Ehhadh) and a marker of proximal tubular injury (KIM-1/ Havcr1). Liver levels of TCVG were significantly correlated with KIM-1/ Havcr1 in kidney, consistent with kidney injury being associated with GSH conjugation. We found that the default uncertainty factor for human variability may be marginally adequate to protect 95%, but not more, of the population for kidney toxicity mediated by PERC. DISCUSSION: Overall, this study demonstrates the utility of the CC mouse population in characterizing metabolism-toxicity interactions and quantifying interindividual variability. Further refinement of the characterization of interindividual variability can be accomplished by incorporating these data into in silico population models both for TK (such as a physiologically based pharmacokinetic model), as well as for toxicodynamic responses. https://doi.org/10.1289/EHP5105.

Population-Based Analysis of DNA Damage and Epigenetic Effects of 1,3-Butadiene in the Mouse.

Metabolism of 1,3-butadiene, a known human and rodent carcinogen, results in formation of reactive epoxides, a key event in its carcinogenicity. Although mice exposed to 1,3-butadiene present DNA adducts in all tested tissues, carcinogenicity is limited to liver, lung, and lymphoid tissues. Previous studies demonstrated that strain- and tissue-specific epigenetic effects in response to 1,3-butadiene exposure may influence susceptibly to DNA damage and serve as a potential mechanism of tissue-specific carcinogenicity. This study aimed to investigate interindividual variability in the effects of 1,3-butadiene using a population-based mouse model. Male mice from 20 Collaborative Cross strains were exposed to 0 or 635 ppm 1,3-butadiene by inhalation (6 h/day, 5 days/week) for 2 weeks. We evaluated DNA damage and epigenetic effects in target (lung and liver) and nontarget (kidney) tissues of 1,3-butadiene-induced carcinogenesis. DNA damage was assessed by measuring N-7-(2,3,4-trihydroxybut-1-yl)-guanine (THB-Gua) adducts. To investigate global histone modification alterations, we evaluated the trimethylation and acetylation of histones H3 and H4 across tissues. Changes in global cytosine DNA methylation were evaluated from the levels of methylation of LINE-1 and SINE B1 retrotransposons. We quantified the degree of variation across strains, deriving a chemical-specific human variability factor to address population variability in carcinogenic risk, which is largely ignored in current cancer risk assessment practice. Quantitative trait locus mapping identified four candidate genes related to chromatin remodeling whose variation was associated with interstrain susceptibility. Overall, this study uses 1,3-butadiene to demonstrate how the Collaborative Cross mouse population can be used to identify the mechanisms for and quantify the degree of interindividual variability in tissue-specific effects that are relevant to chemically induced carcinogenesis.

Polyphasic characterization of four soil-derived phenanthrene-degrading Acidovorax strains and proposal of Acidovorax carolinensis sp. nov.

Four bacterial strains identified as members of the Acidovorax genus were isolated from two geographically distinct but similarly contaminated soils in North Carolina, USA, characterized, and their genomes sequenced. Their 16S rRNA genes were highly similar to those previously recovered during stable-isotope probing (SIP) of one of the soils with the polycyclic aromatic hydrocarbon (PAH) phenanthrene. Heterotrophic growth of all strains occurred with a number of organic acids, as well as phenanthrene, but no other tested PAHs. Optimal growth occurred aerobically under mesophilic temperature, neutral pH, and low salinity conditions. Predominant fatty acids were C16:1ω7c/C16:1ω6c, C16:0, and C18:1ω7c, and were consistent with the genus. Genomic G+C contents ranged from 63.6 to 64.2%. A combination of whole genome comparisons and physiological analyses indicated that these four strains likely represent a single species within the Acidovorax genus. Chromosomal genes for phenanthrene degradation to phthalate were nearly identical to highly conserved regions in phenanthrene-degrading Delftia, Burkholderia, Alcaligenes, and Massilia species in regions flanked by transposable or extrachromosomal elements. The lower degradation pathway for phenanthrene metabolism was inferred by comparisons to described genes and proteins. The novel species Acidovorax carolinensis sp. nov. is proposed, comprising the four strains described in this study with strain NA3T as the type strain (=LMG 30136, =DSM 105008).

Correction: Molecular Subtypes in Head and Neck Cancer Exhibit Distinct Patterns of Chromosomal Gain and Loss of Canonical Cancer Genes.

[This corrects the article DOI: 10.1371/journal.pone.0056823.].