COVID-19 vaccine uptake among children and adolescents in Norway: A comprehensive registry-based cohort study of over 800,000 individuals

Factors related with COVID-19 vaccine uptake in children and adolescents in Norway remain unclear, despite this being useful knowledge for future pandemic preparedness. This study aimed to comprehensively examine individual and familial factors associated with vaccine uptake in children and adolescents in Norway. We utilized nationwide registry-data from various health registries and Statistics Norway, encompassing all children and adolescents living in Norway during the pandemic, until 31-Dec-2022. Vaccine uptake is defined as receiving at least one dose of COVID-19 vaccine. We employed a forward stepwise logistic regression model and a random forest machine-learning algorithm to explore the relationship between vaccine uptake and socio-cultural, demographic, and health-related factors. We included 423,548 5-11-year-olds, 269,830 12-15-year-olds, and 120,854 16-17-year-olds. Vaccine uptake in these three groups was respectively 2.6 %, 73.3 %, and 87.3 %. Factors associated with vaccine uptake varied by age group. In youngest children, immigrant background (Odds-ratio (OR) = 1.58, 95 % confidence interval (CI) (1.14-2.19)), born extremely preterm (OR = 2.38, 95 % CI (1.60-3.54)), having risk of severe COVID-19 (OR = 5.40, 95 % CI (4.69-6.23) and maternal COVID-19 vaccination (OR = 6.34, 95 % CI (5.35-7.53)) were positively associated with vaccine uptake. The latter two factors were also strongly, positively associated with vaccine uptake in 12-15-year-olds, while previous SARS-CoV-2 infection was negatively associated (OR = 0.12, 95 % CI (0.11-0.14). Similar findings were observed in 16-17-year-olds. COVID-19 vaccine uptake differed markedly by age group, and major associated factors included socio-demographics and parental COVID-19 vaccination status, prior SARS-CoV-2 infection, but also being born premature and having moderate or high risk of severe COVID-19.

Urinary phthalate/DINCH metabolites associations with kisspeptin and reproductive hormones in teenagers: A cross-sectional study from the HBM4EU aligned studies

BACKGROUND

Exposure to phthalate/DINCH metabolites can induce human reproductive toxicity, however, their endocrine-disrupting mechanisms are not fully elucidated.

OBJECTIVE

To investigate the association between concentrations of phthalate/DINCH metabolites, serum kisspeptin, and reproductive hormones among European teenagers from three of the HBM4EU Aligned Studies.

METHODS

In 733 Belgian (FLEHS IV study), Slovak (PCB cohort follow-up), and Spanish (BEA study) teenagers, ten phthalate and two DINCH metabolites were measured in urine by high-performance liquid chromatography-tandem mass spectrometry. Serum kisspeptin (kiss54) protein, follicle-stimulating hormone (FSH), total testosterone (TT), estradiol (E2), and sex hormone-binding globulin (SHBG) levels were measured by immunosorbent assays. Free Androgen Index (FAI) was calculated as a proxy of free testosterone. Adjusted sex-stratified linear regression models for individual studies, mixed effect models (LME) accounting for random effects for pooled studies, and g-computation and Bayesian kernel machine regression (BKMR) models for the phthalate/DINCH mixture were performed.

RESULTS

The LME suggested that each IQR increase in ln-transformed levels of several phthalates was associated with lower kisspeptin [MnBP: %change (95%CI): -2.8 (-4.2;-0.4); MEHP: -1.4 (-3.4,0.2)] and higher FSH [∑DINP: 11.8 (-0.6;25.1)] levels in females from pooled studies. G-computation showed that the phthalates/DINCH mixture was associated with lower kisspeptin [-4.28 (-8.07;-0.34)] and higher FSH [22.13 (0.5;48.4)] also in females; BKMR showed similar although non-significant pattern. In males, higher phthalates metabolites [MEHP: -12.22 (-21.09;-1.18); oxo-MEHP: -12.73 (-22.34;-1.93)] were associated with lower TT and FAI, although higher DINCH [OH-MINCH: 16.31 (6.23;27.35), cx-MINCH: 16.80 (7.03;27.46), ∑DINCH: 17.37 (7.26;29.74)] were associated with higher TT levels. No mixture associations were found in males.

CONCLUSION

We observed sex-specific associations between urinary concentrations of phthalate/DINCH metabolites and the panel of selected effect biomarkers (kisspeptin and reproductive hormones). This suggests that exposure to phthalates would be associated with changes in kisspeptin levels, which would affect the HPG axis and thus influence reproductive health. However, further research is needed, particularly for phthalate replacements such as DINCH.

PFAS association with kisspeptin and sex hormones in teenagers of the HBM4EU aligned studies

Exposure to Perfluoroalkyl acids (PFAS) can impair human reproductive function, e.g., by delaying or advancing puberty, although their mechanisms of action are not fully understood. We therefore set out to evaluate the relationship between serum PFAS levels, both individually and as a mixture, on the Hypothalamic-Pituitary-Gonadal (HPG) axis by analyzing serum levels of reproductive hormones and also kisspeptin in European teenagers participating in three of the HBM4EU Aligned Studies. For this purpose, PFAS compounds were measured in 733 teenagers from Belgium (FLEHS IV study), Slovakia (PCB cohort follow-up), and Spain (BEA study) by high performance liquid chromatography-tandem mass spectrometry (HPLC/MS) in laboratories under the HBM4EU quality assurance quality control (QA/QC) program. In the same serum samples, kisspeptin 54 (kiss-54) protein, follicle-stimulating hormone (FSH), total testosterone (TT), estradiol (E2), and sex hormone-binding globulin (SHBG) levels were also measured using immunosorbent assays. Sex-stratified single pollutant linear regression models for separate studies, mixed single pollutant models accounting for random effects for pooled studies, and g-computation and Bayesian kernel machine regression (BKMR) models for the mixture of the three most available (PFNA, PFOA, and PFOS) were fit. PFAS associations with reproductive markers differed according to sex. Each natural log-unit increase of PFOA, PFNA, and PFOS were associated with higher TT [18.41 (6.18; 32.31), 15.60 (7.25; 24.61), 14.68 (6.18; 24.61), respectively] in girls, in the pooled analysis (all studies together). In males, G-computation showed that PFAS mixture was associated with lower FSH levels [-10.51 (-18.81;-1.36)]. The BKMR showed the same patterns observed in G-computation, including a significant increase on male Kiss-54 and SHBG levels. Overall, effect biomarkers may enhance the current epidemiological knowledge regarding the adverse effect of PFAS in human HPG axis, although further research is warranted.

PFAS and Phthalate/DINCH Exposure in Association with Age at Menarche in Teenagers of the HBM4EU Aligned Studies

Early puberty has been found to be associated with adverse health outcomes such as metabolic and cardiovascular diseases and hormone-dependent cancers. The decrease in age at menarche observed during the past decades has been linked to an increased exposure to endocrine-disrupting compounds (EDCs). Evidence for the association between PFAS and phthalate exposure and menarche onset, however, is inconsistent. We studied the association between PFAS and phthalate/DINCH exposure and age at menarche using data of 514 teenagers (12 to 18 years) from four aligned studies of the Human Biomonitoring for Europe initiative (HBM4EU): Riksmaten Adolescents 2016-2017 (Sweden), PCB cohort (follow-up; Slovakia), GerES V-sub (Germany), and FLEHS IV (Belgium). PFAS concentrations were measured in blood, and phthalate/DINCH concentrations in urine. We assessed the role of each individual pollutant within the context of the others, by using different multi-pollutant approaches, adjusting for age, age- and sex-standardized body mass index z-score and household educational level. Exposure to di(2-ethylhexyl) phthalate (DEHP), especially mono(2-ethyl-5-hydroxyhexyl) phthalate (5OH-MEHP), was associated with an earlier age at menarche, with estimates per interquartile fold change in 5OH-MEHP ranging from -0.34 to -0.12 years in the different models. Findings from this study indicated associations between age at menarche and some specific EDCs at concentrations detected in the general European population, but due to the study design (menarche onset preceded the chemical measurements), caution is needed in the interpretation of causality.

Cross-sectional associations between exposure to per- and polyfluoroalkyl substances and body mass index among European teenagers in the HBM4EU aligned studies

Per- and polyfluoroalkyl substances (PFAS) are widespread pollutants that may impact youth adiposity patterns. We investigated cross-sectional associations between PFAS and body mass index (BMI) in teenagers/adolescents across nine European countries within the Human Biomonitoring for Europe (HBM4EU) initiative. We used data from 1957 teenagers (12-18 yrs) that were part of the HBM4EU aligned studies, consisting of nine HBM studies (NEBII, Norway; Riksmaten Adolescents 2016-17, Sweden; PCB cohort (follow-up), Slovakia; SLO CRP, Slovenia; CROME, Greece; BEA, Spain; ESTEBAN, France; FLEHS IV, Belgium; GerES V-sub, Germany). Twelve PFAS were measured in blood, whilst weight and height were measured by field nurse/physician or self-reported in questionnaires. We assessed associations between PFAS and age- and sex-adjusted BMI z-scores using linear and logistic regression adjusted for potential confounders. Random-effects meta-analysis and mixed effects models were used to pool studies. We assessed mixture effects using molar sums of exposure biomarkers with toxicological/structural similarities and quantile g-computation. In all studies, the highest concentrations of PFAS were PFOS (medians ranging from 1.34 to 2.79 μg/L). There was a tendency for negative associations with BMI z-scores for all PFAS (except for PFHxS and PFHpS), which was borderline significant for the molar sum of [PFOA and PFNA] and significant for single PFOA [β-coefficient (95% CI) per interquartile range fold change = -0.06 (-0.17, 0.00) and -0.08 (-0.15, -0.01), respectively]. Mixture assessment indicated similar negative associations of the total mixture of [PFOA, PFNA, PFHxS and PFOS] with BMI z-score, but not all compounds showed associations in the same direction: whilst [PFOA, PFNA and PFOS] were negatively associated, [PFHxS] associated positively with BMI z-score. Our results indicated a tendency for associations of relatively low PFAS concentrations with lower BMI in European teenagers. More prospective research is needed to investigate this potential relationship and its implications for health later in life.

Role of Physiologically Based Kinetic modelling in addressing environmental chemical mixtures - A review

•

The availability and applicability of Physiologically Based Kinetic (PBK) models for mixtures is reviewed.

•

PBK models can support risk assessment of mixtures by incorporating the toxicokinetic processes.

•

Quantitative structure-activity relationship (QSAR) models can be used to fill data gaps in PBK modelling.

•

PBK models for mixtures can be improved by including various types of interactions.

The role of Physiologically Based Kinetic (PBK) modelling in assessing mixture toxicology has been growing for the last three decades. It has been widely used to investigate and address interactions in mixtures. This review describes the current state-of-the-art of PBK models for chemical mixtures and to evaluate the applications of PBK modelling for mixtures with emphasis on their role in chemical risk assessment. A total of 35 mixture PBK models were included after searching web resources (Scopus, PubMed, Web of Science, and Google Scholar), screening for duplicates, and excluding articles based on eligibility criteria. Binary mixtures and volatile organic compounds accounted for two-thirds of the chemical mixtures identified. The most common exposure route and modelled system were found to be inhalation and rats respectively. Twenty two (22) models were for binary mixtures, 5 for ternary mixtures, 3 for quaternary mixtures, and 5 for complex mixtures. Both bottom-up and top-down PBK modelling approaches are described. Whereas bottom-up approaches are based on a series of binary interactions, top-down approaches are based on the lumping of mixture components. Competitive inhibition is the most common type of interaction among the various types of mixtures, and usually becomes a concern at concentrations higher than environmental exposure levels. It leads to reduced biotransformation that either means a decrease in the amount of toxic metabolite formation or an increase in toxic parent chemical accumulation. The consequence is either lower or higher toxicity compared to that estimated for the mixture based on the additivity principle. Therefore, PBK modelling can play a central role in predicting interactions in chemical mixture risk assessment.

Next generation physiologically based kinetic (NG-PBK) models in support of regulatory decision making

The fields of toxicology and chemical risk assessment seek to reduce, and eventually replace, the use of animals for the prediction of toxicity in humans. In this context, physiologically based kinetic (PBK) modelling based on in vitro and in silico kinetic data has the potential to a play significant role in reducing animal testing, by providing a methodology capable of incorporating in vitro human data to facilitate the development of in vitro to in vivo extrapolation of hazard information. In the present article, we discuss the challenges in: 1) applying PBK modelling to support regulatory decision making under the toxicology and risk-assessment paradigm shift towards animal replacement; 2) constructing PBK models without in vivo animal kinetic data, while relying solely on in vitro or in silico methods for model parameterization; and 3) assessing the validity and credibility of PBK models built largely using non-animal data. The strengths, uncertainties, and limitations of PBK models developed using in vitro or in silico data are discussed in an effort to establish a higher degree of confidence in the application of such models in a regulatory context. The article summarises the outcome of an expert workshop hosted by the European Commission Joint Research Centre (EC-JRC) - European Union Reference Laboratory for Alternatives to Animal Testing (EURL ECVAM), on "Physiologically-Based Kinetic modelling in risk assessment - reaching a whole new level in regulatory decision-making" held in Ispra, Italy, in November 2016, along with results from an international survey conducted in 2017 and recently reported activities occurring within the PBK modelling field. The discussions presented herein highlight the potential applications of next generation (NG)-PBK modelling, based on new data streams.