Study populations
Figure 1 shows a flow diagram of the study cohorts. The WLGP incidence cohort consisted of 923,941 individuals contributing 4,391,444 PYAR (mean = 4.8 years, SD = 3.2 years) and the PEDW incidence cohort consisted of 958,603 individuals contributing 4,545,876 PYAR (mean = 4.7 years, SD = 3.2 years). The mortality cohort consisted of 465,242 individuals, contributing 3,746,991 years of person-time (mean = 8.1 years, SD = 3.1 years), of whom 1,416 died during the 10-year follow-up period.
Table 2 summarises the proportion of the incidence cohorts with codes for SUD only, MD only or CC in their history. In the WLGP cohort, 75.4% were NC, 21.8% were MD only, 0.8% were SUD only and 2.0% were CC. 70.4% of individuals with SUD also had a code for MD and 8.4% of individuals with MD also had a code for SUD. In the PEDW cohort, 94.8% were NC, 2.9% were MD only, 0.5% were SUD only and 1.9% were CC. 79.7% of individuals with SUD also had a code for MD and 38.8% of individuals with MD also had a code for SUD.
Table 2
WLGP and PEDW cohorts by condition group
|
|
Total
|
n
|
% (95% CI)
|
WLGP
|
% of total with NC
|
923941
|
696691
|
75.4 (75.3–75.5)
|
|
% of total with MD only
|
923941
|
200981
|
21.8 (21.7–21.8)
|
|
% of total with SUD only
|
923941
|
7778
|
0.8 (0.8–0.9)
|
|
% of total with CC
|
923941
|
18491
|
2.0 (2.0–2.0)
|
|
% of MD with SUD
|
219472
|
18491
|
8.4 (8.3–8.5)
|
|
% of SUD with MD
|
26269
|
18491
|
70.4 (69.8–70.9)
|
PEDW
|
% of total with NC
|
958603
|
908363
|
94.8 (94.7–94.8)
|
|
% of total with MD only
|
958603
|
27985
|
2.9 (2.9-3.0)
|
|
% of total with SUD only
|
958603
|
4515
|
0.5 (0.5–0.5)
|
|
% of total with CC
|
958603
|
17740
|
1.9 (1.8–1.9)
|
|
% of MD with SUD
|
45725
|
17740
|
38.8 (38.4–39.2)
|
|
% of SUD with MD
|
22255
|
17740
|
79.7 (79.2–80.2)
|
Table 3 summarises the condition groups of the 923,941 individuals present in both the WLGP and PEDW incidence cohorts, by sex and across both settings (primary care and hospital admission). Overall, a greater proportion of females than males had a record of MD in either setting (26.0, 95% CI 25.9–26.2 compared with 17.0, 95% CI 16.9–17.1), whereas more males than females had a record for SUD or CC (1.0, 95% CI 1.0-1.1 compared with 0.6%, 95% CI 0.6–0.6 for SUD and 3.9%, 95% CI 3.9-4.0 compared with 3.1%, 95% CI 3.0-3.1 for CC). A higher proportion of males than females had no record of a condition in either setting (78.0%, 95% CI 77.9–78.1 compared with 70.3%, 95% CI 70.2–70.4).
A greater proportion of females than males had only a primary care record with an MD (22.6%, 95% CI 22.4–22.7 compared with 14.9%, 95% CI 14.9–15.0). In the WLGP SUD only and CC groups, the proportion of males with no PEDW record was greater than that for females; 0.8% (95% CI 0.8–0.9) compared with 0.4% (95% CI 0.4–0.4) for SUD only and 1.3% (95% CI 1.3–1.4) compared with 0.8% (95% CI 0.8–0.9) for the CC group. Across both sexes, 86.2% of the WLGP MD only group (173,167 out of 200,891), 74.0% of the WLGP SUD only group (5,754 out of 7,778) and 53.7% of the WLGP CC group (9,933 out of 18,491) had not had a relevant PEDW admission. Of the 696,691 individuals in the WLGP NC group, 5,754 (0.8%) were MD only in PEDW, 1,505 (0.2%) were SUD only in PEDW and 4,148 (0.6%) were CC in PEDW.
Table 3
Comparison of WLGP and PEDW by condition group
|
|
MALE
|
FEMALE
|
TOTAL
|
|
|
n
|
% (95% CI)
|
n
|
% (95% CI)
|
n
|
% (95% CI)
|
Across both settings
|
MD only
|
78415
|
17.0 (16.9–17.1)
|
120408
|
26.0 (25.9–26.2)
|
198823
|
21.5 (21.4–21.6)
|
|
SUD only
|
4809
|
1.0 (1.0-1.1)
|
2719
|
0.6 (0.6–0.6)
|
7528
|
0.8 (0.8–0.8)
|
|
CC
|
18166
|
3.9 (3.9-4.0)
|
14140
|
3.1 (3.0-3.1)
|
32306
|
3.5 (3.5–3.5)
|
|
No condition
|
360096
|
78.0 (77.9–78.1)
|
325188
|
70.3 (70.2–70.4)
|
685284
|
74.2 (74.1–74.3)
|
WLGP: MD only
|
PEDW: MD only
|
6828
|
1.5 (1.4–1.5)
|
13074
|
2.8 (2.8–2.9)
|
19902
|
2.2 (2.1–2.2)
|
|
PEDW: SUD only
|
660
|
0.1 (0.1–0.2)
|
1329
|
0.3 (0.3–0.3)
|
1989
|
0.2 (0.2–0.2)
|
|
PEDW: CC
|
2641
|
0.6 (0.6–0.6)
|
3282
|
0.7 (0.7–0.7)
|
5923
|
0.6 (0.6–0.7)
|
|
PEDW: No condition
|
68792
|
14.9 (14.8–15)
|
104375
|
22.6 (22.4–22.7)
|
173167
|
18.7 (18.7–18.8)
|
WLGP: SUD only
|
PEDW: MD only
|
63
|
0.0 (0.0–0.0)
|
35
|
0.0 (0.0–0.0)
|
98
|
0.0 (0.0–0.0)
|
|
PEDW: SUD only
|
157
|
0.0 (0.0–0.0)
|
112
|
0.0 (0.0–0.0)
|
269
|
0.0 (0.0–0.0)
|
|
PEDW: CC
|
1110
|
0.2 (0.2–0.3)
|
547
|
0.1 (0.1–0.1)
|
1657
|
0.2 (0.2–0.2)
|
|
PEDW: No condition
|
3904
|
0.8 (0.8–0.9)
|
1850
|
0.4 (0.4–0.4)
|
5754
|
0.6 (0.6–0.6)
|
WLGP: CC
|
PEDW: MD only
|
888
|
0.2 (0.2–0.2)
|
1062
|
0.2 (0.2–0.2)
|
1950
|
0.2 (0.2–0.2)
|
|
PEDW: SUD only
|
380
|
0.1 (0.1–0.1)
|
342
|
0.1 (0.1–0.1)
|
722
|
0.1 (0.1–0.1)
|
|
PEDW: CC
|
3352
|
0.7 (0.7–0.8)
|
2534
|
0.5 (0.5–0.6)
|
5886
|
0.6 (0.6–0.7)
|
|
PEDW: No condition
|
6096
|
1.3 (1.3–1.4)
|
3837
|
0.8 (0.8–0.9)
|
9933
|
1.1 (1.1–1.1)
|
WLGP: No condition
|
PEDW: MD only
|
2795
|
0.6 (0.6–0.6)
|
2959
|
0.6 (0.6–0.7)
|
5754
|
0.6 (0.6–0.6)
|
|
PEDW: SUD only
|
748
|
0.2 (0.2–0.2)
|
757
|
0.2 (0.2–0.2)
|
1505
|
0.2 (0.2–0.2)
|
|
PEDW: CC
|
2976
|
0.6 (0.6–0.7)
|
1172
|
0.3 (0.2–0.3)
|
4148
|
0.4 (0.4–0.5)
|
|
PEDW: No condition
|
360096
|
78.0 (77.9–78.1)
|
325188
|
70.3 (70.2–70.4)
|
685284
|
74.2 (74.1–74.3)
|
|
Total
|
461486
|
|
462455
|
|
923941
|
|
Of the 1,416 individuals in the mortality cohort who died during follow-up, 1020 (72.0%) were male and 396 (28.0%) were female. 607 (42.9%) were NC, 417 (29.4%) were MD only, 60 (4.2%) were SUD only and 332 (23.4%) were CC (0.2% of the NC group, 0.3% of MD only, 0.9% of SUD only and 1.2% of CC). Of the 165,835 individuals with MD and/or SUD, 809 (0.5%) died during follow-up.
Incidence
Figure 2 and Fig. 3 summarise trends in CC incidence rate per 1,000 PYAR between 2008 and 2017 by sex, age and WIMD quintile, presented separately for WLGP and PEDW. Table 4 summarises the incidence of CC by sex, age, WIMD and year, including IRRs derived from the Quasi-Poisson regression.
Overall incidence in WLGP significantly reduced over the period (2.5/1,000, 95% CI 2.3–2.6 in 2008 and 2.1/1,000, 95% CI 2.0-2.2 in 2017, IRR = 0.9, 95% CI 0.8-1.0, p = 0.01). Incidence in PEDW was stable (2.3/1,000, 95% CI 2.1–2.4 in 2008 and 2.2/1,000, 95% CI 2.0-2.3 in 2017, IRR = 1.0, 95% CI 0.9–1.1, non-significant).
Incidence for males (WLGP = 2.5/1,000, 95% CI 2.5–2.6; PEDW = 2.4/1,000, 95% CI 2.3–2.4) was significantly higher than for females (WLGP = 2.1/1,000, 95% CI 2.0-2.1; PEDW = 1.9/1,000, 95% CI 1.9-2.0), with IRR = 1.2 (95% CI 1.2–1.3, p = < 0.01) for both WLGP and PEDW. Incidence among females in WLGP (but not PEDW) declined whereas for males it remained stable; incidence from WLGP for females in 2008 was 2.3/1,000 (95% CI 2.1–2.5) and in 2017 was 1.8/1,000 (95% CI 1.6-2.0).
Higher incidence was significantly related to increasing age: incidence in WLGP increased from 0.2/1,000 (95% CI 0.1–0.2) for 11–14 year olds to 3.8/1,000 (95% CI 3.7–3.9) for 22–25 year olds (IRR = 25.1, 95% CI 20.1–31.9, p < 0.01); and in PEDW from 0.6/1,000 (95% CI 0.6–0.6) for 11–14 year olds to 2.8/1,000 (95% CI 2.7–2.9) for 22–25 year olds (IRR = 4.7, 95% CI 4.2–5.2, p < 0.01). The association between higher incidence and increasing age was stronger for primary care than for hospital admissions, with rates in WLGP lower than in PEDW in the youngest age band but higher in the oldest; this was evident in greater IRRs in WLGP between age bands.
Higher incidence was associated with greater deprivation; the lowest incidence rates were among the least deprived quintile (WLGP = 1.1/1,000, 95% CI 1.1–1.2; PEDW = 1.2/1,000, 95% CI 1.2–1.3;) with the highest among the most deprived quintile (WLGP = 3.7/1,000, 95% CI 3.6–3.9; PEDW = 3.2/1,000, 95% CI 3.1–3.3), IRR (WLGP) = 3.3 (95% CI 3.0-3.7) and IRR (PEDW) = 2.6, 95% CI 2.4–2.9, with rates declining in the intermediate quintiles as deprivation reduced. Between 2008 and 2017, the gap between most and least deprived quintiles reduced considerably in WLGP, with a significant reduction in the most deprived quintile and a smaller, non-significant increase in the least deprived quintile; in 2008, incidence was 4.5/1,000 (95% CI 4.1–4.9) in the most deprived quintile and 1.0/1,000 (95% CI 0.8–1.2) in the least deprived quintile; by 2017 incidence was 3.0/1,000 (95% CI 2.7–3.4) in the most deprived quintile and 1.4/1,000 (95% CI 1.1–1.7) in the least deprived quintile. This was not observed to the same extent in PEDW.
Table 4
Incidence of CC and IRR by sex, age band, WIMD and year
|
|
WLGP events
|
PEDW events
|
|
|
Incidence (95%CI)
|
IRR (95%CI)
|
p value
|
Incidence (95%CI)
|
IRR (95%CI)
|
p value
|
Sex
|
Female
|
2.1(2.0-2.1)
|
ref
|
|
1.9(1.9-2)
|
ref
|
|
|
Male
|
2.5(2.5–2.6)
|
1.2(1.2–1.3)
|
0.00
|
2.4(2.3–2.4)
|
1.2(1.2–1.3)
|
0.00
|
Age band
|
11–14
|
0.2(0.1–0.2)
|
ref
|
|
0.6(0.6–0.6)
|
ref
|
|
|
15–17
|
1.2(1.2–1.3)
|
8.1(6.4–10.4)
|
0.00
|
2.2(2.1–2.3)
|
3.6(3.2–4.1)
|
0.00
|
|
18–21
|
3.3(3.2–3.4)
|
22(17.6–27.9)
|
0.00
|
2.8(2.7–2.9)
|
4.8(4.3–5.4)
|
0.00
|
|
22–25
|
3.8(3.7–3.9)
|
25.1(20.1–31.9)
|
0.00
|
2.8(2.7–2.9)
|
4.7(4.2–5.2)
|
0.00
|
WIMD
|
1 (least depr)
|
1.1(1.1–1.2)
|
ref
|
|
1.2(1.2–1.3)
|
ref
|
|
Quintile
|
2
|
1.6(1.5–1.7)
|
1.4(1.3–1.6)
|
0.00
|
1.6(1.6–1.7)
|
1.3(1.2–1.5)
|
0.00
|
|
3
|
2.1(2.0-2.2)
|
1.9(1.7–2.1)
|
0.00
|
2.1(2.0-2.2)
|
1.7(1.6–1.9)
|
0.00
|
|
4
|
2.8(2.6–2.9)
|
2.4(2.2–2.7)
|
0.00
|
2.5(2.4–2.6)
|
2.0(1.8–2.2)
|
0.00
|
|
5 (most depr)
|
3.7(3.6–3.9)
|
3.3(3.0-3.7)
|
0.00
|
3.2(3.1–3.3)
|
2.6(2.4–2.9)
|
0.00
|
|
n/a
|
0.5(0.3–0.7)
|
0.3(0.2–0.6)
|
0.00
|
1.1(0.9–1.4)
|
0.8(0.6-1.0)
|
0.11
|
Year
|
2008
|
2.5(2.3–2.6)
|
ref
|
|
2.3(2.1–2.4)
|
ref
|
|
|
2009
|
2.6(2.4–2.7)
|
1.0(0.9–1.2)
|
0.69
|
2.3(2.2–2.4)
|
1.0(0.9–1.1)
|
0.88
|
|
2010
|
2.5(2.3–2.6)
|
1.0(0.9–1.1)
|
0.67
|
2.1(2.0-2.2)
|
0.9(0.8-1.0)
|
0.20
|
|
2011
|
2.2(2.1–2.4)
|
0.9(0.8-1.0)
|
0.06
|
2.2(2.1–2.3)
|
1.0(0.9–1.1)
|
0.54
|
|
2012
|
2.3(2.2–2.4)
|
0.9(0.8-1.0)
|
0.11
|
2.2(2.0-2.3)
|
0.9(0.8–1.1)
|
0.31
|
|
2013
|
2.4(2.2–2.5)
|
0.9(0.8–1.1)
|
0.34
|
2.3(2.2–2.5)
|
1.0(0.9–1.1)
|
0.82
|
|
2014
|
2.2(2.1–2.3)
|
0.9(0.8-1.0)
|
0.03
|
1.9(1.8-2.0)
|
0.8(0.7–0.9)
|
0.00
|
|
2015
|
2.2(2.1–2.3)
|
0.9(0.8-1.0)
|
0.03
|
2.2(2.1–2.3)
|
1.0(0.9–1.1)
|
0.58
|
|
2016
|
2.1(2.0-2.2)
|
0.8(0.7-1.0)
|
0.01
|
2.0(1.8–2.1)
|
0.9(0.8-1.0)
|
0.02
|
|
2017
|
2.1(2.0-2.2)
|
0.9(0.8-1.0)
|
0.01
|
2.2(2.0-2.3)
|
1.0(0.9–1.1)
|
0.59
|
Mortality
Figure 4 summarises observed unadjusted mortality rates for each condition group per 1,000 PYAR, by sex, age at start of follow-up and WIMD quintile. The highest rate was for individuals with CC (1.4/1,000 PYAR, 95% CI 1.2–1.5), followed by those with SUD only (1.1/1,000, 95% CI 0.9–1.4); these rates were not significantly different but both were significantly higher than rates for MD only (0.4/1,000, 95% CI 0.3–0.4) and for NC (0.3/1,000, 95% CI 0.2–0.3). Rates were significantly higher for males than females for all condition groups except SUD only, and were significantly higher for 18–25 year olds than 11–17 year olds for all condition groups except SUD only. Rates for the most deprived WIMD quintile were higher than any of the other quintiles, but other than in the NC group (most deprived = 0.3/1,000, 95% CI 0.3–0.4, least deprived = 0.2/1,000, 95% CI 0.2–0.2) there were no significant differences by deprivation other than between the most deprived quintile (1.6,1,000, 95% CI 1.4–1.9) and the second least deprived quintile (0.9/1,000, 95% CI 0.6–1.3) in the CC group.
Of 392 deaths among the SUD only and CC groups, we identified six who died in hospital with no prior history of SUD before their final admission. Reclassifying these as NC in the analysis made no significant difference.
Survival
Figures 5–11 show plots of Kaplan-Meier survival curves with p-values derived from Log Rank tests, by condition group, sex, age band and WIMD quintile. Due to risk of statistical disclosure arising from small counts, the curves for SUD only were excluded from Figs. 6–11. To further prevent statistical disclosure, age band and WIMD quintile were collapsed to two levels (11–17 and 18–25; least deprived 60%, or quintiles 1–3 and most deprived 40%, or quintiles 4 and 5).
Survival was significantly different for individuals with CC, NC or MD only, for both males and females (p < 0.0001, Fig. 5). Figures 6–11 show that survival for males was significantly lower than for females in all condition groups and in both age bands at p < 0.0001, and for 11–17 year olds with CC at p < 0.05. The group who were 18–25 at follow-up start had significantly lower survival for all conditions (all significant at p < 0.01) except females with NC where there was no significant difference by age. Results by WIMD group were mixed; survival for both males and females with NC was significantly lower for the more deprived group (females = p < 0.05; males = p < 0.001). Differences in survival between the least and most deprived females with MD only and CC, and between the least and most deprived males with CC were not significant; differences between the least and most deprived males with MD only were significant at p < 0.05.
Figure 12 summarises the results of a Cox regression with death from all cause as the outcome. Results showed that compared to the NC group, the risk of death during the study window was significantly higher for individuals with MD only (HR = 2.7, 95% CI 2.4–3.1, p < 0.001), with SUD only (HR = 4.5, 95% CI 3.4–5.9, p < 0.001) and with CC (HR = 8.7, 95% CI 7.5–10.0, p < 0.001).