Fatty acids and recurrence of major depressive disorder: combined analysis of two Dutch clinical cohorts.

OBJECTIVE: Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acid (PUFA) alterations in patients with major depressive disorder (MDD) have been shown to persist after remission. Whether these alterations are risk factors for MDD recurrence remains unknown. Here, we examined whether fatty acids predict time until MDD recurrence in remitted MDD patients. METHODS: Data were used from remitted MDD patients of the Netherlands Study of Depression and Anxiety (n = 356) and the Depression Evaluation Longitudinal Therapy Assessment studies (n = 118). Associations of FAs with time until MDD recurrence up to 8-year follow-up were analyzed using Cox regression analyses. Study-specific estimates were pooled using mega- and meta-analysis techniques. RESULTS: 27.5% (NESDA) and 56.8% (DELTA) participants had an MDD recurrence. Pooled results showed that no FA was significantly associated with time until MDD recurrence (n-3 PUFAs: hazard ratio (HR) = 1.17, 95% confidence interval (CI) = 0.98-1.41, P = 0.082; n-6 PUFAs: HR = 1.08, 95% CI = 0.84-1.38, P = 0.55). CONCLUSION: In remitted MDD patients, circulating PUFAs were not associated with prospective risk of MDD recurrence. Consequently, circulating PUFAs are unlikely to reflect a vulnerability marker for recurrence, so correcting n-3 PUFA 'deficits' through supplementation does not seem a promising option to prevent MDD recurrence.

Unraveling the genetic architecture of major depressive disorder: merits and pitfalls of the approaches used in genome-wide association studies.

To identify genetic risk loci for major depressive disorder (MDD), two broad study design approaches have been applied: (1) to maximize sample size by combining data from different phenotype assessment modalities (e.g. clinical interview, self-report questionnaires) and (2) to reduce phenotypic heterogeneity through selecting more homogenous MDD subtypes. The value of these strategies has been debated. In this review, we summarize the most recent findings of large genomic studies that applied these approaches, and we highlight the merits and pitfalls of both approaches with particular attention to methodological and psychometric issues. We also discuss the results of analyses that investigated the heterogeneity of MDD. We conclude that both study designs are essential for further research. So far, increasing sample size has led to the identification of a relatively high number of genomic loci linked to depression. However, part of the identified variants may be related to a phenotype common to internalizing disorders and related traits. As such, samples containing detailed clinical information are needed to dissect depression heterogeneity and enable the potential identification of variants specific to a more restricted MDD phenotype. A balanced portfolio reconciling both study design approaches is the optimal approach to progress further in unraveling the genetic architecture of depression.

Metabolic and inflammatory markers: associations with individual depressive symptoms.

BACKGROUND: Literature has shown that obesity, metabolic syndrome and inflammation are associated with depression, however, evidence suggests that these associations are specific to atypical depression. Which of the atypical symptoms are driving associations with obesity-related outcomes and inflammation is unknown. We evaluated associations between individual symptoms of depression (both atypical and non-atypical) and body mass index (BMI), metabolic syndrome components and inflammatory markers. METHODS: We included 808 persons with a current diagnosis of depression participating in the Netherlands Study of Depression and Anxiety (67% female, mean age 41.6 years). Depressive symptoms were derived from the Composite International Diagnostic Interview and the Inventory of Depressive Symptomatology. Univariable and multivariable regression analyses adjusting for sex, age, educational level, depression severity, current smoking, physical activity, anti-inflammatory medication use, and statin use were performed. RESULTS: Increased appetite was positively associated with BMI, number of metabolic syndrome components, waist circumference, C-reactive protein and tumor necrosis factor-α. Decreased appetite was negatively associated with BMI and waist circumference. Psychomotor retardation was positively associated with BMI, high-density lipoprotein cholesterol and triglycerides, and insomnia with number of metabolic syndrome components. CONCLUSION: Increased appetite - in the context of a depressive episode - was the only symptom that was associated with both metabolic as well as inflammatory markers, and could be a key feature of an immuno-metabolic form of depression. This immuno-metabolic depression should be considered in clinical trials evaluating effectiveness of compounds targeting metabolic and inflammatory pathways or lifestyle interventions.

Sit, step, sweat: longitudinal associations between physical activity patterns, anxiety and depression.

BACKGROUND: Physical inactivity has been identified as a risk factor for depression and, less often, as a long-term consequence of depression. Underexplored is whether similar bi-directional longitudinal relationships are observed for anxiety disorders, particularly in relation to three distinct indicators of activity levels - sports participation, general physical activity and sedentary behavior. METHOD: Participants were from the Netherlands Study of Depression and Anxiety (NESDA; N = 2932, 18-65 years old; 57% current anxiety or depressive disorder, 21% remitted disorder, 22% healthy controls). At baseline, 2, 4, and 6 years, participants completed a diagnostic interview and self-report questionnaires assessing psychopathology symptom severity, physical activity indicators, and sociodemographic and health covariates. RESULTS: Consistently across assessment waves, people with anxiety and/or depressive disorders had lower sports participation and general physical activity compared to healthy controls. Greater anxiety or depressive symptoms were associated with lower activity according to all three indicators. Over time, a diagnosis or greater symptom severity at one assessment was associated with poorer sports participation and general physical activity 2 years later. In the opposite direction, only low sports participation was associated with greater symptom severity and increased odds of disorder onset 2 years later. Stronger effects were observed for chronicity, with lower activity according to all indicators increasing the odds of disorder chronicity after 2 years. CONCLUSIONS: Over time, there seems to a mutually reinforcing, bidirectional relationship between psychopathology and lower physical activity, particularly low sports participation. People with anxiety are as adversely affected as those with depression.

Polygenic dissection of major depression clinical heterogeneity.

The molecular mechanisms underlying major depressive disorder (MDD) are largely unknown. Limited success of previous genetics studies may be attributable to heterogeneity of MDD, aggregating biologically different subtypes. We examined the polygenic features of MDD and two common clinical subtypes (typical and atypical) defined by symptom profiles in a large sample of adults with established diagnoses. Data were from 1530 patients of the Netherlands Study of Depression and Anxiety (NESDA) and 1700 controls mainly from the Netherlands Twin Register (NTR). Diagnoses of MDD and its subtypes were based on DSM-IV symptoms. Genetic overlap of MDD and subtypes with psychiatric (MDD, bipolar disorder, schizophrenia) and metabolic (body mass index (BMI), C-reactive protein, triglycerides) traits was evaluated via genomic profile risk scores (GPRS) generated from meta-analysis results of large international consortia. Single nucleotide polymorphism (SNP)-heritability of MDD and subtypes was also estimated. MDD was associated with psychiatric GPRS, while no association was found for GPRS of metabolic traits. MDD subtypes had differential polygenic signatures: typical was strongly associated with schizophrenia GPRS (odds ratio (OR)=1.54, P=7.8e-9), while atypical was additionally associated with BMI (OR=1.29, P=2.7e-4) and triglycerides (OR=1.21, P=0.006) GPRS. Similar results were found when only the highly discriminatory symptoms of appetite/weight were used to define subtypes. SNP-heritability was 32% for MDD, 38% and 43% for subtypes with, respectively, decreased (typical) and increased (atypical) appetite/weight. In conclusion, MDD subtypes are characterized by partially distinct polygenic liabilities and may represent more homogeneous phenotypes. Disentangling MDD heterogeneity may help the psychiatric field moving forward in the search for molecular roots of depression.

Gene expression in major depressive disorder.

The search for genetic variants underlying major depressive disorder (MDD) has not yet provided firm leads to its underlying molecular biology. A complementary approach is to study gene expression in relation to MDD. We measured gene expression in peripheral blood from 1848 subjects from The Netherlands Study of Depression and Anxiety. Subjects were divided into current MDD (N=882), remitted MDD (N=635) and control (N=331) groups. MDD status and gene expression were measured again 2 years later in 414 subjects. The strongest gene expression differences were between the current MDD and control groups (129 genes at false-discovery rate, FDR<0.1). Gene expression differences across MDD status were largely unrelated to antidepressant use, inflammatory status and blood cell counts. Genes associated with MDD were enriched for interleukin-6 (IL-6)-signaling and natural killer (NK) cell pathways. We identified 13 gene expression clusters with specific clusters enriched for genes involved in NK cell activation (downregulated in current MDD, FDR=5.8 × 10(-5)) and IL-6 pathways (upregulated in current MDD, FDR=3.2 × 10(-3)). Longitudinal analyses largely confirmed results observed in the cross-sectional data. Comparisons of gene expression results to the Psychiatric Genomics Consortium (PGC) MDD genome-wide association study results revealed overlap with DVL3. In conclusion, multiple gene expression associations with MDD were identified and suggest a measurable impact of current MDD state on gene expression. Identified genes and gene clusters are enriched with immune pathways previously associated with the etiology of MDD, in line with the immune suppression and immune activation hypothesis of MDD.

Meta-analysis of genome-wide association studies of anxiety disorders.

Anxiety disorders (ADs), namely generalized AD, panic disorder and phobias, are common, etiologically complex conditions with a partially genetic basis. Despite differing on diagnostic definitions based on clinical presentation, ADs likely represent various expressions of an underlying common diathesis of abnormal regulation of basic threat-response systems. We conducted genome-wide association analyses in nine samples of European ancestry from seven large, independent studies. To identify genetic variants contributing to genetic susceptibility shared across interview-generated DSM-based ADs, we applied two phenotypic approaches: (1) comparisons between categorical AD cases and supernormal controls, and (2) quantitative phenotypic factor scores (FS) derived from a multivariate analysis combining information across the clinical phenotypes. We used logistic and linear regression, respectively, to analyze the association between these phenotypes and genome-wide single nucleotide polymorphisms. Meta-analysis for each phenotype combined results across the nine samples for over 18 000 unrelated individuals. Each meta-analysis identified a different genome-wide significant region, with the following markers showing the strongest association: for case-control contrasts, rs1709393 located in an uncharacterized non-coding RNA locus on chromosomal band 3q12.3 (P=1.65 × 10(-8)); for FS, rs1067327 within CAMKMT encoding the calmodulin-lysine N-methyltransferase on chromosomal band 2p21 (P=2.86 × 10(-9)). Independent replication and further exploration of these findings are needed to more fully understand the role of these variants in risk and expression of ADs.

A common polymorphism in the ABCB1 gene is associated with side effects of PGP-dependent antidepressants in a large naturalistic Dutch cohort.

The drug efflux transporter permeability glycoprotein (PGP) and cytochrome P450 (CYP) 2C19 are important for eliminating antidepressants from the brain and body. The ABCB1 gene, encoding for PGP, and CYP2C19 gene have several variants that could influence enzyme function and thereby the effect of PGP- and 2C19-dependent antidepressants. We investigated the association of antidepressant side effect and common genetic variation in 789 antidepressant users. In PGP-dependent antidepressant users, the A-allele of the rs2032588 single-nucleotide polymorphism (SNP) was associated with a lower number of side effects after adjusting for gender, age, dosage and duration of use, (B=-0.44, q=4.6 × 10(-3)). This association was different from and absent in non-PGP-dependent antidepressant users. Other SNP associations as well as an interaction analysis between the rs2032588 SNP and the CYP2C19 SNPs were not statistically significant after adjusting for covariates and multiple comparisons. The association of rs2032588 with antidepressant side effects suggests the involvement of the ABCB1 genotype in the clinical pharmacology of PGP-dependent antidepressants.

Depressive and anxiety disorders and short leukocyte telomere length: mediating effects of metabolic stress and lifestyle factors.

BACKGROUND: Depressive and anxiety disorders are associated with shorter leukocyte telomere length (LTL), an indicator of cellular aging. It is, however, unknown which pathways underlie this association. This study examined the extent to which lifestyle factors and physiological changes such as inflammatory or metabolic alterations mediate the relationship. METHOD: We applied mediation analysis techniques to data from 2750 participants of the Netherlands Study of Depression and Anxiety. LTL was assessed using quantitative polymerase chain reaction. Independent variables were current depressive (30-item Inventory of Depressive Symptoms - Self Report) and anxiety (21-item Beck's Anxiety Inventory) symptoms and presence of a depressive or anxiety disorder diagnosis based on DSM-IV; mediator variables included physiological stress systems, metabolic syndrome components and lifestyle factors. RESULTS: Short LTL was associated with higher symptom severity (B = -2.4, p = 0.002) and current psychiatric diagnosis (B = -63.3, p = 0.024). C-reactive protein, interleukin-6, waist circumference, triglycerides, high-density lipoprotein cholesterol and cigarette smoking were significant mediators in the relationship between psychopathology and LTL. When all significant mediators were included in one model, the effect sizes of the relationships between LTL and symptom severity and current diagnosis were reduced by 36.7 and 32.7%, respectively, and the remaining direct effects were no longer significant. CONCLUSIONS: Pro-inflammatory cytokines, metabolic alterations and cigarette smoking are important mediators of the association between depressive and anxiety disorders and LTL. This calls for future research on intervention programs that take into account lifestyle changes in mental health care settings.

The association between lower educational attainment and depression owing to shared genetic effects? Results in ~25,000 subjects.

An association between lower educational attainment (EA) and an increased risk for depression has been confirmed in various western countries. This study examines whether pleiotropic genetic effects contribute to this association. Therefore, data were analyzed from a total of 9662 major depressive disorder (MDD) cases and 14,949 controls (with no lifetime MDD diagnosis) from the Psychiatric Genomics Consortium with additional Dutch and Estonian data. The association of EA and MDD was assessed with logistic regression in 15,138 individuals indicating a significantly negative association in our sample with an odds ratio for MDD 0.78 (0.75-0.82) per standard deviation increase in EA. With data of 884,105 autosomal common single-nucleotide polymorphisms (SNPs), three methods were applied to test for pleiotropy between MDD and EA: (i) genetic profile risk scores (GPRS) derived from training data for EA (independent meta-analysis on ~120,000 subjects) and MDD (using a 10-fold leave-one-out procedure in the current sample), (ii) bivariate genomic-relationship-matrix restricted maximum likelihood (GREML) and (iii) SNP effect concordance analysis (SECA). With these methods, we found (i) that the EA-GPRS did not predict MDD status, and MDD-GPRS did not predict EA, (ii) a weak negative genetic correlation with bivariate GREML analyses, but this correlation was not consistently significant, (iii) no evidence for concordance of MDD and EA SNP effects with SECA analysis. To conclude, our study confirms an association of lower EA and MDD risk, but this association was not because of measurable pleiotropic genetic effects, which suggests that environmental factors could be involved, for example, socioeconomic status.

The association between low vitamin D and depressive disorders.

It has been hypothesized that hypovitaminosis D is associated with depression but epidemiological evidence is limited. We investigated the association between depressive disorders and related clinical characteristics with blood concentrations of 25-hydroxyvitamin D [25(OH)D] in a large cohort. The sample consisted of participants (aged 18-65 years) from the Netherlands Study of Depression and Anxiety (NESDA) with a current (N=1102) or remitted (N=790) depressive disorder (major depressive disorder, dysthymia) defined according to DSM-IV criteria, and healthy controls (N=494). Serum levels of 25(OH)D measured and analyzed in multivariate analyses adjusting for sociodemographics, sunlight, urbanization, lifestyle and health. Of the sample, 33.6% had deficient or insufficient serum 25(OH)D (<50 nmol l(-1)). As compared with controls, lower 25(OH)D levels were found in participants with current depression (P=0.001, Cohen's d=0.21), particularly in those with the most severe symptoms (P=0.001, Cohen's d=0.44). In currently depressed persons, 25(OH)D was inversely associated with symptom severity (ß=-0.19, s.e.=0.07, P=0.003) suggesting a dose-response gradient, and with risk (relative risk=0.90, 95% confidence interval=0.82-0.99, P=0.03) of having a depressive disorders at 2-year follow-up. This large cohort study indicates that low levels of 25(OH)D were associated to the presence and severity of depressive disorder suggesting that hypovitaminosis D may represent an underlying biological vulnerability for depression. Future studies should elucidate whether-the highly prevalent-hypovitaminosis D could be cost-effectively treated as part of preventive or treatment interventions for depression.

Relationship between vitamin D status and left ventricular geometry in a healthy population: results from the Baltimore Longitudinal Study of Aging.

OBJECTIVES: The effects of vitamin D on the heart have been studied in patients with cardiac disease, but not in healthy persons. We investigated the relation between vitamin D status and left ventricular (LV) structure and function in community-dwelling subjects without heart disease. DESIGN: The relationship between concentrations of 25-hydroxyvitamin D [25(OH)D], a marker of vitamin D reserve, and LV transthoracic echocardiography measures was analysed in 711 participants in the Baltimore Longitudinal Study of Aging who were without cardiac disease. RESULTS: Mean 25(OH)D in the study population was 32.3 ± 11.4 ng mL(-1) ; only 15.5% of subjects had moderate or severe vitamin D deficiency [25(OH)D < 20 ng mL(-1) ]. Adjusting for age, body mass index, cardiovascular disease risk factors, physical activity, calcium and parathyroid hormone, 25(OH)D was positively correlated with LV thickness (ß 0.095, SE 0.039, P < 0.05) and LV mass index (ß 7.5, SE 2.6, P < 0.01). A significant nonlinear relation between 25(OH)D and LV concentric remodelling was observed. LV remodelling was more likely in participants with 25(OH)D levels <30 ng mL(-1) [odds ratio (OR) 1.24; 95% confidence interval (CI) 0.83-1.85] or ≥38 ng mL(-1) (OR 1.73; 95% CI 1.13-2.65), compared with those with 30-37 ng mL(-1) 25(OH)D. Consistently, LV relative wall thickness was significantly lower (P for trend=0.05), and LV diastolic internal diameter index (P for trend<0.05) and end-diastolic volume index (P for trend<0.05) were significantly higher in subjects with 30-37 ng mL(-1) 25(OH)D compared to the rest of the study population. There was a significant interaction between 25(OH)D and hypertension on the risk of LV hypertrophy (P < 0.05). CONCLUSIONS: In a population-based sample of predominantly vitamin D-sufficient subjects without heart disease, LV geometry was most favourable at intermediate 25(OH)D concentrations.

The relationship between visceral adiposity and left ventricular diastolic function: results from the Baltimore Longitudinal Study of Aging.

BACKGROUND AND AIMS: It is unclear whether subcutaneous and visceral fat are differentially correlated to the decline in left ventricular (LV) diastolic function with aging. This study sought to examine the hypothesis that age-related changes in the regional fat distribution account for changes in LV diastolic function and to explore potential mediators of this association. METHODS AND RESULTS: In this cross-sectional study, we evaluated 843 participants of the Baltimore Longitudinal Study of Aging with echocardiogram, dual-energy X-ray absorptiometry (DEXA), abdominal computed tomography (CT) and blood tests performed at the same visit. LV diastolic function was assessed by parameters of LV relaxation (E/A ratio, Em and Em/Am ratio) and LV filling pressures (E/Em ratio). Total body fat was computed by DEXA, while visceral and subcutaneous fat were determined from abdominal CT. In multivariate models adjusted for demographics, cardiovascular risk factors, antihypertensive medications, physical activity and LV mass, both visceral and subcutaneous fat were associated with LV diastolic dysfunction. When both measures of adiposity were simultaneously included in the same model, only visceral fat was significantly associated with LV diastolic dysfunction. Triglycerides and sex-hormone binding globulin, but not adiponectin and leptin, were found to be significant mediators of the relationship between visceral fat and LV diastolic function, explaining 28-47% of the association. Bootstrapping analyses confirmed the significance of these findings. CONCLUSIONS: Increased visceral adiposity is associated with LV diastolic dysfunction, possibly through a metabolic pathway involving blood lipids and ectopic fat accumulation rather than adipokines.

Efficacy of celecoxib add-on treatment for immuno-metabolic depression: Protocol of the INFLAMED double-blind placebo-controlled randomized controlled trial.

Introduction: As the role of (neuro)inflammation in depression pathophysiology is emerging, augmentation of antidepressant treatments with anti-inflammatory drugs have shown beneficial results, but not consistently across all studies. Inconsistencies may be due to depression biological and clinical heterogeneity. Immuno-Metabolic Depression (IMD) has been put forward as a form of depression characterized by the clustering of low-grade inflammation, metabolic dysregulations and atypical, energy-related symptoms (overeating, weight gain, hypersomnia, fatigue and leaden paralysis). IMD features are present in ∼30% of patients with Major Depressive Disorder (MDD). By selecting these specific patients, directly targeting inflammation may reduce depressive symptoms. Methods: and analysis INFLAMED is a double-blind randomized controlled trial. 140 MDD patients with IMD characteristics (MDD with Inventory of Depressive Symptomatology (IDS) ≥ 26, IDS atypical, energy related symptoms ≥6, C-Reactive Protein (CRP) > 1 mg/L) will receive either 400 mg celecoxib per day or matching placebo for a period of 12 weeks. Biological, physical and interview data will be collected after 2, 6 and 12 weeks of starting the intervention. Questionnaires will be sent out bi-weekly during the study period. The main study outcome is the IDS (30-item self-report) total score during 12-week follow-up. Secondary study outcomes include response, remission, adverse side effects, symptom profiles (atypical, energy-related symptoms), fatigue, food craving, sleep, anxiety symptoms, functioning, pain, and optionally, microbiome composition. Explorative analyses will be performed on the role of CRP, IL-6, TNF-α, cholesterol, triglycerides, glucose, BMI, waist and hip circumference. Ethics and dissemination: This protocol has been approved by the Medical Ethics Review Board of the Amsterdam UMC, location VUmc (2022.0015) on 2-6-2022, as well as by the competent authority in The Netherlands: CCMO, on 3-8-2022. Registration details: Trail registration numbers NCT05415397, EudraCT 2021-003850-21.

Pathophysiology of major depressive disorder: mechanisms involved in etiology are not associated with clinical progression.

Meta-analyses support the involvement of different pathophysiological mechanisms (inflammation, hypothalamic-pituitary (HPA)-axis, neurotrophic growth and vitamin D) in major depressive disorder (MDD). However, it remains unknown whether dysregulations in these mechanisms are more pronounced when MDD progresses toward multiple episodes and/or chronicity. We hypothesized that four central pathophysiological mechanisms of MDD are not only involved in etiology, but also associated with clinical disease progression. Therefore, we expected to find increasingly more dysregulation across consecutive stages of MDD progression. The sample from the Netherlands Study of Depression and Anxiety (18-65 years) consisted of 230 controls and 2333 participants assigned to a clinical staging model categorizing MDD in eight stages (0, 1A, 1B, 2, 3A, 3B, 3C and 4), from familial risk at MDD (stage 0) to chronic MDD (stage 4). Analyses of covariance examined whether pathophysiological mechanism markers (interleukin (IL)-6, C-reactive protein (CRP), cortisol, brain-derived neurotrophic factor and vitamin D) showed a linear trend across controls, those at risk for MDD (stages 0, 1A and 1B), and those with full-threshold MDD (stages 2, 3A, 3B, 3C and 4). Subsequently, pathophysiological differences across separate stages within those at risk and with full-threshold MDD were examined. A linear increase of inflammatory markers (CRP P=0.026; IL-6 P=0.090), cortisol (P=0.025) and decrease of vitamin D (P<0.001) was found across the entire sample (for example, from controls to those at risk and those with full-threshold MDD). Significant trends of dysregulations across stages were present in analyses focusing on at-risk individuals (IL-6 P=0.050; cortisol P=0.008; vitamin D P<0.001); however, no linear trends were found in dysregulations for any of the mechanisms across more progressive stages of full-threshold MDD. Our results support that the examined pathophysiological mechanisms are involved in MDD's etiology. These same mechanisms, however, are less important in clinical progression from first to later MDD episodes and toward chronicity.

Determinants of plasma 25-hydroxyvitamin D levels in healthy adults in the Netherlands.

BACKGROUND: Vitamin D plays a key role in maintaining skeletal health, but is also related to various non-skeletal health issues. Several determinants have been identified that influence blood plasma levels of 25-hydroxyvitamin D (25(OH) D), often in specific patients or elderly populations. This paper aims to replicate these findings in a healthy population. METHODS: Plasma levels of 25(OH)D were measured using tandem mass spectrometry. We examined the cross-sectional association of sociodemographic, health, lifestyle and sampling characteristics with 25(OH)D in a group of 539 adults, who were healthy control subjects in the NESDA study in the Netherlands (latitude 52 °N). RESULTS: Mean 25(OH)D levels were 68.0 (± 27.2) nmol/l. Levels under 50 nmol/l occurred in 27% of the population; 40% reached levels above 75 nmol/l. Women had higher levels than men, and the use of oral contraceptives showed a significant positive association among females. Subjects with non-European ancestry had dramatically lower 25(OH) D levels. Other factors that were negatively associated were body mass index and the renal estimated glomerular filtration rate (eGFR). Meteorological data replaced season as a significant determinant. Moderate alcohol consumption and sports showed a positive association, while physical activity and the hepatic marker gamma-glutamyl transferase did not. Our results disconfirm the influence of age in this population of under 65 year olds. CONCLUSION: Insufficient 25(OH)D levels were common in a healthy population. The set of eight variables that were significant in a multiple regression model (sex, ancestry, oral contraceptives, eGFR, BMI, sports, alcohol, sunshine) explained 29.5% of the variance.