The association between ambient UVB dose and ANCA-associated vasculitis relapse and onset
Arthritis Research & Therapy volume 24, Article number: 147 (2022)
The aetiology of ANCA-associated vasculitis (AAV) and triggers of relapse are poorly understood. Vitamin D (vitD) is an important immunomodulator, potentially responsible for the observed latitudinal differences between granulomatous and non-granulomatous AAV phenotypes. A narrow ultraviolet B spectrum induces vitD synthesis (vitD-UVB) via the skin. We hypothesised that prolonged periods of low ambient UVB (and by extension vitD deficiency) are associated with the granulomatous form of the disease and an increased risk of AAV relapse.
Patients with AAV recruited to the Irish Rare Kidney Disease (RKD) (n = 439) and UKIVAS (n = 1961) registries were studied. Exposure variables comprised latitude and measures of ambient vitD-UVB, including cumulative weighted UVB dose (CW-D-UVB), a well-validated vitD proxy. An n-of-1 study design was used to examine the relapse risk using only the RKD dataset. Multi-level models and logistic regression were used to examine the effect of predictors on AAV relapse risk, phenotype and serotype.
Residential latitude was positively correlated (OR 1.41, 95% CI 1.14–1.74, p = 0.002) and average vitD-UVB negatively correlated (0.82, 0.70–0.99, p = 0.04) with relapse risk, with a stronger effect when restricting to winter measurements (0.71, 0.57–0.89, p = 0.002). However, these associations were not restricted to granulomatous phenotypes. We observed no clear relationship between latitude, vitD-UVB or CW-D-UVB and AAV phenotype or serotype.
Our findings suggest that low winter ambient UVB and prolonged vitD status contribute to AAV relapse risk across all phenotypes. However, the development of a granulomatous phenotype does not appear to be directly vitD-mediated. Further research is needed to determine whether sufficient vitD status would reduce relapse propensity in AAV.
Anti-neutrophil cytoplasm antibody (ANCA)-associated vasculitis (AAV) is a rare multi-system autoimmune disorder, with an annual incidence of 13–20 cases per million .
It manifests as one of three main phenotypes: granulomatosis with polyangiitis (GPA), microscopic polyangiitis (MPA) or eosinophilic granulomatosis with polyangiitis (EGPA). Characteristic autoantibodies directed against myeloperoxidase (MPO) or proteinase-3 (PR3) molecules determine the ANCA serotype. It follows a relapsing and remitting course, resulting in cumulative tissue damage. While the precise triggers of AAV onset and relapse are unknown, evidence supports a complex interaction of genetic susceptibility loci , epigenetics  and environmental factors .
A potential effect of ultraviolet (UV) radiation on disease phenotype and disease activity has been proposed, as marked variations in GPA:MPA incidence ratios occur with latitude [5,6,7,8]. Increasing latitude (and hence decreasing UV)  is associated with an increased incidence of GPA and EGPA subtypes [5, 8] and of anti-PR3 positivity . The magnitude of association may be strongest for winter UV levels , because it is these that particularly limit cutaneous vitamin D (vitD) synthesis  and the strongest latitudinal gradient in UV is observed in winter . A similar inverse correlation between UV radiation exposure and incidence is observed in other autoimmune diseases (ADs), for example, multiple sclerosis [12, 13] and type 1 diabetes . In contrast, for the MPA phenotype, an increasing incidence with decreasing latitude was observed [5, 6], although Gatenby et al. could not replicate this relationship (nor did they identify an association with UV radiation) .
vitD has immunomodulatory effects, including suppression of Th1 and Th17 cells and upregulation of Treg cells, resulting in the reduction of inflammatory cytokines and a more tolerogenic environment [4, 14,15,16]. Th1, Th17 and Treg cells are crucial to granuloma formation, characteristically seen in GPA and EGPA, but not in MPA, where Th2 cells and antibody-mediated mechanisms predominate . Gatenby et al. postulated that this vitD-mediated UV effect on granulomas was responsible for the strong inverse relationship between GPA (but not MPA) incidence and UV radiation, but this hypothesis has not been tested.
Several small studies demonstrate that serum 25-hydroxyvitamin D (25OHD) levels are lower in AAV patients than in healthy controls [18,19,20]. Low 25OHD levels have been associated with increased disease activity in other ADs [21,22,23], with suggestive data in AAV . To our knowledge, there are no large studies examining the influence of UV radiation or vitD on AAV disease activity, as distinct from ANCA phenotype or serotype. Ambient UVB dose at wavelengths that induce cutaneous vitD3 synthesis (vitD-UVB) is the key determinant of an individual’s vitD status. Repeated 25OHD assessment is costly and not routinely performed. vitD levels also drop with inflammation limiting the value of 25OHD measurement at relapse. Therefore, we chose to use the cumulative-weighted ambient UVB dose (CW-D-UVB)—a well-validated predictor of an individual’s vitD status [24,25,26,27]—as a surrogate.
The aim of this study is, firstly, to examine the association between measures of ambient UV radiation, UVB-predicted vitD status (CW-D-UVB) and AAV relapse, using a novel n-of-1 case-control design. Secondly, we aim to explore the relationship between these variables and ANCA phenotype and serotype at diagnosis. We hypothesise that prolonged periods of low ambient UVB (and by extension vitD deficiency) are associated with an increased risk of AAV relapse in GPA/EGPA, but not MPA, and an increased risk of a GPA/EGPA phenotype at diagnosis.
UKIVAS cohort 1 comprised patients from the United Kingdom and Ireland Vasculitis Rare Disease Group (UKIVAS) Registry (https://ukivas.ndorms.ox.ac.uk/), established in 2009. Patients were included in this study if they had a diagnosis of definite AAV with documented diagnosis date and location, defined by relevant clinical features of GPA, MPA or EGPA , with either positive anti-MPO or PR3 serology  and/or diagnostic histopathology. Patients with secondary vasculitis or dual anti-glomerular basement membrane disease were excluded .
RKD cohort 2 comprised non-overlapping patients from the national Irish Rare Kidney Disease (RKD) Registry and Biobank (https://www.tcd.ie/medicine/thkc/research/rare.php), established in 2012. The same inclusion and exclusion criteria were applied. For the AAV relapse component, participants were also required to have > 1 documented period of remission prior to the commencement of the analysis on 18 August 2020. Participant selection is outlined in Additional file 1: Figs. S1 and S2. The AAV diagnosis study included both cohorts, while the AAV relapse study included RKD only, as granular longitudinal follow-up was not recorded in UKIVAS. Ethical approval for the RKD and UKIVAS components was granted by the Tallaght University Research Ethics Committee (ref 2019-08 List 29 (07)) and the NHS National Research Ethics Service (REC 10/H1102/77), respectively. All participants provided written informed consent.
Additional file 1: Fig. S3 shows the participant locations at diagnosis. We employed a smartphone application for the RKD cohort (patientMpower)  to capture participants’ location by GPS . If unavailable at the time of assessment, the electoral division of the participant’s residence or their hospital was used (ranked in preference order). Data linkage was performed using a Resource Description Framework (RDF) semantic web model . Participants’ residential postcodes informed the location coordinates for the UKIVAS cohort.
Ultraviolet B (UVB) data source, measures of ambient vitD-UVB and cumulative-weighted UVB dose (CW-D-UVB)
Daily UVB data at wavelengths specific for vitD production (vitD-UVB) were obtained from the Tropospheric Emission Monitoring Internet Service (TEMIS) database [34, 35]. To generate UVB-predicted vitamin D status, a specific cumulative-weighted UVB dose (CW-D-UVB) was calculated for each participant, determined by high-resolution location and date, based on an algorithm proposed by Kelly et al.  and later validated [25,26,27]. Briefly, a daily vitD-UVB dose over 135 days before a given date was used with a decay function assuming a vitD half-life of 35 days (Fig. 1). Further detail is provided in the supplementary materials.
The average (2004–2019) annual, average winter and preceding winter vitD-UVB and CW-D-UVB were derived for each participant. Winter measures were specifically included given the stronger effect sizes previously noted .
The study was split into AAV relapse and diagnosis components, further described in Fig. 2 and supplementary materials. Briefly, relapse was studied using an n-of-1 design where each participant acts as their own control, eliminating confounding by time-independent factors such as gender, skin type and occupation and genotype. Case (relapse) and control (remission) windows were defined using an algorithm that considered the effect of residual disease activity and the CW-D-UVB interval (135 days). The diagnosis was studied using the estimated date of symptom onset. Relapse was defined as the return of signs and/or symptoms of active vasculitis (BVAS > 0), with supporting laboratory or histopathological evidence, escalation in immunosuppression and a clinical response to the same. Relapse events were adjudicated by a committee (JS and ML), blinded to the individual’s UV data. Remission was defined as the absence of signs, symptoms and laboratory evidence of vasculitis activity. Follow-up continued until the last documented visit or death.
Patient and Public Involvement
Patient and Public Involvement is detailed in supplementary materials.
A multi-level model (lme4 package, glmer function)  was applied to investigate the association between latitude/measures of ambient vitD-UVB/CW-D-UVB and AAV relapse. Covariates included age, gender, AAV phenotype, ANCA serotype and treatment. A random effect was included to account for repeated measures and the varying relapse risk between individuals . The disease phenotype was included as an interaction term, to study the (hypothesised) effect of exposure variables on disease subtype relapse risk. Sensitivity analysis restricted to (non-) granulomatous phenotypes was also performed. The power calculation is detailed in supplementary materials.
Logistic regression was used to examine the relationships between AAV phenotype and serotype (outcome variables), latitude, measures of ambient vitD-UVB and CW-D-UVB (examined separately due to collinearity). Covariates included age, gender and ethnicity. We performed a sensitivity analysis restricted to Whites only, given their more substantial vitD synthesis in the skin.
Odds ratios (ORs) and 95% confidence intervals (CIs) were reported. p < 0.05 was considered statistically significant. All statistical analyses were performed using R (version 4.0.4).
UKIVAS cohort 1 (N = 1961, Table 1) and RKD cohort 2 (N = 439, Table 2) were similar. An excess of non-White and GPA participants in UKIVAS were the main differences. Induction and maintenance treatment are detailed in Additional file 1: Table S1. The Median follow-up time in RKD was 58.3 months (IQR 32.1–138.1). A total of 343 (78.1%) patients did not relapse; there were 135 relapse events in 96 patients.
The annual peak and nadir of CW-D-UVB were observed in August and February, respectively, lagging behind those of vitD-UVB by 2 months thus mimicking 25OHD seasonal fluctuations (Fig. 1 and Additional file 1: Fig. S4). The daily average vitD-UVB and CW-D-UVB were strongly correlated with latitude (Additional file 1: Fig. S5a-b), as expected. However, latitude alone was weakly correlated with vitD-UVB or CW-vitD-UVB on a particular day (e.g. diagnosis), because the impact of the season is much greater than latitude. This correlation remained weak even when stratified by season due to large variation in daily UVB (Additional file 1: Fig. S5c-g), thus strengthening the rationale for using CW-D-UVB as vitD proxy, over latitude or season alone.
Impact of latitude, measures of ambient vitD-UVB and cumulative-weighted UVB (CW-D-UVB) on AAV relapse
Using multi-level logistic regression analysis, latitude was associated with AAV relapse risk (OR 1.41, 95% CI 1.14–1.74, Table 3 and Fig. 3a). There was also a significant inverse relationship between relapse risk and average winter vitD-UVB (0.71, 0.57–0.89, Figs. 3b and 4) and annual vitD-UVB (0.82, 0.70–0.99, Fig. 3d), with a slightly larger effect size in MPA versus not-MPA (Additional file 1: Table S2). However, no relationship was observed with preceding winter vitD-UVB (0.90, 0.74–1.10, Table 3 and Additional file 1: Fig. S6b). Similarly, we noted a significant inverse relationship between relapse risk and both average (0.82, 0.70–0.99, Fig. 3c) and preceding (0.81, 0.66–0.99, Additional file 1: Fig. S6c) winter CW-D-UVB. The median CW-D-UVB at symptom onset was not significantly associated with relapse, overall or when stratified by AAV phenotype and ANCA serotype (1.06, 0.88–1.28, Additional file 1: Fig. S6a). The addition of an interaction term between AAV phenotype and any exposure variable did not alter the findings. Additional independent risk factors for relapse included GPA/EGPA subtype (1.78, 1.03–3.05), younger age (0.75, 0.61–0.92) and being off immunosuppression therapy (2.65, 1.70–4.14) or on glucocorticoid monotherapy (1.85, 1.02–3.35, Additional file 1: Table S3).
Impact of latitude, measures of ambient vitD-UVB and cumulative-weighted UVB (CW-D-UVB) on AAV phenotype and serotype at diagnosis
We then investigated the effect of UVB exposure on AAV phenotype and serotype. Approximately half of the cases were diagnosed in winter/spring, with no demonstrable seasonal effect (Tables 1 and 2). UKIVAS and RKD participants reside at a relatively narrow latitudinal range, between 49.4–60.3° and 51.9–55.3°, respectively. Analyses of potential predictors of AAV phenotype and serotype are summarised in Tables 4 and Additional file 1: Table S4, respectively, and in Additional file 1: Figs. S7 and S8. We observed no significant associations between latitude (adjusted p = 0.39), average annual vitD-UVB (p = 0.65), average winter vitD-UVB (p = 0.46) or CW-D-UVB at symptom onset (p = 0.66) and being diagnosed with MPA or being MPO-ANCA positive. The findings were consistent in the sensitivity analysis when restricted to White participants (Additional file 1: Table S5). While non-significant, we found that vitD-UVB status at diagnosis (CW-D-UVB) was consistently below the annual average at the participants’ location (Additional file 1: Fig. S9).
For the first time, we examined the relationships between latitude, measures of ambient vitD-UVB and UVB-predicted vitD status and AAV relapse risk. In keeping with other ADs, we demonstrated that low ambient UVB and, by extension, vitD deficiency are associated with AAV relapse, in a genetically homogeneous Irish cohort, after adjustment for key confounders. However, in this large registry-based study that included 2400 AAV patients from an ethnically and latitudinally restricted Irish and British cohort, we found no significant association between latitude, measures of ambient vitD-UVB or UVB-predicted vitD status (CW-D-UVB) and AAV phenotype or ANCA serotype.
While vitD deficiency has consistently been associated with increased disease activity in other ADs [21, 22], it has not previously been explored in depth in AAV. The risk of AAV relapse was positively associated with residential latitude and negatively with average vitD-UVB. The magnitude of the effect was greatest for average winter vitD-UVB, akin to the findings of Gatenby et al. . This is biologically plausible as low winter exposure to vitD-UVB limits cutaneous vitD synthesis, particularly at more northerly latitudes . We observed a similar negative correlation between average winter CW-D-UVB and relapse risk. For each 1 kJ/m2 increase in average winter CW-D-UVB, which equates to a 27.2 nmol/L rise in serum 25OHD , the odds of relapse decreased by 26%. Clinically, this is similar to moving from vitD ‘deficiency’ (< 25 nmol/L) to ‘sufficiency’ (> 50 nmol/L) .
We then explored narrower time intervals to pinpoint the at-risk window. While the average CW-D-UVB over the preceding winter (to the event) remained independently predictive of AAV relapse, the average vitD-UVB over the same period was not. The vitD-UVB nadir is in December, while CW-D-UVB (and serum 25OHD levels) continues to fall over ‘winter’, with a nadir 2 months later, so the effect of the latter may thus have been lost by restricting the analysis to winter. Importantly, we found no relationship between UVB-predicted vitD status (CW-D-UVB) at relapse onset and relapse risk. Together, these suggest that, an individual’s vitD status over a more prolonged period, particularly during winter, rather than their status peri-relapse, renders them susceptible to relapse and that multiple periods of low UVB and vitD deficiency might be contributing towards disease susceptibility. There was no significant interaction between AAV phenotype and any exploratory variable, suggesting there was no phenotype-specific effect, contrary to our hypothesis.
Drawing on incidence studies from 10 countries, Gatenby et al.  observed that GPA (but not MPA) incidence positively and negatively correlated with latitude (35.30° S to 70° N) and ambient UV radiation, respectively . Consistent with this and other series [5,6,7], we observed a subtle (non-significant) increased odds of MPA (vs. GPA/EGPA) with decreasing latitude and increasing average vitD-UVB. It is possible that this diagnosis study component, using standard logistic regression, was underpowered due to the narrow latitudinal range in our cohort. However, contrary to suggestions that phenotype-specific UV effects are mediated by vitD effects on granuloma formation , we observed no relationship between UVB-predicted vitD status at symptom onset and AAV phenotype. As with AAV relapse, it may be an individual’s vitD status over a longer period, rather than immediately preceding disease onset, that is important. Critically, the outcome was MPA (vs. non-MPA) in this diagnosis study component. In our case-only cohorts, we could not directly interrogate whether UV radiation or UVB-predicted vitD status are implicated in overall AAV risk (i.e. AAV vs. non-AAV control). It is likely that vitD is important in de novo AAV onset but may not be responsible for defining the (non-) granulomatous phenotype.
The classification of AAV by ANCA serotype is increasingly preferred due to its enhanced predictive value in terms of relapse and renal survival [43, 44]. We found no significant association between latitude, average vitD-UVB or CW-D-UVB and ANCA serotype at diagnosis, although the direction of trends between latitude (negative)/vitD-UVB (positive) and MPO positivity was similar to that of MPA. This diverges from a multi-centre AAV glomerulonephritis registry-based cohort study that found anti-PR3 positivity was positively associated with latitude and negatively with UV radiation . This study was comparable in size (n = 1408) to our cohort but covered a wider latitudinal range (Europe and the USA: 35.8° N to 69.6° N) and included a greater ethnic mix. In keeping with our findings, the significant association between ANCA serotype and latitude/UV radiation was lost when the analysis was restricted to centres in North and Central Europe, focusing on a more genetically homogeneous population. The observed variation may therefore be due to genetic or other differences between populations (± their interaction with UVB) as an alternative explanation for reported geographic patterns, for example, the distribution of the frequency of the HLA-DPB*0401 allele .
Our study is the first to investigate the effect of CW-D-UVB, a precise and well-validated vitD proxy, in AAV at an individual level, with high spatial and temporal resolution. This contrasts with the less refined vitD proxies used previously at a population level, such as latitude, sunshine hours or rainfall [8, 9]. The AAV relapse component of our study is the first of its kind to explore the role of ambient measures of vitD-UVB and UVB-predicted vitD on AAV relapse. For this, we exploited the benefits of the novel n-of-1 study design, primed for use in rare diseases , by providing increased power to detect associations by expanding control events. The case-crossover design allowed us to account for many unmeasured determinants of vitD status, such as skin type, lifestyle or SPF use, by comparing UVB measures in active and remission states within the same individual. We were also able to leverage detailed longitudinal immunosuppressive treatment data, enabling us to account for this as a time-varying confounder.
The limitations of our study must also be considered. Our predominantly White cohort originates from a relatively narrow Irish/UK latitudinal gradient. Consequently, restricted latitudinal and genetic variability may contribute to a type II error. Using ambient UVB dose rather than vitD has advantages and disadvantages. While CW-D-UVB as a vitD proxy has been validated in the UK biobank population of almost half a million participants , a community-dwelling Irish cohort > 60 years  and a RCT of inflammatory bowel disease , it has not been validated in AAV patients. Other factors beyond UVB contribute to vitD status: e.g. vitD supplementation might have been high in our cohort (due to its use in renal insufficiency and frequent over-the-counter use following strong public health campaigns), or disease-associated inflammation might have depleted vitD, potentially weakening the link between UVB and vitD status. We did not have data on vitD supplementation, sun exposure, time spent outside, clothing, SPF use, angle from the sun and many other factors—many of which would be impossible to obtain. However, we assumed each participant’s habits were consistent in the application of the multilevel model and the n-of-1 study design controls for such within-individual factors. Sporadic measurement of serum 25OHD concentration, the best marker of vitD status is of limited value, because it is strongly affected by season and disease activity  (reverse causation) and is not representative of longer-term trajectories. Therefore, a prospective cohort study examining the effect of vitD status on AAV onset is impossible (due to disease rarity) and a case-control study would be biassed by reverse causality (inflammation decreases vitD). Therefore, our study design using CW-D-UVB is valuable in the context of these limitations.
In conclusion, we observed a strong inverse relationship between AAV relapse risk and individualised measures of vitD-UVB. A maximal effect was observed in winter. Prolonged low ambient UVB likely renders an individual vulnerable to a subsequent ‘hit’, which ultimately triggers disease activation. We found no evidence limiting this to granulomatous phenotypes. Our large study did not find any significant association between latitude, measures of vitD-UVB nor UVB-predicted vitamin D status and ANCA phenotype or serotype, albeit in an ethnically restricted cohort, across a narrow latitudinal window. It is plausible that prolonged vitD deficiency predisposes an individual to de novo AAV onset, but it does not appear to be important in the differentiation of AAV phenotype or serotype. Our observations have implications for the design of future vitD studies involving all relapsing-remitting diseases, not just AAV. We recommend a prospective cohort study to measure serial serum (or bloodspot) 25OHD at 3-monthly intervals to examine the trajectory of vitD in AAV relapse more accurately, focusing on average or winter vitD status. A RCT of vitD supplementation in AAV, targeting ‘sufficient’ 25OHD levels (>50nmol/L), with relapse as the primary outcome may then be warranted. vitD supplementation is generally well tolerated, safe and cheap to administer, but routine prescription should be avoided prior to definitive confirmation given the extra pill burden, unnecessary cost and potential for vitD toxicity with resultant hypercalcaemia.
Availability of data and materials
Due to the rarity of ANCA-associated vasculitis, coupled with the identifiable nature of the data, raw data must remain confidential and cannot be freely shared. We would invite any potential research collaborations or data requests through the corresponding author on reasonable request.
Anti-neutrophil cytoplasm antibody
Small vessel vasculitis
Granulomatosis with polyangiitis
Eosinophilic granulomatosis with polyangiitis
Ambient UVB dose at wavelengths that induce cutaneous vitD3 synthesis
Cumulative-weighted ambient UVB dose
United Kingdom and Ireland Vasculitis Rare Disease Group
- RKD Registry and Biobank:
Rare Kidney Disease Registry and Biobank
Resource Description Framework
Tropospheric Emission Monitoring Internet Service
- 95% CI:
95% confidence interval
Akaike Information Criterion
Not statistically significant
Randomised controlled trial
Watts RA, Mahr A, Mohammad AJ, Gatenby P, Basu N, Flores-Suárez LF. Classification, epidemiology and clinical subgrouping of antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis. Nephrol Dial Transplant. 2015;30:i14–22.
Lyons PA, Rayner TF, Trivedi S, Holle JU, Watts RA, Jayne DRW, et al. Genetically distinct subsets within ANCA-associated vasculitis. N Engl J Med. 2012;367(3):214–23.
Oates T, Salama AD. Epigenetic modifications in ANCA-associated vasculitis: potential for insights into disease pathogenesis and prediction of outcome? J Am Soc Nephrol. 2017;28(4):1011–3.
Scott J, Hartnett J, Mockler D, Little MA. Environmental risk factors associated with ANCA associated vasculitis: a systematic mapping review. Autoimmun Rev. 2020;19(11):102660.
Watts RA, Gonzalez-Gay MA, Lane SE, Garcia-Porrua C, Bentham G, Scott DG. Geoepidemiology of systemic vasculitis: comparison of the incidence in two regions of Europe. Ann Rheum Dis. 2001;60(2):170–2.
Watts RA, Lane SE, Scott DG, Koldingsnes W, Nossent H, Gonzalez-Gay MA, et al. Epidemiology of vasculitis in Europe. Ann Rheum Dis. 2001;60(12):1156–7.
O’Donnell JL, Stevanovic VR, Frampton C, Stamp LK, Chapman PT. Wegener’s granulomatosis in New Zealand: evidence for a latitude-dependent incidence gradient. Intern Med J. 2007;37(4):242–6.
Gatenby PA, Lucas RM, Engelsen O, Ponsonby AL, Clements M. Antineutrophil cytoplasmic antibody-associated vasculitides: could geographic patterns be explained by ambient ultraviolet radiation? Arthritis Rheum. 2009;61(10):1417–24.
Weiner M, Bjørneklett R, Hrusková Z, Mackinnon B, Poulton CJ, Sindelar L, et al. Proteinase-3 and myeloperoxidase serotype in relation to demographic factors and geographic distribution in anti-neutrophil cytoplasmic antibody-associated glomerulonephritis. Nephrol Dial Transplant. 2019;34(2):301–8.
Jablonski NG, Chaplin G. The evolution of human skin coloration. J Hum Evol. 2000;39(1):57–106.
Frederick JE, Slusser JR, Bigelow DS. Annual and interannual behavior of solar ultraviolet irradiance revealed by broadband measurements. Photochem Photobiol. 2000;72(4):488–96.
Gatenby P, Lucas R, Swaminathan A. Vitamin D deficiency and risk for rheumatic diseases: an update. Curr Opin Rheumatol. 2013;25(2):184–91.
Ramagopalan SV, Handel AE, Giovannoni G, Rutherford Siegel S, Ebers GC, Chaplin G. Relationship of UV exposure to prevalence of multiple sclerosis in England. Neurology. 2011;76(16):1410–4.
O’Connell K, Sulaimani J, Basdeo SA, Kinsella K, Jordan S, Kenny O, et al. Effects of vitamin D(3) in clinically isolated syndrome and healthy control participants: a double-blind randomised controlled trial. Mult Scler J Exp Transl Clin. 2017;3(3):2055217317727296.
Fisher SA, Rahimzadeh M, Brierley C, Gration B, Doree C, Kimber CE, et al. The role of vitamin D in increasing circulating T regulatory cell numbers and modulating T regulatory cell phenotypes in patients with inflammatory disease or in healthy volunteers: a systematic review. PLoS One. 2019;14(9):e0222313.
Penna G, Amuchastegui S, Giarratana N, Daniel KC, Vulcano M, Sozzani S, et al. 1,25-Dihydroxyvitamin D3 selectively modulates tolerogenic properties in myeloid but not plasmacytoid dendritic cells. J Immunol. 2007;178(1):145–53.
Kallenberg CG. Antineutrophil cytoplasmic autoantibody-associated small-vessel vasculitis. Curr Opin Rheumatol. 2007;19(1):17–24.
Zimba O. 338. VITAMIN D DEFICIENCY IN PATIENTS WITH GRANULOMATOSIS WITH POLYANGIITIS. Rheumatology. 2017;56(Issue suppl_2):kex062.340. https://doi.org/10.1093/rheumatology/kex062.340.
Korkmaz FN, Ozen G, Unal AU, Can M, Tuglular S, Direskeneli H. AB0530 25-hydroxyvitamin D3 deficiency: prevalance and its impact on disease activity in small-medium vessel vasculitis. Ann Rheum Dis. 2016;75(Suppl 2):1086.
Kalsch AI, Peters A, Buhl B, Breedijk A, Prem K, Schmitt WH, et al. Retinoid X receptor beta polymorphisms do not explain functional differences in vitamins D and A response in antineutrophil cytoplasmic antibody associated vasculitis patients. Autoimmunity. 2009;42(5):467–74.
Simpson S Jr, Taylor B, Blizzard L, Ponsonby AL, Pittas F, Tremlett H, et al. Higher 25-hydroxyvitamin D is associated with lower relapse risk in multiple sclerosis. Ann Neurol. 2010;68(2):193–203.
Laursen JH, Søndergaard HB, Sørensen PS, Sellebjerg F, Oturai AB. Vitamin D supplementation reduces relapse rate in relapsing-remitting multiple sclerosis patients treated with natalizumab. Mult Scler Relat Disord. 2016;10:169–73.
Kemna MJ, Tervaert JWC, Broen K. Timmermans SAMEG, Van Paassen P, Damoiseaux JGMC: Seasonal influence on the risk of relapse at a rise of antineutrophil cytoplasmic antibodies in vasculitis patients with renal involvement. J Rheumatol. 2017;44(4):473–81.
Kelly D, Theodoratou E, Farrington SM, Fraser R, Campbell H, Dunlop MG, et al. The contributions of adjusted ambient ultraviolet B radiation at place of residence and other determinants to serum 25-hydroxyvitamin D concentrations. Br J Dermatol. 2016;174(5):1068–78.
O’Sullivan F, Raftery T, van Weele M, van Geffen J, McNamara D, O’Morain C, et al. Sunshine is an important determinant of vitamin D status even among high-dose supplement users: secondary analysis of a randomized controlled trial in Crohn’s disease patients. Photochem Photobiol. 2019;95(4):1060–7.
O’Sullivan F, Laird E, Kelly D, van Geffen J, van Weele M, McNulty H, et al. Ambient UVB dose and sun enjoyment are important predictors of vitamin D status in an older population. J Nutr. 2017;147(5):858–68.
Li X, et al. Genetically-predicted vitamin D status, ambient UVB during the pandemic and COVID-19 risk in UK Biobank: Mendelian randomisation study. 2020;medRxiv 2020.08.18.20177691. [pre-print]. https://doi.org/10.1101/2020.08.18.20177691.
Jennette JC, Falk RJ, Bacon PA, Basu N, Cid MC, Ferrario F, et al. 2012 Revised International Chapel Hill Consensus Conference Nomenclature of Vasculitides. Arthritis Rheum. 2012;65(1):1–11.
Bossuyt X, Cohen Tervaert JW, Arimura Y, Blockmans D, Flores-Suarez LF, Guillevin L, et al. Position paper: revised 2017 international consensus on testing of ANCAs in granulomatosis with polyangiitis and microscopic polyangiitis. Nat Rev Rheumatol. 2017;13(11):683–92.
Watts R, Lane S, Hanslik T, Hauser T, Hellmich B, Koldingsnes W, et al. Development and validation of a consensus methodology for the classification of the ANCA-associated vasculitides and polyarteritis nodosa for epidemiological studies. Ann Rheum Dis. 2007;66(2):222–7.
patientMpower is a smartphone application. https://apps.apple.com/ie/app/patientmpower/id1077787395. developed in 2019.
Trinity Health Kidney Centre: AVERT [https://www.tcd.ie/medicine/thkc/avert/]
Reddy BP, Houlding B, Hederman L, et al. Data linkage in medical science using the resource description framework: the AVERT model [version 2; peer review: 2 approved]. HRB Open Res. 2019;1:20. https://doi.org/10.12688/hrbopenres.12851.2.
Van Geffen J, Van Weele M, Allaart M, Van der AR. Royal Netherlands Meteorological Institute (KNMI). Dataset. 2017. https://doi.org/10.21944/temis-uv-oper-v2. Accessed August 2021.
TEMIS v2 [http://www.temis.nl/uvradiation/UVdose.html; version 2.0]
Houben E, Groenland SL, van der Heijden JW, Voskuyl AE, Doodeman HJ, Penne EL. Relation between duration of the prodromal phase and renal damage in ANCA-associated vasculitis. BMC Nephrol. 2017;18(1):378.
Shamseer L, Sampson M, Bukutu C, Schmid CH, Nikles J, Tate R, et al. CONSORT extension for reporting N-of-1 trials (CENT) 2015: explanation and elaboration. BMJ. 2015;350:h1793.
Bates D, Mächler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. Journal of Statistical Software. 2015;67(1):1–48. https://doi.org/10.18637/jss.v067.i01.
Kirkwood BRSJ. Essential medical statistics. 2nd ed. Oxford: Blackwell Science; 2003.
Fox JWS. An R companion to applied regression. 3rd ed. Thousand Oaks: Sage; 2019.
Fox J, Weisberg S. Visualizing Fit and Lack of Fit in Complex Regression Models with Predictor Effect Plots and Partial Residuals. Journal of Statistical Software. 2018;87(9):1–27. https://doi.org/10.18637/jss.v087.i09.
Ross AC, Manson JE, Abrams SA, Aloia JF, Brannon PM, Clinton SK, et al. The 2011 report on dietary reference intakes for calcium and vitamin D from the Institute of Medicine: what clinicians need to know. J Clin Endocrinol Metab. 2011;96(1):53–8.
Lionaki S, Blyth ER, Hogan SL, Hu Y, Senior BA, Jennette CE, et al. Classification of antineutrophil cytoplasmic autoantibody vasculitides: the role of antineutrophil cytoplasmic autoantibody specificity for myeloperoxidase or proteinase 3 in disease recognition and prognosis. Arthritis Rheum. 2012;64(10):3452–62.
Mohammad AJ, Segelmark M. A population-based study showing better renal prognosis for proteinase 3 antineutrophil cytoplasmic antibody (ANCA)-associated nephritis versus myeloperoxidase ANCA-associated nephritis. J Rheumatol. 2014;41(7):1366–73.
Watts RA, MacGregor AJ, Mackie SL. HLA allele variation as a potential explanation for the geographical distribution of granulomatosis with polyangiitis. Rheumatology. 2014;54(2):359–62.
Maclure M. The case-crossover design: a method for studying transient effects on the risk of acute events. Am J Epidemiol. 1991;133(2):144–53.
Reid D, Toole BJ, Knox S, Talwar D, Harten J, O’Reilly DS, et al. The relation between acute changes in the systemic inflammatory response and plasma 25-hydroxyvitamin D concentrations after elective knee arthroplasty. Am J Clin Nutr. 2011;93(5):1006–11.
The authors would like to thank all participants who are enrolled in the UKIVAS and Irish Rare Kidney Disease registries, without whom this research would not be possible. We also extend our sincere gratitude to all clinicians across the UK and Ireland who actively recruit AAV patients to these registries to enable ongoing research.
Members of the RKD consortium:
Tallaght Hospital: Mark Little*, Peter Lavin**, Catherine Wall**, George Mellotte**, Jennifer Scott**, Ted Fitzgerald**, Hannah O’Keefe**, Rachel Dilworth**, Pamela O’Neill# and Vicki Carr#
St James’ Hospital: Niall Conlon*, Brenda Griffin*, Donal Sexton* and Caroline Kosgei#
Mater University Hospital: Yvonne O’Meara*, Eoghan White** and Stephen Mahony**
St Vincent’s University Hospital: Eamonn Molloy* and John Holian**
Galway University Hospital: Matt Griffin*, David Lappin**, Conor Judge**, Sarah Cormican**, Blathnaid O’Connell** and Michelle Clince**
Limerick University Hospital: Liam Casserly*
Cork University Hospital: Michael Clarkson*, Michelle O’Shaughnessy**, Alyssa Verrelli**, Sinead Stoeman**, Fergus Daly**, Laura Slattery** and Aisling Murphy#
Beaumont Hospital: Mark Little*, Declan De Freitas**, Peter Conlon**, Mark Denton**, Carol Treanor**, Colm Magee**, Conall O Seaghdha**, Paul O’Hara**, Susan McGrath**, Brona Moloney**, Dean Moore**, Dearbhla Kelly**, Mary McCarthy**, Tamara Wanigasekera**, Ayanfeoluwa Obilana**, Claire Kennedy**, Dervla Connaughton**, Mark Canney**, Limy Wong** and Sarah Moran**
*Principal investigator; **co-investigator; #study coordinator
Dr. Jennifer Scott is a Wellcome-HRB Irish Clinical Academic Training (ICAT) Fellow, and this work was performed within the Irish Clinical Academic Training (ICAT) Programme, supported by the Wellcome Trust and the Health Research Board (Grant Number 203930/B/16/Z), the Health Service Executive, National Doctors Training and Planning and the Health and Social Care, Research and Development Division, Northern Ireland. MAL received funding from Health Research Board/Irish Nephrology Society (MRCG-2016-12) and Science Foundation Ireland (13/RC/2106_P2 and 11/YI/B2093). ANG and EH received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 813545. The funders were not involved in any part of the study design, analysis or writing of the manuscript.
Ethics approval and consent to participate
Ethical approval for the RKD and UKIVAS components was granted by the Tallaght University Research Ethics Committee (ref 2019-08 List 29 (07)) and the NHS National Research Ethics Service (REC 10/H1102/77), respectively. All participants provided written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flowchart of participant selection for the UKIVAS cohort (1). Fig. S2. Flowchart of participant selection for the RKD cohort (2). Fig. S3. Participant locations at diagnosis for i). UKIVAS and ii). RKD cohort. Shading represents the absolute number of participants recruited to the study in that region, at the time of their diagnosis. Tendency to cluster around vasculitis referral centres is observed. Fig. S4. Distribution of vitD-UVB (kJ/m2) (a and b) and CW-D-UVB (kJ/m2) (c and d) by month and season: CW-D-UVB peaks in late summer and nadirs in winter. Fig. S5. The relationship between latitude (degrees) and a). annual daily average vitD-UVB (kJ/m2, Pearson correlation coefficient (corr) -0.98, p <0.001), b). annual daily average CW-D-UVB (kJ/m2, corr -0.96, p <0.001). CW-D-UVB at diagnosis c). not adjusted for season (kJ/m2, corr -0.12, p <0.001), and d-g). stratified by season (corr winter -0.25, spring -0.15, summer -0.47, autumn -0.22, p <0.001 for all). Fig. S6. The relationship between latitude (degrees) and a). annual daily average vitD-UVB (kJ/m2, Pearson correlation coefficient (corr) -0.98, p <0.001), b). annual daily average CW-D-UVB (kJ/m2, corr -0.96, p <0.001). CW-D-UVB at diagnosis c). not adjusted for season (kJ/m2, corr -0.12, p <0.001), and d-g). stratified by season (corr winter -0.25, spring -0.15, summer -0.47, autumn -0.22, p <0.001 for all). Fig. S7. AAV diagnosis: a). Latitude (degrees), b). average annual vitD-UVB (kJ/m2), c). average winter vitD-UVB (kJ/m2) and d). CW-D-UVB (kJ/m2) stratified by AAV phenotype (GPA and EGPA (“Not MPA”) vs. MPA), at diagnosis. ANCA-associated vasculitis (AAV), Microscopic polyangiitis (MPA), Ambient UVB dose at wavelengths than induce vitD synthesis (vitD-UVB), Cumulative-weighted UVB dose (CW-D-UVB). Fig. S8. AAV diagnosis: Latitude (degrees), average annual vitD-UVB (kJ/m2), average winter vitD-UVB (kJ/m2) and CW-D-UVB (kJ/m2) stratified by AAV serotype (a-d: MPO vs. Not MPO), at diagnosis. ANCA-associated vasculitis (AAV), Myeloperoxidase (MPO), Cumulative-weighted UVB dose (CW-D-UVB). Fig. S9. Ratio of cumulative-weighted UVB at diagnosis relative to the average value at the participant’s location: a). cohort 1 & 2 and stratified by b). AAV phenotype and c). AAV serotype. ANCA-associated vasculitis (AAV), Microscopic polyangiitis (MPA), Myeloperoxidase (MPO), Cumulative-weighted UVB dose (CW-D-UVB), Not statistically significant (ns). If a true association between UVB-predicted vitD status (CW-D-UVB) and AAV phenotype/serotype at diagnosis exists, one would expect that the ratio of CW-D-UVB at diagnosis relative to the average CW-D-UVB at that location would be significantly <1 for non-MPA disease subtypes. A one-sided Wilcoxon rank sum test was used to evaluate if the ratio was significantly different from 1. While the median ratio trended below 1 in the pooled analysis and when stratified by AAV phenotype and serotype, this was not statistically different from 1. Table S1. Induction and maintenance treatment of RKD cohort (2). Intravenous (IV). Table S2. Multi-level model investigating the factors associated with AAV relapse risk, stratified by AAV phenotype (sensitivity analysis). N (individuals) differs from N (observations) as multiple observations (remission +/- relapse) per individual were included, according to each participant’s disease course. The OR refers to the probability of having an AAV relapse (relative to remission). Model 1 investigates the effect of average winter (Dec-Feb) vitD-UVB (2004-2019), adjusted for age at diagnosis, gender, AAV phenotype, ANCA serotype and treatment, restricted to MPA phenotype. Model 2 investigates the effect of average winter (Dec-Feb) vitD-UVB (2004-2019), adjusted for age at diagnosis, gender, AAV phenotype, ANCA serotype and treatment, restricted to Not MPA phenotype. Cumulative-weighted UVB dose (CW-D-UVB), Standard deviation (SD), Microscopic polyangiitis (MPA), Myeloperoxidase (MPO), Odds ratio (OR), 95% Confidence interval (95% CI). Table S3. Multi-level model investigating the association between CW-D-UVB and the risk of AAV relapse (with treatment breakdown). N (individuals) = 439, N (observations) = 2080 as multiple observations (remission +/- relapse) per individual were included according to each participant’s disease course. The OR refers to the probability of having an AAV relapse (relative to remission). Cumulative-weighted UVB dose (CW-D-UVB), Standard deviation (SD), Microscopic polyangiitis (MPA), Myeloperoxidase (MPO), Odds ratio (OR), 95% Confidence interval (95% CI), Azathioprine (AZA), Mycophenolate mofetil (MMF), Methotrexate (MTX). Table S4. Uni- and multivariable logistic regression analysis of factors associated with AAV serotype at diagnosis, in the combined UKIVAS and RKD cohort. Note: OR (95% CI, p). The OR refers to the probability of having MPO-positivity (ref: not MPO-ANCA) at diagnosis. Model 1 investigates the effect of latitude, adjusted for age at diagnosis, gender and ethnicity (observations 1260, 1053 missing ANCA serology, AIC 1626.9). Model 2 investigates the effect of average annual vitD-UVB, adjusted for age at diagnosis, gender and ethnicity (observations 1260, 1053 missing ANCA serology, AIC 1626.8). Model 3 investigates the effect of average winter (Dec-Feb) vitD-UVB (2004-2019), adjusted for age at diagnosis, gender and ethnicity (observations 1260, 1053 missing ANCA serology, AIC 1626.5). Model 4 investigates the effect of CW-D-UVB at symptom onset, adjusted for age at diagnosis, gender and ethnicity (observations 1260, 1053 missing ANCA serology, AIC 1627.3). Ambient UVB dose at wavelengths than induce vitD synthesis (vitD-UVB), Cumulative-weighted UVB dose (CW-D-UVB), Standard deviation (SD), Myeloperoxidase (MPO), Akaike information criterion (AIC), Odds ratio (OR), 95% Confidence interval (95% CI). Table S5. Uni- and multivariable logistic regression analysis of factors associated with AAV phenotype in the combined UKIVAS and RKD cohort, restricted to White participants (sensitivity analysis). Note: OR (95% CI, p). Observations = 2142 (1 missing age). The OR refers to the probability of having MPA-AAV (ref: not MPA) at diagnosis. Model 1 investigates the effect of latitude, adjusted for age at diagnosis and gender. Model 2 investigates the effect of average annual vitD-UVB, adjusted for age at diagnosis and gender. Model 3 investigates the effect of average winter (Dec-Feb) vitD-UVB (2004-2019), adjusted for age at diagnosis and gender. Model 4 investigates the effect of CW-D-UVB at symptom onset, adjusted for age at diagnosis and gender. Cumulative-weighted UVB dose (CW-D-UVB), Standard deviation (SD), Microscopic polyangiitis (MPA), ANCA-associated vasculitis (AAV), Akaike information criterion (AIC), Odds ratio (OR), 95% Confidence interval (95% CI). Supplementary Methods.
About this article
Cite this article
Scott, J., Havyarimana, E., Navarro-Gallinad, A. et al. The association between ambient UVB dose and ANCA-associated vasculitis relapse and onset. Arthritis Res Ther 24, 147 (2022). https://doi.org/10.1186/s13075-022-02834-6