Two-Stage Machine Learning-Based Approach to Predict Points of Departure for Human Noncancer and Developmental/Reproductive Effects.

Chemical points of departure (PODs) for critical health effects are crucial for evaluating and managing human health risks and impacts from exposure. However, PODs are unavailable for most chemicals in commerce due to a lack of in vivo toxicity data. We therefore developed a two-stage machine learning (ML) framework to predict human-equivalent PODs for oral exposure to organic chemicals based on chemical structure. Utilizing ML-based predictions for structural/physical/chemical/toxicological properties from OPERA 2.9 as features (Stage 1), ML models using random forest regression were trained with human-equivalent PODs derived from in vivo data sets for general noncancer effects (n = 1,791) and reproductive/developmental effects (n = 2,228), with robust cross-validation for feature selection and estimating generalization errors (Stage 2). These two-stage models accurately predicted PODs for both effect categories with cross-validation-based root-mean-squared errors less than an order of magnitude. We then applied one or both models to 34,046 chemicals expected to be in the environment, revealing several thousand chemicals of moderate concern and several hundred chemicals of high concern for health effects at estimated median population exposure levels. Further application can expand by orders of magnitude the coverage of organic chemicals that can be evaluated for their human health risks and impacts.

Control of false discoveries in grouped hypothesis testing for eQTL data.

BACKGROUND: Expression quantitative trait locus (eQTL) analysis aims to detect the genetic variants that influence the expression of one or more genes. Gene-level eQTL testing forms a natural grouped-hypothesis testing strategy with clear biological importance. Methods to control family-wise error rate or false discovery rate for group testing have been proposed earlier, but may not be powerful or easily apply to eQTL data, for which certain structured alternatives may be defensible and may enable the researcher to avoid overly conservative approaches. RESULTS: In an empirical Bayesian setting, we propose a new method to control the false discovery rate (FDR) for grouped hypotheses. Here, each gene forms a group, with SNPs annotated to the gene corresponding to individual hypotheses. The heterogeneity of effect sizes in different groups is considered by the introduction of a random effects component. Our method, entitled Random Effects model and testing procedure for Group-level FDR control (REG-FDR), assumes a model for alternative hypotheses for the eQTL data and controls the FDR by adaptive thresholding. As a convenient alternate approach, we also propose Z-REG-FDR, an approximate version of REG-FDR, that uses only Z-statistics of association between genotype and expression for each gene-SNP pair. The performance of Z-REG-FDR is evaluated using both simulated and real data. Simulations demonstrate that Z-REG-FDR performs similarly to REG-FDR, but with much improved computational speed. CONCLUSION: Our results demonstrate that the Z-REG-FDR method performs favorably compared to other methods in terms of statistical power and control of FDR. It can be of great practical use for grouped hypothesis testing for eQTL analysis or similar problems in statistical genomics due to its fast computation and ability to be fit using only summary data.

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.

Risk-based prioritization of PFAS using phenotypic and transcriptomic data from human induced pluripotent stem cell-derived hepatocytes and cardiomyocytes.

Per- and polyfluoroalkyl substances (PFAS) are chemicals with important applications; they are persistent in the environment and may pose human health hazards. Regulatory agencies are considering restrictions and bans of PFAS; however, little data exists for informed decisions. Several prioritization strategies were proposed for evaluation of potential hazards of PFAS. Structure-based grouping could expedite the selection of PFAS for testing; still, the hypothesis that structure-effect relationships exist for PFAS requires confirmation. We tested 26 structurally diverse PFAS from 8 groups using human-induced pluripotent stem cell-derived hepatocytes and cardiomyocytes, and tested concentration-response effects on cell function and gene expression. Few phenotypic effects were observed in hepatocytes, but negative chronotropy was observed for 8 of the 26 PFAS. Substance- and cell type-dependent transcriptomic changes were more prominent but lacked substantial group-specific effects. In hepatocytes, we found up-regulation of stress-related and extracellular matrix organization pathways, and down-regulation of fat metabolism. In cardiomyocytes, contractility-related pathways were most affected. We derived phenotypic and transcriptomic points of departure and compared them to predicted PFAS exposures. The conservative estimates for bioactivity and exposure were used to derive bioactivity-to-exposure ratio (BER) for each PFAS, most (23 of 26) PFAS had BER > 1. Overall, these data suggests that structure-based grouping of PFAS may not be sufficient to predict their biological effects. Testing of individual PFAS may be needed for scientific-based decision-making. Our proposed strategy of using two human cell types and considering phenotypic and transcriptomic effects, combined with dose-response analysis and calculation of BER, may be used for PFAS prioritization.

Per- and polyfluoroalkyl substances (PFAS) are man-made chemicals used in many products. However, most of these substances have not been tested for safety, and concerns exist that they may be harmful to human health and/or the environment. This study aimed to use human cell-based models to investigate if some of the PFAS may exhibit hazardous properties and if similarities among substances are observed. Few effects were observed in liver cells, but a decrease in beating frequency was observed in heart cells for some PFAS. Gene expression changes were substance- and cell type-dependent. We did not find convincing structure-based similarities among PFAS; this suggests that testing of individual PFAS may be necessary in the future to inform health decisions. Overall, this study showed that a test strategy of using two human cell types, from liver and heart, may inform PFAS prioritization without a need for testing in animals.

Hazard and risk characterization of 56 structurally diverse PFAS using a targeted battery of broad coverage assays using six human cell types.

Per- and poly-fluoroalkyl substances (PFAS) are extensively used in commerce leading to their prevalence in the environment. Due to their chemical stability, PFAS are considered to be persistent and bioaccumulative; they are frequently detected in both the environment and humans. Because of this, PFAS as a class (composed of hundreds to thousands of chemicals) are contaminants of very high concern. Little information is available for the vast majority of PFAS, and regulatory agencies lack safety data to determine whether exposure limits or restrictions are needed. Cell-based assays are a pragmatic approach to inform decision-makers on potential health hazards; therefore, we hypothesized that a targeted battery of human in vitro assays can be used to determine whether there are structure-bioactivity relationships for PFAS, and to characterize potential risks by comparing bioactivity (points of departure) to exposure estimates. We tested 56 PFAS from 8 structure-based subclasses in concentration response (0.1-100 µM) using six human cell types selected from target organs with suggested adverse effects of PFAS - human induced pluripotent stem cell (iPSC)-derived hepatocytes, neurons, and cardiomyocytes, primary human hepatocytes, endothelial and HepG2 cells. While many compounds were without effect; certain PFAS demonstrated cell-specific activity highlighting the necessity of using a compendium of in vitro models to identify potential hazards. No class-specific groupings were evident except for some chain length- and structure-related trends. In addition, margins of exposure (MOE) were derived using empirical and predicted exposure data. Conservative MOE calculations showed that most tested PFAS had a MOE in the 1-100 range; ∼20% of PFAS had MOE<1, providing tiered priorities for further studies. Overall, we show that a compendium of human cell-based models can be used to derive bioactivity estimates for a range of PFAS, enabling comparisons with human biomonitoring data. Furthermore, we emphasize that establishing structure-bioactivity relationships may be challenging for the tested PFAS.

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.

Analysis of reproducibility and robustness of a renal proximal tubule microphysiological system OrganoPlate 3-lane 40 for in vitro studies of drug transport and toxicity.

Microphysiological systems are an emerging area of in vitro drug development, and their independent evaluation is important for wide adoption and use. The primary goal of this study was to test reproducibility and robustness of a renal proximal tubule microphysiological system, OrganoPlate 3-lane 40, as an in vitro model for drug transport and toxicity studies. This microfluidic model was compared with static multiwell cultures and tested using several human renal proximal tubule epithelial cell (RPTEC) types. The model was characterized in terms of the functional transport for various tubule-specific proteins, epithelial permeability of small molecules (cisplatin, tenofovir, and perfluorooctanoic acid) versus large molecules (fluorescent dextrans, 60-150 kDa), and gene expression response to a nephrotoxic xenobiotic. The advantages offered by OrganoPlate 3-lane 40 as compared with multiwell cultures are the presence of media flow, albeit intermittent, and increased throughput compared with other microfluidic models. However, OrganoPlate 3-lane 40 model appeared to offer only limited (eg, MRP-mediated transport) advantages in terms of either gene expression or functional transport when compared with the multiwell plate culture conditions. Although OrganoPlate 3-lane 40 can be used to study cellular uptake and direct toxic effects of small molecules, it may have limited utility for drug transport studies. Overall, this study offers refined experimental protocols and comprehensive comparative data on the function of RPETCs in traditional multiwell culture and microfluidic OrganoPlate 3-lane 40, information that will be invaluable for the prospective end-users of in vitro models of the human proximal tubule.

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.

Integrative Chemical-Biological Grouping of Complex High Production Volume Substances from Lower Olefin Manufacturing Streams.

Human cell-based test methods can be used to evaluate potential hazards of mixtures and products of petroleum refining ("unknown or variable composition, complex reaction products, or biological materials" substances, UVCBs). Analyses of bioactivity and detailed chemical characterization of petroleum UVCBs were used separately for grouping these substances; a combination of the approaches has not been undertaken. Therefore, we used a case example of representative high production volume categories of petroleum UVCBs, 25 lower olefin substances from low benzene naphtha and resin oils categories, to determine whether existing manufacturing-based category grouping can be supported. We collected two types of data: nontarget ion mobility spectrometry-mass spectrometry of both neat substances and their organic extracts and in vitro bioactivity of the organic extracts in five human cell types: umbilical vein endothelial cells and induced pluripotent stem cell-derived hepatocytes, endothelial cells, neurons, and cardiomyocytes. We found that while similarity in composition and bioactivity can be observed for some substances, existing categories are largely heterogeneous. Strong relationships between composition and bioactivity were observed, and individual constituents that determine these associations were identified. Overall, this study showed a promising approach that combines chemical composition and bioactivity data to better characterize the variability within manufacturing categories of petroleum UVCBs.

Simultaneous modeling of multivariate heterogeneous responses and heteroskedasticity via a two-stage composite likelihood.

Multivariate heterogeneous responses and heteroskedasticity have attracted increasing attention in recent years. In genome-wide association studies, effective simultaneous modeling of multiple phenotypes would improve statistical power and interpretability. However, a flexible common modeling system for heterogeneous data types can pose computational difficulties. Here we build upon a previous method for multivariate probit estimation using a two-stage composite likelihood that exhibits favorable computational time while retaining attractive parameter estimation properties. We extend this approach to incorporate multivariate responses of heterogeneous data types (binary and continuous), and possible heteroskedasticity. Although the approach has wide applications, it would be particularly useful for genomics, precision medicine, or individual biomedical prediction. Using a genomics example, we explore statistical power and confirm that the approach performs well for hypothesis testing and coverage percentages under a wide variety of settings. The approach has the potential to better leverage genomics data and provide interpretable inference for pleiotropy, in which a locus is associated with multiple traits.

A tiered testing strategy based on in vitro phenotypic and transcriptomic data for selecting representative petroleum UVCBs for toxicity evaluation in vivo.

Hazard evaluation of substances of "unknown or variable composition, complex reaction products and biological materials" (UVCBs) remains a major challenge in regulatory science because their chemical composition is difficult to ascertain. Petroleum substances are representative UVCBs and human cell-based data have been previously used to substantiate their groupings for regulatory submissions. We hypothesized that a combination of phenotypic and transcriptomic data could be integrated to make decisions as to selection of group-representative worst-case petroleum UVCBs for subsequent toxicity evaluation in vivo. We used data obtained from 141 substances from 16 manufacturing categories previously tested in 6 human cell types (induced pluripotent stem cell [iPSC]-derived hepatocytes, cardiomyocytes, neurons, and endothelial cells, and MCF7 and A375 cell lines). Benchmark doses for gene-substance combinations were calculated, and both transcriptomic and phenotype-derived points of departure (PODs) were obtained. Correlation analysis and machine learning were used to assess associations between phenotypic and transcriptional PODs and to determine the most informative cell types and assays, thus representing a cost-effective integrated testing strategy. We found that 2 cell types-iPSC-derived-hepatocytes and -cardiomyocytes-contributed the most informative and protective PODs and may be used to inform selection of representative petroleum UVCBs for further toxicity evaluation in vivo. Overall, although the use of new approach methodologies to prioritize UVCBs has not been widely adopted, our study proposes a tiered testing strategy based on iPSC-derived hepatocytes and cardiomyocytes to inform selection of representative worst-case petroleum UVCBs from each manufacturing category for further toxicity evaluation in vivo.

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.

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.

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.

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.

A Population-Based Human In Vitro Approach to Quantify Inter-Individual Variability in Responses to Chemical Mixtures.

Human cell-based population-wide in vitro models have been proposed as a strategy to derive chemical-specific estimates of inter-individual variability; however, the utility of this approach has not yet been tested for cumulative exposures in mixtures. This study aimed to test defined mixtures and their individual components and determine whether adverse effects of the mixtures were likely to be more variable in a population than those of the individual chemicals. The in vitro model comprised 146 human lymphoblastoid cell lines from four diverse subpopulations of European and African descent. Cells were exposed, in concentration−response, to 42 chemicals from diverse classes of environmental pollutants; in addition, eight defined mixtures were prepared from these chemicals using several exposure- or hazard-based scenarios. Points of departure for cytotoxicity were derived using Bayesian concentration−response modeling and population variability was quantified in the form of a toxicodynamic variability factor (TDVF). We found that 28 chemicals and all mixtures exhibited concentration−response cytotoxicity, enabling calculation of the TDVF. The median TDVF across test substances, for both individual chemicals or defined mixtures, ranged from a default assumption (101/2) of toxicodynamic variability in human population to >10. The data also provide a proof of principle for single-variant genome-wide association mapping for toxicity of the chemicals and mixtures, although replication would be necessary due to statistical power limitations with the current sample size. This study demonstrates the feasibility of using a set of human lymphoblastoid cell lines as an in vitro model to quantify the extent of inter-individual variability in hazardous properties of both individual chemicals and mixtures. The data show that population variability of the mixtures is unlikely to exceed that of the most variable component, and that similarity in genome-wide associations among components may be used to accrue additional evidence for grouping of constituents in a mixture for cumulative assessments.

TwinEQTL: ultrafast and powerful association analysis for eQTL and GWAS in twin studies.

We develop a computationally efficient alternative, TwinEQTL, to a linear mixed-effects model for twin genome-wide association study data. Instead of analyzing all twin samples together with linear mixed-effects model, TwinEQTL first splits twin samples into 2 independent groups on which multiple linear regression analysis can be validly performed separately, followed by an appropriate meta-analysis-like approach to combine the 2 nonindependent test results. Through mathematical derivations, we prove the validity of TwinEQTL algorithm and show that the correlation between 2 dependent test statistics at each single-nucleotide polymorphism is independent of its minor allele frequency. Thus, the correlation is constant across all single-nucleotide polymorphisms. Through simulations, we show empirically that TwinEQTL has well controlled type I error with negligible power loss compared with the gold-standard linear mixed-effects models. To accommodate expression quantitative loci analysis with twin subjects, we further implement TwinEQTL into an R package with much improved computational efficiency. Our approaches provide a significant leap in terms of computing speed for genome-wide association study and expression quantitative loci analysis with twin samples.

Model systems and organisms for addressing inter- and intra-species variability in risk assessment.

Addressing inter- and intra-species differences in potential hazardous effects of chemicals remains a long-standing challenge in human health risk assessment that is typically addressed heuristically through use of 10-fold default "uncertainty" or "safety" factors. Although it has long been recognized that chemical-specific data would be preferable to replace the "defaults," only recently have there emerged experimental model systems and organisms with the potential to experimentally quantify the population variability in both toxicokinetics and toxicodynamics for specific chemicals. Progress is most evident in the use of population in vitro human cell-based models and population in vivo mouse models. Multiple case studies were published in the past 10-15 years that clearly demonstrate the utility of such models to derive data with direct application to quantifying variability at hazard identification, exposure-response assessment, and mechanistic understanding of toxicity steps of traditional risk assessments. Here, we review recent efforts to develop fit-for-purpose approaches utilizing these novel population-based in vitro and in vivo models in the context of risk assessment. We also describe key challenges and opportunities to broadening application of population-based experimental approaches. We conclude that population-based models are now beginning to realize their potential to address long-standing data gaps in inter- and intra-species variability.

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.