Benchmark dose (BMD) modeling: current practice, issues, and challenges.

Benchmark dose (BMD) modeling is now the state of the science for determining the point of departure for risk assessment. Key advantages include the fact that the modeling takes account of all of the data for a particular effect from a particular experiment, increased consistency, and better accounting for statistical uncertainties. Despite these strong advantages, disagreements remain as to several specific aspects of the modeling, including differences in the recommendations of the US Environmental Protection Agency (US EPA) and the European Food Safety Authority (EFSA). Differences exist in the choice of the benchmark response (BMR) for continuous data, the use of unrestricted models, and the mathematical models used; these can lead to differences in the final BMDL. It is important to take confidence in the model into account in choosing the BMDL, rather than simply choosing the lowest value. The field is moving in the direction of model averaging, which will avoid many of the challenges of choosing a single best model when the underlying biology does not suggest one, but additional research would be useful into methods of incorporating biological considerations into the weights used in the averaging. Additional research is also needed regarding the interplay between the BMR and the UF to ensure appropriate use for studies supporting a lower BMR than default values, such as for epidemiology data. Addressing these issues will aid in harmonizing methods and moving the field of risk assessment forward.

Human relevance of rodent liver tumors: Key insights from a Toxicology Forum workshop on nongenotoxic modes of action.

The Toxicology Forum sponsored a workshop in October 2016, on the human relevance of rodent liver tumors occurring via nongenotoxic modes of action (MOAs). The workshop focused on two nuclear receptor-mediated MOAs (Constitutive Androstane Receptor (CAR) and Peroxisome Proliferator Activated Receptor-alpha (PPARα), and on cytotoxicity. The goal of the meeting was to review the state of the science to (1) identify areas of consensus and differences, data gaps and research needs; (2) identify reasons for inconsistencies in current regulatory positions; and (3) consider what data are needed to demonstrate a specific MOA, and when additional research is needed to rule out alternative possibilities. Implications for quantitative risk assessment approaches were discussed, as were implications of not considering MOA and dose in hazard characterization and labeling schemes. Most, but not all, participants considered the CAR and PPARα MOAs as not relevant to humans based on quantitative and qualitative differences. In contrast, cytotoxicity is clearly relevant to humans, but a threshold applies. Questions remain for all three MOAs concerning what data are necessary to determine the MOA and to what extent it is necessary to exclude other MOAs.

Derivation of an oral toxicity reference value for nickel.

Nickel (Ni) is in the earth's crust and can be found in environmental compartments such as water, soil, and air, as well as food. This paper presents an assessment of the oral nickel toxicity data in support of non-cancer health-based oral exposure limits or toxicity reference values (TRVs). This paper derives TRVs for three populations of interest: adults, toddlers, and people who have been dermally sensitized to nickel. The adult/lifetime TRV of 20 µg Ni/kg-day is based on post-implantation loss/perinatal mortality in a 2-generation reproductive study in rats. Several recent assessments by regulatory agencies have used the same study and endpoint, but the dose-response modeling conducted here was more appropriate for the study design. Toxicokinetic data from rats and humans indicate that the applied uncertainty factors are very conservative. Because the endpoint relates to fetal exposure and is not relevant to toddlers, a toddler TRV was derived based on decreased body weight in young rats; this TRV was also 20 µg Ni/kg-day. A separate TRV of 4 µg Ni/kg in addition to Ni in food was derived for protection of nickel-sensitized populations from flare-up of dermatitis, based on studies of single exposures in humans under conditions that maximize oral absorption.

Framework for human health risk assessment of non-cancer effects resulting from short-duration and intermittent exposures to chemicals.

Durations of exposure to chemicals, whether for single, repeated or intermittent periods, may vary from those upon which most guidance values are normally based. Because it is presently not feasible to conduct toxicity studies or develop toxicity reference values (TRVs) specific to each scenario of interest, methods are needed to address these various durations, drawing as much as possible on existing TRVs. A working framework was developed to address the potential for non-cancer effects resulting from continuous short-duration and intermittent exposures to chemicals. The framework presents an integrated, tiered approach that assists the user in identifying when existing TRVs can be applied directly, and the adaptations needed to assess the acceptability of short-duration or intermittent exposure scenarios. Descriptions of when and how toxicokinetic and toxicodynamic aspects need to be taken into consideration are also presented. The framework incorporates the use of TRVs based on exposure periods as similar as possible to the "actual" exposure periods and application of dose averaging under limited, specified conditions. This framework has been developed to aid in improving the scientific basis for the evaluation of short-duration and intermittent exposures in a variety of settings. Copyright © 2016 John Wiley & Sons, Ltd.

Advancing Risk Analysis for Nanoscale Materials: Report from an International Workshop on the Role of Alternative Testing Strategies for Advancement.

The Society for Risk Analysis (SRA) has a history of bringing thought leadership to topics of emerging risk. In September 2014, the SRA Emerging Nanoscale Materials Specialty Group convened an international workshop to examine the use of alternative testing strategies (ATS) for manufactured nanomaterials (NM) from a risk analysis perspective. Experts in NM environmental health and safety, human health, ecotoxicology, regulatory compliance, risk analysis, and ATS evaluated and discussed the state of the science for in vitro and other alternatives to traditional toxicology testing for NM. Based on this review, experts recommended immediate and near-term actions that would advance ATS use in NM risk assessment. Three focal areas-human health, ecological health, and exposure considerations-shaped deliberations about information needs, priorities, and the next steps required to increase confidence in and use of ATS in NM risk assessment. The deliberations revealed that ATS are now being used for screening, and that, in the near term, ATS could be developed for use in read-across or categorization decision making within certain regulatory frameworks. Participants recognized that leadership is required from within the scientific community to address basic challenges, including standardizing materials, protocols, techniques and reporting, and designing experiments relevant to real-world conditions, as well as coordination and sharing of large-scale collaborations and data. Experts agreed that it will be critical to include experimental parameters that can support the development of adverse outcome pathways. Numerous other insightful ideas for investment in ATS emerged throughout the discussions and are further highlighted in this article.

Bayesian evaluation of a physiologically-based pharmacokinetic (PBPK) model of long-term kinetics of metal nanoparticles in rats.

Biomathematical modeling quantitatively describes the disposition of metal nanoparticles in lungs and other organs of rats. In a preliminary model, adjustable parameters were calibrated to each of three data sets using a deterministic approach, with optimal values varying among the different data sets. In the current effort, Bayesian population analysis using Markov chain Monte Carlo (MCMC) simulation was used to recalibrate the model while improving assessments of parameter variability and uncertainty. The previously-developed model structure and some physiological parameter values were modified to improve physiological realism. The data from one of the three previously-identified studies and from two other studies were used for model calibration. The data from the one study that adequately characterized mass balance were used to generate parameter distributions. When data from a second study of the same nanomaterial (iridium) were added, the level of agreement was still acceptable. Addition of another data set (for silver nanoparticles) led to substantially lower precision in parameter estimates and large discrepancies between the model predictions and experimental data for silver nanoparticles. Additional toxicokinetic data are needed to further evaluate the model structure and performance and to reduce uncertainty in the kinetic processes governing in vivo disposition of metal nanoparticles.

Systems Biology and Biomarkers of Early Effects for Occupational Exposure Limit Setting.

In a recent National Research Council document, new strategies for risk assessment were described to enable more accurate and quicker assessments. This report suggested that evaluating individual responses through increased use of bio-monitoring could improve dose-response estimations. Identification of specific biomarkers may be useful for diagnostics or risk prediction as they have the potential to improve exposure assessments. This paper discusses systems biology, biomarkers of effect, and computational toxicology approaches and their relevance to the occupational exposure limit setting process. The systems biology approach evaluates the integration of biological processes and how disruption of these processes by chemicals or other hazards affects disease outcomes. This type of approach could provide information used in delineating the mode of action of the response or toxicity, and may be useful to define the low adverse and no adverse effect levels. Biomarkers of effect are changes measured in biological systems and are considered to be preclinical in nature. Advances in computational methods and experimental -omics methods that allow the simultaneous measurement of families of macromolecules such as DNA, RNA, and proteins in a single analysis have made these systems approaches feasible for broad application. The utility of the information for risk assessments from -omics approaches has shown promise and can provide information on mode of action and dose-response relationships. As these techniques evolve, estimation of internal dose and response biomarkers will be a critical test of these new technologies for application in risk assessment strategies. While proof of concept studies have been conducted that provide evidence of their value, challenges with standardization and harmonization still need to be overcome before these methods are used routinely.

Advancing human health risk assessment: integrating recent advisory committee recommendations.

Over the last dozen years, many national and international expert groups have considered specific improvements to risk assessment. Many of their stated recommendations are mutually supportive, but others appear conflicting, at least in an initial assessment. This review identifies areas of consensus and difference and recommends a practical, biology-centric course forward, which includes: (1) incorporating a clear problem formulation at the outset of the assessment with a level of complexity that is appropriate for informing the relevant risk management decision; (2) using toxicokinetics and toxicodynamic information to develop Chemical Specific Adjustment Factors (CSAF); (3) using mode of action (MOA) information and an understanding of the relevant biology as the key, central organizing principle for the risk assessment; (4) integrating MOA information into dose-response assessments using existing guidelines for non-cancer and cancer assessments; (5) using a tiered, iterative approach developed by the World Health Organization/International Programme on Chemical Safety (WHO/IPCS) as a scientifically robust, fit-for-purpose approach for risk assessment of combined exposures (chemical mixtures); and (6) applying all of this knowledge to enable interpretation of human biomonitoring data in a risk context. While scientifically based defaults will remain important and useful when data on CSAF or MOA to refine an assessment are absent or insufficient, assessments should always strive to use these data. The use of available 21st century knowledge of biological processes, clinical findings, chemical interactions, and dose-response at the molecular, cellular, organ and organism levels will minimize the need for extrapolation and reliance on default approaches.

Application of Markov chain Monte Carlo analysis to biomathematical modeling of respirable dust in US and UK coal miners.

A biomathematical model was previously developed to describe the long-term clearance and retention of particles in the lungs of coal miners. The model structure was evaluated and parameters were estimated in two data sets, one from the United States and one from the United Kingdom. The three-compartment model structure consists of deposition of inhaled particles in the alveolar region, competing processes of either clearance from the alveolar region or translocation to the lung interstitial region, and very slow, irreversible sequestration of interstitialized material in the lung-associated lymph nodes. Point estimates of model parameter values were estimated separately for the two data sets. In the current effort, Bayesian population analysis using Markov chain Monte Carlo simulation was used to recalibrate the model while improving assessments of parameter variability and uncertainty. When model parameters were calibrated simultaneously to the two data sets, agreement between the derived parameters for the two groups was very good, and the central tendency values were similar to those derived from the deterministic approach. These findings are relevant to the proposed update of the ICRP human respiratory tract model with revisions to the alveolar-interstitial region based on this long-term particle clearance and retention model.

A framework for fit-for-purpose dose response assessment.

The NRC report Science and Decisions: Advancing Risk Assessment made several recommendations to improve chemical risk assessment, with a focus on in-depth chronic dose-response assessments conducted by the U.S. Environmental Protection Agency. The recommendations addressed two broad elements: improving technical analysis and utility for decision making. To advance the discussions in the NRC report, in three multi-stakeholder workshops organized by the Alliance for Risk Assessment, available and evolving risk assessment methodologies were considered through the development and application of case studies. A key product was a framework (http://www.allianceforrisk.org/Workshop/Framework/ProblemFormulation.html) to guide risk assessors and managers to various dose-response assessment methods relevant to a range of decision contexts ranging from priority setting to full assessment, as illustrated by case studies. It is designed to facilitate selection of appropriate methodology for a variety of problem formulations and includes a variety of methods with supporting case studies, for areas flagged specifically by the NRC committee for consideration--e.g., susceptible sub-populations, population variability and background. The framewok contributes to organization and communication about methodologies for incorporating increasingly biologically informed and chemical specific knowledge into dose-response analysis, which is considered critical in evolving fit-for-purpose assessment to address relevant problem formulations.

Development of a framework for risk assessment of dietary carcinogens.

Although there are a number of guidance documents and frameworks for evaluation of carcinogenicity, none of the current methods fully reflects the state of the science. Common limitations include the absence of dose-response assessment and not considering the impact of differing exposure patterns (e.g., intermittent, high peaks vs. lower, continuous exposures). To address these issues, we have developed a framework for risk assessment of dietary carcinogens. This framework includes an enhanced approach for weight of evidence (WOE) evaluation for genetic toxicology data, with a focus on evaluating studies based on the most recent testing guidance to determine whether a chemical is a mutagen. Included alongside our framework is a discussion of resources for evaluating tissue dose and the temporal pattern of internal dose, taking into account the chemical's toxicokinetics. The framework then integrates the mode of action (MOA) and associated dose metric category with the exposure data to identify the appropriate approach(es) to low-dose extrapolation and level of concern associated with the exposure scenario. This framework provides risk managers with additional flexibility in risk management and risk communication options, beyond the binary choice of linear low-dose extrapolation vs. application of uncertainty factors.

Considerations for application of benchmark dose modeling in radiation research: workshop highlights.

BACKGROUND: Exposure to different forms of ionizing radiation occurs in diverse occupational, medical, and environmental settings. Improving the accuracy of the estimated health risks associated with exposure is therefore, essential for protecting the public, particularly as it relates to chronic low dose exposures. A key aspect to understanding health risks is precise and accurate modeling of the dose-response relationship. Toward this vision, benchmark dose (BMD) modeling may be a suitable approach for consideration in the radiation field. BMD modeling is already extensively used for chemical hazard assessments and is considered statistically preferable to identifying low and no observed adverse effects levels. BMD modeling involves fitting mathematical models to dose-response data for a relevant biological endpoint and identifying a point of departure (the BMD, or its lower bound). Recent examples in chemical toxicology show that when applied to molecular endpoints (e.g. genotoxic and transcriptional endpoints), BMDs correlate to points of departure for more apical endpoints such as phenotypic changes (e.g. adverse effects) of interest to regulatory decisions. This use of BMD modeling may be valuable to explore in the radiation field, specifically in combination with adverse outcome pathways, and may facilitate better interpretation of relevant in vivo and in vitro dose-response data. To advance this application, a workshop was organized on June 3rd, 2022, in Ottawa, Ontario that brought together BMD experts in chemical toxicology and the radiation scientific community of researchers, regulators, and policy-makers. The workshop's objective was to introduce radiation scientists to BMD modeling and its practical application using case examples from the chemical toxicity field and demonstrate the BMDExpress software using a radiation dataset. Discussions focused on the BMD approach, the importance of experimental design, regulatory applications, its use in supporting the development of adverse outcome pathways, and specific radiation-relevant examples. CONCLUSIONS: Although further deliberations are needed to advance the use of BMD modeling in the radiation field, these initial discussions and partnerships highlight some key steps to guide future undertakings related to new experimental work.

The Long Goodbye: Finally Moving on from the Relative Potency Approach to a Mixtures Approach for Polycyclic Aromatic Hydrocarbons (PAHs).

For the past several decades, a relative potency approach has been used to estimate the human health risks from exposure to polycyclic aromatic hydrocarbon (PAH) mixtures. Risk estimates are derived using potency equivalence factors (PEFs; also called relative potency factors [RPFs]), based on the ratio of selected PAHs to benzo[a]pyrene (BaP), expressed qualitatively by orders of magnitude. To quantify PEFs for 18 selected carcinogenic PAHs, a systematic approach with a priori and dose response criteria was developed, building on draft work by the US EPA in 2010 and its review by US EPA Science Advisory Board (SAB) in 2011. An exhaustive search for carcinogenicity studies that included both target PAHs and BaP with environmentally relevant exposure routes found only 48 animal bioassay datasets (mostly pre-1992 based on skin painting). Only eight datasets provided adequate low-response data, and of these only four datasets were appropriate for modeling to estimate PEFs; only benzo[b]fluoranthene and cyclopenta[c,d]pyrene had a PEF that could be quantified. Thus, current knowledge of PAH carcinogenicity is insufficient to support quantitative PEFs for PAH mixtures. This highlights the long-acknowledged need for an interdisciplinary approach to estimate risks from PAH mixtures. Use of alternative and short-term toxicity testing methods, improved mixture characterization, understanding the fate and bioavailability of PAH mixtures, and understanding exposure route-related differences in carcinogenicity are discussed as ways to improve the understanding of the risks of PAHs.

Linear low-dose extrapolation for noncancer heath effects is the exception, not the rule.

The nature of the exposure-response relationship has a profound influence on risk analyses. Several arguments have been proffered as to why all exposure-response relationships for both cancer and noncarcinogenic endpoints should be assumed to be linear at low doses. We focused on three arguments that have been put forth for noncarcinogens. First, the general "additivity-to-background" argument proposes that if an agent enhances an already existing disease-causing process, then even small exposures increase disease incidence in a linear manner. This only holds if it is related to a specific mode of action that has nonuniversal properties-properties that would not be expected for most noncancer effects. Second, the "heterogeneity in the population" argument states that variations in sensitivity among members of the target population tend to "flatten out and linearize" the exposure-response curve, but this actually only tends to broaden, not linearize, the dose-response relationship. Third, it has been argued that a review of epidemiological evidence shows linear or no-threshold effects at low exposures in humans, despite nonlinear exposure-response in the experimental dose range in animal testing for similar endpoints. It is more likely that this is attributable to exposure measurement error rather than a true nonthreshold association. Assuming that every chemical is toxic at high exposures and linear at low exposures does not comport to modern-day scientific knowledge of biology. There is no compelling evidence-based justification for a general low-exposure linearity; rather, case-specific mechanistic arguments are needed.

Impact of updated BMD modeling methods on perchlorate and chlorate assessments of human health hazard.

Several exposure limits for perchlorate have been developed based on an early key event, inhibition of radioactive iodide uptake (RAIU) by the thyroid. These assessments have used a variety of definitions of the point of departure. The current assessment revisited the modeling for inhibition of RAIU, using state of the science methods. Bayesian hierarchical modeling was used to account for the repeated measures on the same individuals in the key dataset, and the underlying Beta distribution used for the modeling correctly reflected the bounding of RAIU between 0 and 1. We defined the BMR as a point value of 8% RAIU (rather than a change in RAIU), based on descriptions in the medical literature that RAIU below this value is considered abnormal. Because a definition of the BMR based on the mean response would correspond to about 50% of the population with a response below the BMR at the benchmark dose, we used a hybrid definition of the BMR. That is, the BMD was defined as the dose at which it was estimated that there would be a 10% extra risk in the population of having RAIU of 8% or lower. The resulting point of departure based on the BMDL was 0.03 mg/kg-day.

Bayesian hierarchical evaluation of dose-response for peanut allergy in clinical trial screening.

Risk-based labeling based on the minimal eliciting doses (EDs) in sensitized populations is a potential replacement for precautionary allergen labeling of food allergens. We estimated the dose-response distribution for peanut allergen using data from double-blind placebo-controlled food challenges (DBPCFCs) conducted in the US at multiple sites, testing a population believed to be similar to the general U.S. food allergic population. Our final (placebo-adjusted) dataset included 548 challenges of 481 subjects. Bayesian hierarchical analysis facilitated model fitting, and accounted for variability associated with various levels of data organization. The data are best described using a complex hierarchical structure that accounts for inter-individual variability and variability across study locations or substudies. Bayesian model averaging could simultaneously consider the fit of multiple models, but the Weibull model dominated so strongly that model averaging was not needed. The ED01 and ED05 (and 95% credible intervals) are 0.052 (0.021, 0.13) and 0.49 (0.22, 0.97) mg peanut protein, respectively. Accounting for challenges with severe reactions at the LOAEL, by using the dose prior to the LOAEL as the new LOAEL, the ED01 drops to 0.029 (0.014, 0.074) mg peanut protein. Our results could aid in establishing improved food labeling guidelines in the management of food allergies.

Evaluation of concentration-response options for diacetyl in support of occupational risk assessment.

The current emphasis on occupational exposures to diacetyl has led to new research on its effects. We evaluated whether the data are sufficient to support a transition from a hazard-based risk management approach to a quantitative occupational risk assessment approach, characterized by developing a health-based occupational exposure limit (OEL). Inhalation health effects data were evaluated and issues and uncertainties related to occupational risk assessment needs were identified. A systematic hazard characterization, supported by both the toxicology and epidemiology literature, showed that the respiratory tract effects of diacetyl are the primary end points of relevance for developing an OEL. In an effort to provide a systematic approach for the analysis of the issues that need to be considered in developing an occupational risk assessment for diacetyl, a potential OEL was derived. A concentration-response assessment was completed using tracheobronchial effects in mice as the critical effect. The resulting benchmark concentration (lower bound estimate or BMCL) was adjusted to a human equivalent concentration of 1.8 ppm. A composite uncertainty factor of 10 was recommended to account for extrapolation from an adjusted BMCL from an animal study and for human variability in sensitivity and taking into account other uncertainties in the overall database. The resulting OEL recommendation of 0.2 ppm as a time-weighted average (TWA) was supported by the current occupational epidemiology literature. This evaluation showed that a health-based OEL value can be derived for diacetyl with moderate to high confidence.

A Bayesian network model for biomarker-based dose response.

A Bayesian network model was developed to integrate diverse types of data to conduct an exposure-dose-response assessment for benzene-induced acute myeloid leukemia (AML). The network approach was used to evaluate and compare individual biomarkers and quantitatively link the biomarkers along the exposure-disease continuum. The network was used to perform the biomarker-based dose-response analysis, and various other approaches to the dose-response analysis were conducted for comparison. The network-derived benchmark concentration was approximately an order of magnitude lower than that from the usual exposure concentration versus response approach, which suggests that the presence of more information in the low-dose region (where changes in biomarkers are detectable but effects on AML mortality are not) helps inform the description of the AML response at lower exposures. This work provides a quantitative approach for linking changes in biomarkers of effect both to exposure information and to changes in disease response. Such linkage can provide a scientifically valid point of departure that incorporates precursor dose-response information without being dependent on the difficult issue of a definition of adversity for precursors.

Evaluation of human relevance and mode of action for tunica vaginalis mesotheliomas resulting from oral exposure to acrylamide.

The human relevance and mode of action of acrylamide-related tunica vaginalis mesotheliomas (TVMs), a tumor of the scrotum, was evaluated based on the available data on acrylamide and general biology considerations. TVMs are found almost exclusively in F344 rats, suggesting an association with the hormonal milieu unique to F344s, and suggesting an association with Leydig cell tumors (LCTs), which occur in F344 rats at a very high incidence. These hypotheses are biologically plausible, but direct data on acrylamide were lacking for several key events; some of the gaps could be addressed based on other biology information. The data were not sufficient to identify a single definitive MOA. Multiple MOAs may apply, and some contribution from mutagenicity is plausible, along with a likely influence from LCTs or from the same hormonal changes that result in higher LCT incidence in F344 rats. Other MOAs, such as oxidative stress, may also apply. The data reviewed are not sufficient to distinguish between a causal relationship between LCTs and TVMs, and the hypothesis that these tumor types reflect a response to some shared influence (e.g., hormonal milieu of the F344 rat). Some of the plausible MOAs are not relevant to humans, while others are. In light of the very low incidence of TVMs in humans and the MOA data reviewed, the most appropriate upper bound estimate of the risk of acrylamide-related TVMs in humans is below de minimis levels.

Meta-regression analysis of the effects of dietary cholesterol intake on LDL and HDL cholesterol.

Background: Elevated low-density lipoprotein (LDL) cholesterol is a major risk factor for cardiovascular disease. Dietary guidance recommends reducing saturated fatty acid, trans fatty acid, and cholesterol intakes to reduce circulating LDL cholesterol. Cholesterol intake may also affect high-density lipoprotein (HDL)-cholesterol concentrations, but its impact has not been fully quantified. Objectives: The aims of this study were to investigate the dose-response relation between changes in dietary cholesterol intake and changes in lipoprotein-cholesterol markers for cardiovascular disease risk and to provide a reference for clinicians on how changes in dietary cholesterol intake affect circulating cholesterol concentrations, after accounting for intakes of fatty acids. Methods: We used a Bayesian approach to meta-regression analysis, which uses Markov chain Monte Carlo techniques, to assess the relation between the change in dietary cholesterol (adjusted for dietary fatty acids) and changes in LDL and HDL cholesterol based on the use of data from randomized dietary intervention trials. Results: Fifty-five studies (2652 subjects) were included in the analysis. The nonlinear Michaelis-Menten (MM) and Hill models best described the data across the full spectrum of dietary cholesterol changes studied (0-1500 mg/d). Mean predicted changes in LDL cholesterol for an increase of 100 mg dietary cholesterol/d were 1.90, 4.46, and 4.58 mg/dL for the linear, nonlinear MM, and Hill models, respectively. Conclusions: The change in dietary cholesterol was positively associated with the change in LDL-cholesterol concentration. The linear and MM models indicate that the change in dietary cholesterol is modestly inversely related to the change in circulating HDL-cholesterol concentrations in men but is positively related in women. The clinical implications of HDL-cholesterol changes associated with dietary cholesterol remain uncertain.