Practical guidance to evaluate in vitro dermal absorption studies for pesticide registration: An industry perspective.

While there are some regulatory assessment criteria available on how to generally evaluate dermal absorption (DA) studies for risk assessment purposes, practical guidance and examples are lacking. The current manuscript highlights the challenges in interpretating data from in vitro assays and proposes holistic data-based assessment strategies from an industry perspective. Inflexible decision criteria may be inadequate for real data and may lead to irrelevant DA estimates. We recommend the use of mean values for reasonably conservative DA estimates from in vitro studies. In cases where additional conservatism is needed, e.g., due to non-robust data and acute exposure scenarios, the upper 95% confidence interval of the mean may be appropriate. It is critical to review the data for potential outliers and we provide some example cases and strategies to identify aberrant responses. Some regional regulatory authorities require the evaluation of stratum corneum (SC) residue, but here, as a very simple pro-rata approach, we propose to review whether the predicted post 24-h absorption flux exceeds the predicted elimination flux by desquamation because otherwise it is not possible for the SC residue to contribute to systemic dose. Overall, the adjustment of DA estimates due to mass balance (normalization) is not recommended.

Expanding the toxicologist's statistical toolbox: Using effect size estimation and dose-response modelling for holistic assessments instead of generic testing.

It is tempting to base (eco-)toxicological assay evaluation solely on statistical significance tests. The approach is stringent, objective and facilitates binary decisions. However, tests according to null hypothesis statistical testing (NHST) are thought experiments that rely heavily on assumptions. The generic and unreflected application of statistical tests has been called "mindless" by Gigerenzer. While statistical tests have an appropriate application domain, the present work investigates how unreflected testing may affect toxicological assessments. Dunnett multiple-comparison and Williams trend testing and their compatibility intervals are compared with dose-response-modelling in case studies, where data do not follow textbook behavior, nor behave as expected from a toxicological point of view. In such cases, toxicological assessments based only on p-values may be biased and biological evaluations based on plausibility may be prioritized. If confidence in a negative assay outcome cannot be established, further data may be needed for a robust toxicological assessment.

Using historical control data in bioassays for regulatory toxicology.

Historical control data (HCD) consist of pooled control group responses from bioassays. These data must be collected and are often used or reported in regulatory toxicology studies for multiple purposes: as quality assurance for the test system, to help identify toxicological effects and their effect-size relevance and to address the statistical multiple comparison problem. The current manuscript reviews the various classical and potential new approaches for using HCD. Issues in current practice are identified and recommendations for improved use and discussion are provided. Furthermore, stakeholders are invited to discuss whether it is necessary to consider uncertainty when using HCD formally and statistically in toxicological discussions and whether binary inclusion/exclusion criteria for HCD should be revised to a tiered information contribution to assessments. Overall, the critical value of HCD in toxicological bioassays is highlighted when used in a weight-of-evidence assessment.

A Nonmechanistic Parametric Modeling Approach for Benchmark Dose Estimation of Event-Time Data.

We propose benchmark dose estimation for event-time data, using a two-step approach. This approach avoids estimation of complex models and has been previously shown to give robust results for summarizing relevant parameters for risk assessment. In the first step, the probability of the event of interest to occur (in a certain time interval) is described as a function of time, resulting in an event-time model; such a model is fitted allowing an individual curve for each dose, and relevant estimates are extracted. In the second step, a dose-response model is fitted to the estimates of t50 obtained from the event-time model in the first step. Given a predefined benchmark response, the benchmark dose is then estimated from the resulting model. This novel approach is demonstrated in two examples. Our application of the time-to-event model showed a gain in power compared to the traditional analysis of end-of-study summary data.

A Review of Recent Advances in Benchmark Dose Methodology.

In this review, recent methodological developments for the benchmark dose (BMD) methodology are summarized. Specifically, we introduce the advances for the main steps in BMD derivation: selecting the procedure for defining a BMD from a predefined benchmark response (BMR), setting a BMR, selecting a dose-response model, and estimating the corresponding BMD lower limit (BMDL). Although the last decade has shown major progress in the development of BMD methodology, there is still room for improvement. Remaining challenges are the implementation of new statistical methods in user-friendly software and the lack of consensus about how to derive the BMDL.

A comparison of approaches for simultaneous inference of fixed effects for multiple outcomes using linear mixed models.

Longitudinal studies with multiple outcomes often pose challenges for the statistical analysis. A joint model including all outcomes has the advantage of incorporating the simultaneous behavior but is often difficult to fit due to computational challenges. We consider 2 alternative approaches to quantify and assess the loss in efficiency as compared with joint modelling when evaluating fixed effects. The first approach is pairwise fitting of pseudolikelihood functions for pairs of outcomes. The second approach recovers correlations between parameter estimates across multiple marginal linear mixed models. The methods are evaluated in terms of a data example both from a study on the effects of milk protein on health in young adolescents and in an extensive simulation study. We find that the 2 alternatives give similar results in settings where an exchangeability condition is met, but otherwise, pairwise fitting shows a larger loss in efficiency than the marginal models approach. Using an alternative to the joint modelling strategy will lead to some but not necessarily a large loss of efficiency for small sample sizes.

Factors associated with non-participation in cervical cancer screening - A nationwide study of nearly half a million women in Denmark.

Cervical cancer occurs most often in under-screened women. In this nationwide register study, we described differences in sociodemographic characteristics between passive and active non-participants and examined socio-demographic characteristics, reproductive history, and mental and physical health as potential determinants for passive non-participation compared with participation in the Danish cervical cancer screening program. Screening history in women aged 23-49 years invited for cervical cancer screening in 2008-2009 was retrieved from the Danish Pathology Databank with information about dates of invitation and unsubscription. We identified participants (n = 402,984), active non-participants (n = 10,251) and passive non-participants (n = 63,435) within four years following baseline invitation and retrieved data about the study population from high-quality registries. We examined differences in socio-demographic characteristics of passive and active non-participants, and used multiple logistic regression analyses to identify potential determinants of passive non-participation. We found that active and passive non-participants differed in relation to socio-demography. When compared with screening participants, the odds of passive non-participation was increased in women who originated from less developed countries; were unmarried; had basic education or low income; had four or more children; smoked during pregnancy; had multiple induced abortions; or had a history of obesity, intoxicant abuse or schizophrenia or other psychoses. In conclusion, in this nationwide, prospective, population-based study, differences in socio-demographic characteristics between passive and active non-participants were found. Furthermore, sociodemography, reproductive history, and mental and physical health were determinants for passive non-participation. Addressing inequalities in screening attendance may help to further decrease the incidence of and mortality from cervical cancer.

Cigarette smoking is associated with adverse survival among women with ovarian cancer: Results from a pooled analysis of 19 studies.

Cigarette smoking is associated with an increased risk of developing mucinous ovarian tumors but whether it is associated with ovarian cancer survival overall or for the different histotypes is unestablished. Furthermore, it is unknown whether the association between cigarette smoking and survival differs according to strata of ovarian cancer stage at diagnosis. In a large pooled analysis, we evaluated the association between various measures of cigarette smoking and survival among women with epithelial ovarian cancer. We obtained data from 19 case-control studies in the Ovarian Cancer Association Consortium (OCAC), including 9,114 women diagnosed with ovarian cancer. Cox regression models were used to estimate adjusted study-specific hazard ratios (HRs), which were combined into pooled hazard ratios (pHR) with corresponding 95% confidence intervals (CIs) under random effects models. Overall, 5,149 (57%) women died during a median follow-up period of 7.0 years. Among women diagnosed with ovarian cancer, both current (pHR = 1.17, 95% CI: 1.08-1.28) and former smokers (pHR = 1.10, 95% CI: 1.02-1.18) had worse survival compared with never smoking women. In histotype-stratified analyses, associations were observed for mucinous (current smoking: pHR = 1.91, 95% CI: 1.01-3.65) and serous histotypes (current smoking: pHR = 1.11, 95% CI: 1.00-1.23; former smoking: pHR = 1.12, 95% CI: 1.04-1.20). Further, our results suggested that current smoking has a greater impact on survival among women with localized than disseminated disease. The identification of cigarette smoking as a modifiable factor associated with survival has potential clinical importance as a focus area to improve ovarian cancer prognosis.

Acute disruption of glucagon secretion or action does not improve glucose tolerance in an insulin-deficient mouse model of diabetes.

AIMS/HYPOTHESIS: Normal glucose metabolism depends on pancreatic secretion of insulin and glucagon. The bihormonal hypothesis states that while lack of insulin leads to glucose underutilisation, glucagon excess is the principal factor in diabetic glucose overproduction. A recent study reported that streptozotocin-treated glucagon receptor knockout mice have normal glucose tolerance. We investigated the impact of acute disruption of glucagon secretin or action in a mouse model of severe diabetes by three different approaches: (1) alpha cell elimination; (2) glucagon immunoneutralisation; and (3) glucagon receptor antagonism, in order to evaluate the effect of these on glucose tolerance. METHODS: Severe diabetes was induced in transgenic and wild-type mice by streptozotocin. Glucose metabolism was investigated using OGTT in transgenic mice with the human diphtheria toxin receptor expressed in proglucagon producing cells allowing for diphtheria toxin (DT)-induced alpha cell ablation and in mice treated with either a specific high affinity glucagon antibody or a specific glucagon receptor antagonist. RESULTS: Near-total alpha cell elimination was induced in transgenic mice upon DT administration and resulted in a massive decrease in pancreatic glucagon content. Oral glucose tolerance in diabetic mice was neither affected by glucagon immunoneutralisation, glucagon receptor antagonism, nor alpha cell removal, but did not deteriorate further compared with mice with intact alpha cell mass. CONCLUSIONS/INTERPRETATION: Disruption of glucagon action/secretion did not improve glucose tolerance in diabetic mice. Near-total alpha cell elimination may have prevented further deterioration. Our findings support insulin lack as the major factor underlying hyperglycaemia in beta cell-deficient diabetes.

Factors associated with type-specific persistence of high-risk human papillomavirus infection: A population-based study.

Persistent genital infection with high-risk (HR) human papillomavirus (HPV) is a prerequisite for cervical cancer development. The aim of this study was to identify factors associated with type-specific persistence of HR HPV infections. From a population-based cohort of 40,399 women participating in cervical cancer screening established during 2002-2005, we selected all HR HPV-positive women (N = 7,778). During follow-up (2005-2008), we collected cervical samples from these women and tested them for HPV DNA to determine type-specific HR HPV persistence in the interval 1-4.5 years after enrolment. Data on hospitalisations, prescriptions and socioeconomic factors were obtained from nationwide registers. Women with abnormal cytology at baseline or who had undergone conisation during follow-up were excluded. Factors associated with persistence were identified by logistic regression analysis. The overall rate of HR HPV persistence was 31.4%. The risk for persistence was significantly increased among women with a previous episode of genital warts (OR, 1.35; 95% CI, 1.04-1.74), current use of oral contraceptives (OR, 1.35; 95% CI, 1.13-1.63) or use of systemic glucocorticoids (OR, 2.04; 95% CI, 1.16-3.56). The number of pregnancies or births or use of a hormonal intrauterine device, hormonal therapy or nonsteroidal anti-inflammatory drugs was not associated with risk for HR HPV persistence. A history of genital warts and current use of oral contraceptives or systemic glucocorticoids increased the risk, potentially indicating a decreased immune response to HPV infection. These findings suggest that host immune response characteristics are important in HR HPV persistence and consequently in cervical cancer development.

Evaluation of multi-outcome longitudinal studies.

Evaluation of intervention effects on multiple outcomes is a common scenario in clinical studies. In longitudinal studies, such evaluation is a challenge if one wishes to adequately capture simultaneous data behavior. In this situation, a common approach is to analyze each outcome separately. As a result, multiple statistical statements describing the intervention effect need to be reported and an adjustment for multiple testing is necessary. This is typically done by means of the Bonferroni procedure, which does not take into account the correlation between outcomes, thus resulting in overly conservative conclusions. We propose an alternative approach for multiplicity adjustment that incorporates dependence between outcomes, resulting in an appreciably less conservative evaluation. The ability of the proposed method to control the familywise error rate is evaluated in a simulation study, and the applicability of the method is demonstrated in two examples from the literature.

The impact of early growth patterns and infant feeding on body composition at 3 years of age.

Early excessive weight gain is positively associated with later obesity, and yet the effect of weight gain during specific periods and the impact of infant feeding practices are debated. The objective of the present study was to examine the impact of weight gain in periods of early childhood on body composition at 3 years, and whether infant feeding modified the relationship between early growth and body composition at 3 years. We studied 233 children from the prospective cohort study, SKOT (in Danish: Småbørns Kost og Trivsel). Birth weight z-scores (BWZ) and change in weight-for-age z-scores (WAZ) from 0 to 5, 5 to 9, 9 to 18 and 18 to 36 months were analysed for relations with body composition (anthropometry and bioelectrical impedance) at 3 years by multivariate regression analysis. BWZ and change in WAZ from 0 to 5 months were positively associated with BMI, fat mass index (FMI) and fat-free mass index (FFMI) at 3 years. Full breastfeeding for 6 months (compared to less than 1 month) eliminated the effect of early growth (P = 0.01). Full breastfeeding for 6 months (compared to less than 1 month) also eliminated the positive relation between BWZ and FMI (P = 0.009). No effect modification of infant feeding was found for FFMI. In conclusion, high birth weight and rapid growth from 0 to 5 months were associated with increased FMI and FFMI at 3 years. Longer duration of full breastfeeding reduced the effect of birth weight and early weight gain on fat mass.

Simultaneous inference for model averaging of derived parameters.

Model averaging is a useful approach for capturing uncertainty due to model selection. Currently, this uncertainty is often quantified by means of approximations that do not easily extend to simultaneous inference. Moreover, in practice there is a need for both model averaging and simultaneous inference for derived parameters calculated in an after-fitting step. We propose a method for obtaining asymptotically correct standard errors for one or several model-averaged estimates of derived parameters and for obtaining simultaneous confidence intervals that asymptotically control the family-wise Type I error rate. The performance of the method in terms of coverage is evaluated using a simulation study and the applicability of the method is demonstrated by means of three concrete examples.

Validity of anthropometric measurements to assess body composition, including muscle mass, in 3-year-old children from the SKOT cohort.

Nutritional status of children is commonly assessed by anthropometry both in under and overnutrition. The link between anthropometry and body fat, the body compartment most affected by overnutrition, is well known, but the link with muscle mass, the body compartment most depleted in undernutrition, associated with infections, remains unknown. In this study, we examined the relationship between common anthropometric indices and body composition measured by dual-energy X-ray absorptiometry (DEXA) in a sample of 121 healthy 3-year-old Danish children. Appendicular (arms and legs) lean mass was used to estimate muscle mass. Overall, anthropometric measures were more effective to measure absolute size of fat, lean and muscle mass than their relative sizes. Proportion of the variance explained by anthropometry was 79% for lean mass, 76% for fat mass and 74% for muscle mass. For fat mass and lean mass expressed as percentage of total body mass, this proportion was 51% and 66%, respectively; and for muscle mass as percentage of lean mass it was 34%. All the best reduced multivariate models included weight, skinfold and gender except the model estimating the proportion of muscle mass in lean body mass, which included only mid-upper arm circumference and subscapular skinfold. The power of height in the weight-to-height ratio to determine fat mass proportion was 1.71 with a 95% confidence interval (0.83-2.60) including the value of 2 used in body mass index (BMI). Limitations of anthropometry to assess body composition, and especially for muscle mass as a proportion of lean mass, should be acknowledged.

Modelling pesticide degradation and leaching in conservation agriculture: Effect of no-till and mulching.

No-till and mulching are typical management operations in conservation agriculture (CA). To model pesticide degradation and leaching under a CA scenario, as compared to a conventional-tillage scenario (CT), the mulch module of the agro-hydrological model Daisy was extended. A Daisy soil column was parameterized with measurements of topsoil, mulch, and a realistic subsoil, and tested against published experimental data of pesticide fate in laboratory soil columns covered by mulch. Uncertainty and sensitivity analyses of the new Daisy version were conducted for a series of weather, soil, pesticide, and mulch parameters, using 4939 Monte Carlo simulations under each scenario. Results showed that there was no systematic difference in pesticide leaching from the topsoil (to the subsoil and directly to drains via drain-connected biopores) between CA and CT, but pesticide degradation and sorption were significantly different; degradation in the mulch and uppermost soil surface layer (0-3.5 cm) was larger in CA while degradation was larger in CT when considering the whole topsoil (0-30 cm). This difference for the whole topsoil could be explained by pesticide interception in CA in the part of the mulch not in direct contact with the soil where degradation is assumed not to occur. The sensitivity analysis highlighted non-influential parameters and seven parameters out of twenty-five to be better estimated to improve the accuracy of the predictions.

Use of antibiotics during pregnancy increases the risk of asthma in early childhood.

OBJECTIVES: To investigate the hypothesis that mother's use of antibiotics in pregnancy could influence asthma and eczema in early life. STUDY DESIGN: Subjects were included from the Copenhagen Prospective Study on Asthma in Childhood cohort of children born of mothers with asthma (N = 411). Severe asthma exacerbations and eczema were diagnosed by research unit physicians. Replication was sought in children from the Danish National Birth Cohort (N = 30 675). Asthma outcomes were hospitalization and use of inhaled corticosteroids. Eczema was defined by an algorithm developed from cases of clinically verified eczema. All children were followed to age 5 years in a cohort study design. RESULTS: The Copenhagen Prospective Study on Asthma in Childhood data showed increased risk of asthma exacerbation (hazard ratio 1.98 [95% CI 1.08-3.63]) if mothers had used antibiotics during third trimester. The Danish National Birth Cohort confirmed increased risk of asthma hospitalization (hazard ratio 1.17 [1.00-1.36]), and inhaled corticosteroids (1.18 [1.10-1.27]) in the children if mothers used antibiotics any time during pregnancy. In the subgroup of mothers using antibiotics for nonrespiratory infection, the children also had increased risk of asthma. CONCLUSION: We found increased risk of asthma associated with maternal antibiotic use in a clinical study of a birth cohort with increased risk of asthma and replicated this finding in an unselected national birth cohort, and in a subgroup using antibiotics for nonrespiratory infections. This supports a role for bacterial ecology in pre- or perinatal life for the development of asthma.

Higher protein diets consumed ad libitum improve cardiovascular risk markers in children of overweight parents from eight European countries.

Dietary strategies to improve early cardiovascular markers in overweight children are needed. We investigated the effect of dietary protein and glycemic index (GI) on cardiovascular markers and metabolic syndrome (MetS) scores in 5- to 18-y-old children of overweight/obese parents from 8 European centers. Families were randomized to 1 of 5 diets consumed ad libitum: high protein (HP) or low protein (LP) combined with high GI (HGI) or low GI (LGI), or a control diet. At 6 centers, families received dietary instruction (instruction centers); at 2 centers, free foods were also provided (supermarket centers). Diet, anthropometry, blood pressure, and serum cardiovascular markers (lipid profile, glucose regulation, and inflammation) were measured in 253 children at baseline, 1 mo, and/or 6 mo. Protein intake was higher in the HP groups (19.9 ± 1.3% energy) than in the LP groups at 6 mo (16.8 ± 1.2% energy) (P = 0.001). The GI was 4.0 points lower (95% CI: 2.1, 6.1) in the LGI compared with the HGI groups (P < 0.001). In the supermarket centers, the HP and LP groups differed more in protein intake than did the groups in the instruction centers (P = 0.009), indicating better compliance. The HP diets evoked a 2.7-cm (95% CI: 0.9, 5.1) smaller waist circumference and a 0.25-mmol/L (95% CI: 0.09, 0.41) lower serum LDL cholesterol compared with the LP diets at 6 mo (P < 0.007). In a separate supermarket center analysis, the HP compared with LP diets reduced waist circumference (P = 0.004), blood pressure (P < 0.01), serum insulin (P = 0.013), and homeostasis model of assessment-insulin resistance (P = 0.016). In the instruction centers, the HP compared with the LP diets reduced LDL cholesterol (P = 0.004). No consistent effect of GI was seen and the MetS scores were not affected. In conclusion, increased protein intake improved cardiovascular markers in high-risk children, particularly in those undergoing most intensive intervention.

Interaction between asthma and lung function growth in early life.

RATIONALE: The causal direction between asthma and lung function deficit is unknown, but important for the focus of preventive measures and research into the origins of asthma. OBJECTIVES: To analyze the interaction between lung function development and asthma from birth to 7 years of age. METHODS: The Copenhagen Prospective Studies on Asthma in Childhood is a prospective clinical study of a birth cohort of 411 at-risk children. Spirometry was completed in 403 (98%) neonates and again by age 7 in 317 children (77%). MEASUREMENTS AND MAIN RESULTS: Neonatal spirometry and bronchial responsiveness to methacholine was measured during sedation by forced flow-volume measurements. Asthma was diagnosed prospectively from daily diary cards and clinic visits every 6 months. Children with asthma by age 7 (14%) already had a significant airflow deficit as neonates (forced expiratory flow at 50% of vital capacity second in neonates reduced by 0.34 z score by 1 mo; P = 0.03). This deficit progressed significantly during early childhood (forced expiratory flow at 0.5 seconds in neonates at age 7 reduced by 0.82 z score by age 7; P < 0.0001), suggesting that approximately 40% of the airflow deficit associated with asthma is present at birth, whereas 60% develops with clinical disease. Environmental tobacco exposure, but not allergic sensitization, also hampered airflow growth. Bronchial responsiveness to methacholine in the neonates was associated with the development of asthma (P = 0.01). CONCLUSIONS: Children developing asthma by age 7 had a lung function deficit and increased bronchial responsiveness as neonates. This lung function deficit progressed to age 7. Therefore, research into the origins and prevention of asthma should consider early life before and after birth.

Using R in Regulatory Toxicology.

Statistical analyses are an essential part of regulatory toxicological evaluations. While projects would be ideally monitored by both toxicologists and statisticians, this is often not possible in practice. Hence, toxicologists should be trained in some common statistical approaches but also need a tool for statistical evaluations. Due to transparency needed in regulatory processes and standard tests that can be evaluated with template approaches, the freely available open-source statistical software R may be suitable. R is a well-established software in the statistical community. The principal input method is via software code, which is both benefit and weakness of the tool. It is increasingly used by regulating authorities globally and can be easily extended by software packages, e.g., for new statistical functions and features. This manuscript outlines how R can be used in regulatory toxicology, allowing toxicologists to perform all regulatory required data evaluations in a single software solution. Practical applications are shown in case studies on simulated and experimental data. The examples cover a) Dunnett testing of treatment groups against a common control and in relation to a biological relevance threshold, assessing the test's assumptions and plotting the results; b) dose-response analysis and benchmark dose derivation for chronic kidney inflammation as a function of Pyridine; and c) graphical/exploratory data analysis of previously published developmental neurotoxicity data for Chlorpyrifos.

Benchmark dose modelling in regulatory ecotoxicology, a potential tool in pest management.

For several authorities, benchmark dose (BMD) methodology has become the recommended approach by which to derive reference values for risk assessment. However, in practice, the BMD approach is not standard use in risk assessment for pesticides where the no observed adverse effect level, lowest observed adverse effect level and effective dose (ED50 or EDx ) prevail. Regression-based BMD and the benchmark dose lower confidence limit (BMDL) have several advantages, such as utilizing more information from the generated data and being less dependent on tested dose levels. However, the BMD approach requires some degree of expert knowledge for defining an appropriate risk level for estimating the BMD and using more sophisticated statistical methods to calculate BMD and BMDL. The BMD approach is one way to move away from p value-based binary decision-making towards putting the weight on effect sizes. We review the advantages and disadvantages of focusing on the BMD approach for risk assessment of pesticides. Further, we discuss potential applications in efficacy trials for pest management purposes. © 2021 Society of Chemical Industry.