Mortality in a cohort of 3.1 million children, adolescents and young adults.

BACKGROUND: Many youth deaths occur in the first year of life, from prematurity and anomalies. Detailing mortality after age 1 year may differentially guide preventive strategies in children, adolescents and young adults.
METHODS: A cohort study in Ontario, Canada comprised 3 139 698 children born from 1990 to 2016. Adjusted HR (aHR) for death between 1 and 24 years were generated, comparing demographic variables and parity.
RESULTS: After a median of 13.7 years of follow-up, 6930 deaths occurred between ages 1 and 24 years (incidence rate 17.0 per 100 000 person-years), peaking at age 23 years (43.7 per 100 000). The aHR for death was higher among males than females (1.44, 95% CI 1.37 to 1.51), rural versus urban areas (1.48, 95% CI 1.39 to 1.58), lowest versus highest income areas (1.39, 95% CI 1.29 to 1.51) and at parity 1 (1.16, 95% CI 1.10 to 1.23), parity 2 (1.34, 95% CI 1.23 to 1.45), parity 3+ (1.96, 95% CI 1.74 to 2.21), each relative to a child without an older sibling. Among males, the proportion of deaths due to injury jumped from 30% before age 15 years to 65% thereafter, and in females, from 28% to 51%. Intentional self-harm/assault explained 11% of injury-related deaths among males before age 15 years, and 20% thereafter, with respective figures of 18% and 17% for females. Deaths outside of hospital increased with age, from 35% at age 1 year, to 66% at age 22 years.
CONCLUSION: There is a heightened susceptibility of dying starting at age 15 years, especially among males, from injury, and arising outside of hospital. © Author(s) (or their employer(s)) 2020. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Introduction

The death of a child, adolescent or young adult is an untimely event, whose tragic loss is experienced by the parents, siblings, friends and caregivers of the deceased.1 Across the life course, males and persons of lower socioeconomic status are especially at increased risk of death,2 as are those affected by mental illness.3 Although rates of infant and youth mortality have been in decline for decades worldwide, ethnic disparities continue to exist in high-income countries, such as the USA, where American Indian and black youths have much higher rates of death compared with Asian, Hispanic and white youth.2 Even so, mortality rates in young and middle-aged white adults seem to have stagnated in the USA, especially in rural areas, due to increases in unintentional drug poisonings and suicide.4 5 Recent Canadian data observed a resurgence in adult deaths due to the opioid epidemic.6 However, relatively little work has focused on patterns of mortality in children and emerging adults.

A large proportion of childhood deaths occur in the first year of life—largely from causes which are not readily preventable such as a lethal congenital or chromosomal anomaly, prematurity or sudden infant death syndrome.7 8 Therefore, analysing childhood mortality after 1 year of age might provide insight about risk factors and causes of death—due to injury and otherwise—that are preventable, especially since the median age at which a child can walk unassisted is 12 months.9

Linkage of mortality to other health and demographic data may offer new opportunities for ‘precision public health’ approaches among identified higher risk groups. Hence, the current study capitalised on multiple linked health and demographic linked data from Ontario—Canada’s most populous province, which provides universal access to all primary and acute healthcare services. A novel birth cohort was created to explore associations of maternal immigration, time-varying income status, sibling order and cause of death among youths aged 1–24 years, further delineated by 1 year age groups.

Methods

A population-based birth cohort study was completed within the entire province of Ontario, Canada. Included were children born in Ontario between 1 April 1990 and 31 March 2016, and who were alive at age 366 days after birth. Other reasons for study exclusion are outlined in online supplementary files 1 and 2.

All hospital births were identified using administrative health databases, described elsewhere7 and detailed in online supplementary file 1. These datasets were linked using unique encoded identifiers and analysed at ICES. The linked datasets include the Ontario Office of the Registrar General Vital Statistics Deaths, which comprises mandatory completed death certificate information, coded according to the International Classification of Diseases (ICD), ninth revision, as well as the Ministry of Health’s Registered Persons Database, which contains vital status and sociodemographic information for all individuals ever eligible for the Ontario Health Insurance Plan. Maternal world region of origin was identified using the Immigration, Refugees and Citizenship Canada’s Permanent Resident Database, which contains data on country of citizenship for immigrants to Canada from 1985 onwards. Women not linked to this database were classified as Canadian-born/long-term resident. Linkage rates with the healthcare identifier for both datasets are high.10 Income quintile and rural residence were from Statistics Canada census data. The use of data in this project was authorised under section 45 of Ontario’s Personal Health Information Protection Act, which does not require review by a Research Ethics Board or informed consent.

Data analyses

The characteristics of the male and female children included were presented as either a proportion, median (IQR) or mean (SD), and contrasted by a standardised difference for each. For variables with missingness, which was always under 1%, we reported the proportion that were ‘missing’.

Mortality rates (per 100 000 children) and 95% CIs were calculated by 1 year age groups, starting at age 366 days after birth, and ending at the 24th completed year of age, while accounting for the changing number of children in the cohort at risk within each 1 year age interval. Mortality rates and unadjusted rate ratios (RR) at each 1 year age increment were then presented by sex (males vs females), rural vs urban residence (updated at each year of age increment) and income quintile (lowest (Q1) vs highest (Q5)) (updated at each year of age increment).

HRs for child death between ages 1 and 24 years were generated using marginal Cox regression, comparing 1) males versus females; 2) rural versus urban residence, time-varying at each year of age for a given child; 3) residential income Q1, Q2, Q3 or Q4 versus Q5, time-varying at each year of age for a given child and 4) maternal world region of origin (Caribbean, East Asia/Pacific, Hispanic America, Middle East/North Africa, South Asia, sub-Saharan Africa, Western Nations/Europe vs Canada/long-term resident/missing). All HRs were adjusted for maternal factors (World region of origin, age and liveborn parity (ie, number of prior live births) at the child’s birth) and child factors (year of birth, sex, preterm birth <37 weeks’ gestation, a congenital or chromosomal anomaly diagnosed before age 1 year—as well as time-varying residential income quintile and rurality, updated at each year of age). World region analyses were restricted to children residing in an urban area at age 1 year, since most immigrants to Canada settle within a city, where there is generally consistent access to primary healthcare. In each analysis, the proportion of missingness was described, and a corresponding adjusted HR was generated for that missing group.

As they are often preventable,11 the proportion (95% CI) of all deaths due to injury (defined below) was calculated by 1 year age groups, for males and females, separately. The trend in the proportion of injury-related deaths was assessed by age using a non-parametric Jonckheere-Terpstra test.12 All deaths due to injury were then presented by injury subtype: intentional (self-harm (ICD-9 E950-E959) or assault/homicide (ICD-9 E960-E969)) versus unintentional (transport accident (ICD-9 E800-E849), accidental poisoning/overdose (ICD-9 E850-E854) or other/intention unknown (ICD-9 E855-E949, E970-E999 and 800–999). We assumed that a death occurring outside of a hospital might be a proxy for an unanticipated death, such as by suicide or opioid overdose,13 and that this would be more common in the adolescent years. Hence, within each 1 year age group, we also calculated the proportion (95% CI) of all deaths occurring outside of hospital (ie, outside of an emergency department, inpatient setting, same-day surgery setting or psychiatric hospital) (online supplementary file 1). Therein, unadjusted RRs were calculated comparing the proportion of out-of-hospital deaths among males versus females within each 1 year age increment.

There is a documented higher risk of infant and childhood mortality with rising parity,14 15 but less is know about that risk in adolescence, especially deaths from injury. Accordingly, HRs for any death, as well as deaths due to injury or non-injury, were calculated by maternal parity (ie, the number of preceding liveborn siblings—1, 2 or 3+ vs 0). HRs were adjusted for the same covariates as in preceding marginal Cox regression models, except for liveborn parity.

Finally, in an effort to minimise unmeasured confounding due to early environmental or genetic factors,16 a sibling analysis was done comparing the risk of death among males versus females aged 1–24 years. From all available children who survived to age 1 year, and who had at least one sibling of the opposite sex also born in the study period, we randomly selected one boy and one girl sibling from an eligible family. Marginal Cox regression accounted for the matching process, and HRs were adjusted for maternal age and parity at the child’s birth, as well as for prematurity <37 weeks’ gestation and a congenital or chromosomal anomaly up to age 365 days.

This study is reported according to the Strengthening the Reporting of Observational Studies in Epidemiology guideline (online supplementary file 3).

Results

There were 3 400 353 births in the study period, of which 3 139 698 children (92.4%) were alive at age 1 year and residing in Ontario. Other reasons for exclusion are shown in online supplementary file 2. The characteristics of the children, and their mothers, are listed in table 1. There were 1 609 644 males and 1 530 054 females alive at age 366 days. About 7% of all children were born preterm, and 8.2% of males and 6.3% of females were affected by a congenital or chromosomal anomaly (table 1). Of all children, 11% were residing in a rural area at age 1 year.

Table 1

Characteristics of the children included in the cohort, and that of their mothers, contrasted by standardised differences

Characteristic	Male children(n=1 609 644)	Female children(n=1 530 054)	Standardised difference
Of the mother at the birth
Mean (SD) age, years	29.6 (5.4)	29.6 (5.4)	0.00
Median (IQR) number of previous live births	1.0 (0.0–1.0)	1.0 (0.0–1.0)	0.00
No previous live birth	790 081 (49.1)	751 990 (49.1)	0.00
Any previous live birth	564 351 (35.1)	536 816 (35.1)	0.00
Missing information on previous live birth	3782 (0.2)	3624 (0.2)	0.00
Maternal world region of origin
Canada/long-term resident/missing*	1 257 056 (78.1)	1 196 246 (78.2)	0.00
Sub-Saharan Africa	26 199 (1.7)	27 224 (1.7)	0.00
Caribbean	21 602 (1.4)	22 439 (1.4)	0.00
Middle East/North Africa	33 720 (2.1)	32 508 (2.1)	0.00
South Asia	98 887 (6.1)	93 168 (6.1)	0.00
Hispanic America	28 318 (1.8)	27 123 (1.8)	0.00
Western Nations/Europe	59 502 (3.7)	56 185 (3.7)	0.00
East Asia/Pacific	82 410 (5.1)	76 941 (5.0)	0.00
Of the child
Mean (SD) birth weight, g	3440 (585)	3322 (555)	0.21
Mean (SD) gestational age at birth, weeks	38.0 (1.8)	38.0 (1.8)	0.02
<28	3734 (0.2)	3418 (0.2)	0.00
28–31	10 089 (0.6)	8377 (0.6)	0.01
32–36	109 153 (6.8)	92 890 (6.1)	0.03
37+	1 486 668 (92.4)	1 425 369 (93.2)	0.03
Multiple birth	46 496 (2.9)	45 786 (3.0)	0.01
Any congenital or chromosomal anomaly diagnosed before age 1 year	131 134 (8.2)	96 264 (6.3)	0.07
Income quintile (Q) at age 1 year
Q1	350 902 (21.8)	334 349 (21.9)	0.00
Q2	319 496 (19.9)	303 669 (19.9)	0.00
Q3	324 957 (20.2)	309 245 (20.2)	0.00
Q4	331 985 (20.6)	314 879 (20.6)	0.00
Q5	273 672 (17.0)	259 521 (17.0)	0.00
Missing	8632 (0.54)	8391 (0.6)	0.00
Residence at age 1 year
Urban	1 429 301 (88.8)	1 358 953 (88.8)	0.00
Rural	178 576 (11.1)	169 416 (11.1)	0.00
Missing	1767 (0.1)	1685 (0.1)	0.00
Median (IQR) age at end of study follow-up, year	13.7 (7.7–20.4)	13.7 (7.7–20.4)	0.00

All data are presented as a number (%), unless otherwise indicated.

*Includes 170 women with missing information on immigrant country of origin.

The median (IQR) duration of follow-up in the study period was 13.7 (7.7–20.4) years, and there were a total of 40 733 595 person-years of follow-up. There were 6930 deaths between ages 1 and 24 years, for an overall incidence rate of 17.0 (95% CI 16.6 to 17.4) deaths per 100 000 person-years. The age-specific mortality rate declined from 30.8 per 100 000 at age 1 year to a low of 7.3 per 100 000 at age 10 years (figure 1), and rose thereafter to a peak of 43.7 (95% CI 37.7 to 49.7) per 100 000 at age 23 years.

Figure 1

Risk of death among a cohort of 3 139 698 Canadian children, adolescents and young adults aged 1–24 years. Shown on the Y-axis are age-specific mortality rates and 95% CIs, and below the X-axis the number at risk and the number of deaths within each age group.

Of all 6930 deaths at ages 1–24 years, 4237 (61.1%) were among males (20.3 per 100 000 person-years) and 2693 (38.9%) among females (13.6 per 100 000 person-years)—an adjusted HR of 1.44 (95% CI 1.37 to 1.51) (figure 2). Among all fatalities between age 1 and 24 years, the median (IQR) age at death was 12.4 (3.9–19.0) years for males and 8.2 (3.0–16.6) years for females. The unadjusted RR of death was higher among males than females at ages up to 3 years (figure 2), and again by age 14 years, where there was a notable rising separation in the risk of death, being 2.84 (95% CI 1.93 to 4.20) times higher among males than females at age 24 years.

Figure 2

Risk of death among male (blue) and female (red) Canadian children, adolescents and young adults aged 1–24 years. Shown on the Y-axis are age-specific mortality rates and 95% CIs, and below the X-axis are rate ratios for mortality comparing males vs females within each age group.

Compared with those within an urban area, children and youth residing in a rural area had an overall unadjusted HR of death of 1.60 (95% CI 1.50 to 1.70), and an adjusted HR of 1.48 (95% CI 1.39 to 1.58) (figure 3). At most ages, there was a higher risk of death among those living in a rural area (figure 3). There was a tendency for higher mortality with declining residential income (table 2). Compared with a child residing in the highest area income Q, those in the lowest income Q had an adjusted HR of 1.39 (95% CI 1.29 to 1.51) for death at age 1–24 years. At both younger and later ages, the unadjusted RR was significantly higher in Q1 than Q5 areas, rising to 1.97 (95% CI 1.26 to 3.08) by age 23 years (figure 4).

Figure 3

Risk of death among rural (blue) and urban (red) residing Canadian children, adolescents and young adults aged 1–24 years. Shown on the Y-axis are age-specific mortality rates and 95% CIs, and below the X-axis are age-specific rate ratios for mortality comparing those residing in rural vs urban areas.

Table 2

Death of a child, adolescent or young adult between ages 1 and 24 years in relation to their residential income quintile, time varying at each year of age

Residential income quintile (Q)	Number of deaths(incidence rate per 100 000 person-years)	HR (95% CI)
		Unadjusted	Adjusted*†
Q5 (highest)	1248 (14.7)	1.00 (referent)	1.00 (referent)
Q4	1310 (15.0)	1.04 (0.96 to 1.12)	1.02 (0.94 to 1.10)
Q3	1306 (15.9)	1.09 (1.01 to 1.18)	1.05 (0.97 to 1.14)
Q2	1319 (17.4)	1.18 (1.09 to 1.28)	1.12 (1.03 to 1.21)
Q1 (lowest)	1674 (22.4)	1.51 (1.41 to 1.63)	1.39 (1.29 to 1.51)

*Adjusted for maternal world region of origin, age and liveborn parity at the child’s birth, and the child’s year of birth, sex, preterm birth <37 weeks’ gestation, a congenital or chromosomal anomaly diagnosed before age 1 year and time-varying rurality.

†17 023 participants had missing information about residential income quintile. Their corresponding adjusted HR was 1.58 (95% CI 1.20 to 2.09).

Figure 4

Risk of death among Canadian children, adolescents and young adults aged 1–24 years residing in the lowest (blue) and highest (red) area income quintiles (Q). Shown on the Y-axis are age-specific mortality rates and 95% CIs, and below the X-axis are age-specific rate ratios for mortality comparing those residing in the lowest (Q1) vs highest (Q5) areas.

Compared with the children of a Canadian-born mother (15.4 per 100 000 person-years), the mortality rate was slightly higher among those whose mother was from sub-Saharan Africa (16.0 per 100 000 person-years), but not significantly so (adjusted HR 0.98, 95% CI 0.81 to 1.19) (table 3). Those whose mother emigrated from South Asia, Hispanic America, other Western Nations or Europe and East Asia/Pacific were at a significantly lower risk for death.

Table 3

Death of a child, adolescent or young adult between ages 1 and 24 years in relation to the child’s mother’s world region of origin, relative to Canadian-born women/long-term residents

Maternal world region of origin	Number of deaths(incidence rate per 100 000 person-years)	HR (95% CI)
		Unadjusted	Adjusted*†
Canada/long-term resident (n=2 111 160)	4696 (15.4)	1.00 (referent)	1.00 (referent)
Sub-Saharan Africa (n=53 171)	106 (16.0)	1.11 (0.92 to 1.35)	0.98 (0.81 to 1.19)
Caribbean (n=43 803)	87 (14.7)	1.00 (0.81 to 1.24)	0.89 (0.73 to 1.10)
Middle East/North Africa (n=65 942)	90 (12.2)	0.88 (0.71 to 1.09)	0.89 (0.72 to 1.10)
South Asia (n=191 263)	262 (12.0)	0.88 (0.77 to 1.00)	0.87 (0.76 to 0.99)
Hispanic America (n=54 080)	77 (11.3)	0.78 (0.62 to 0.98)	0.79 (0.63 to 0.98)
Western Nations/Europe (n=110 636)	143 (10.2)	0.70 (0.59 to 0.83)	0.73 (0.62 to 0.86)
East Asia/Pacific (n=158 250)	189 (10.1)	0.71 (0.62 to 0.83)	0.78 (0.67 to 0.90)

This analysis was restricted to children residing in an urban area at age 366 days after birth.

*Adjusted for maternal age and liveborn parity at the child’s birth, and the child’s year of birth, sex, preterm birth <37 weeks’ gestation, a congenital or chromosomal anomaly diagnosed before age 1 year and time-varying residential income quintile and rurality.

†170 participants had missing information about maternal world region of origin. They were assigned to the Canada/long-term resident category.

A total of 5648 deaths in the cohort aged 1–24 years had a documented cause of death. In males, the proportion of deaths due to injury jumped at about age 15 years, from a mean of 30% before, to 65% thereafter (p<0.001) (online supplementary file 4a). For females, there was a less pronounced rise at age 15 years, with a mean rate of 28% from age 1 to 14 years, and a mean rate of 51% from 15 to 24 years (p<0.001) (online supplementary file 4b). The greatest proportion of deaths due to injury was at age 19 years for males (77.2%) and at age 20 years for females (62.0%). Among males, intentional causes (self-harm or assault/homicide) explained 11% of injury-related deaths before age 15 years, and 20% of deaths from age 15 to 24 years (online supplementary file 4a). Accidental poisoning/overdose was the dominant cause of death due to injury in both male age groups (58% and 43%, respectively), followed by transport accident (31%). In females, intentional causes explained 18% of injury-related deaths before age 15 years, and 17% of deaths from age 15 to 24 years (online supplementary file 4b). Like males, accidental poisoning/overdose accounted for the largest proportion of injury-related deaths (46% and 44%, respectively), followed by transport accident (35% and 33%, respectively). Of all deaths in each age group, the proportion occurring outside of hospital increased with age, especially after age 15 years (figure 5).

Figure 5

Proportion of all deaths occurring outside of a hospital setting, among Canadian children, adolescents and young adults aged 1–24 years. Shown on the Y-axis are age-specific proportions and 95% CIs for deaths occurring outside of hospital, and below the X-axis are rate ratios for mortality comparing males vs females within each age group.

As the number of preceding liveborn siblings increased, so did a child’s overall risk of dying from any cause (table 4, upper). These HRs were more pronounced for deaths due to injury (table 4, lower) than from non-injury causes (table 4, middle).

Table 4

Death of a child between ages 1 and 24 years in relation to the number of liveborn siblings preceding the birth of the current child

	Number of preceding liveborn siblings	Number of deaths(incidence rate per 100 000 person-years)	HR (95% CI)
			Unadjusted	Adjusted*†
Death due to any cause	0 (n=1 542 071)	3490 (16.5)	1.00 (referent)	1.00 (referent)
	1 (n=1 101 167)	2353 (16.7)	1.05 (1.00 to 1.11)	1.16 (1.10 to 1.23)
	2 (n=349 689)	699 (17.5)	1.17 (1.08 to 1.27)	1.34 (1.23 to 1.45)
	3+ (n=139 365)	329 (23.9)	1.68 (1.50 to 1.89)	1.96 (1.74 to 2.21)
Death due to non-injury	0 (n=1 542 071)	1727 (8.2)	1.00 (referent)	1.00 (referent)
	1 (n=1 101 167)	1104 (7.8)	0.95 (0.88 to 1.03)	1.08 (1.00 to 1.17)
	2 (n=349 689)	339 (8.5)	1.03 (0.91 to 1.15)	1.22 (1.08 to 1.37)
	3+ (n=139 365)	151 (11.0)	1.30 (1.10 to 1.54)	1.62 (1.36 to 1.92)
Death due to injury	0 (n=1 542 071)	1111 (5.3)	1.00 (referent)	1.00 (referent)
	1 (n=1 101 167)	837 (5.9)	1.20 (1.10 to 1.31)	1.39 (1.27 to 1.52)
	2 (n=349 689)	219 (5.5)	1.24 (1.07 to 1.43)	1.60 (1.38 to 1.86)
	3+ (n=139 365)	105 (7.6)	1.96 (1.59 to 2.42)	2.77 (2.23 to 3.45)

Shown are deaths from any cause (upper), non-injury causes (middle) and injury causes (lower).

*Adjusted for maternal world region of origin, age at the child’s birth and the child’s year of birth, sex, preterm birth <37 weeks’ gestation, a congenital or chromosomal anomaly diagnosed before age 1 year and time-varying residential income quintile and rurality.

†7406 participants had missing information about the number of preceding liveborn siblings. Their corresponding adjusted HRs were 2.47 (95% CI 1.90 to 3.21) for death due to any cause, 3.70 (95% CI 2.73 to 5.02) for death due to non-injury and 1.50 (95% CI 0.81 to 2.78) for death due to injury.

In the 1:1 matched sibling analysis, there were 980 deaths among 565 588 sisters (13.3 per 100 000 person-years) and 1485 deaths among matched brothers (20.2 per 100 000 person-years)—an unadjusted HR of 1.52 (95% CI 1.40 to 1.65) and an adjusted HR of 1.46 (95% CI 1.35 to 1.58).

Discussion

In this large population-based cohort of >3.1 million live births—representing about 35% of Canadian live births between 1990 and 2016, a J-shaped relation was observed between age and all-cause mortality. Death in childhood and young adulthood was also more likely in rural and low-income areas, and among children whose mother originated from Canada, sub-Saharan Africa, the Caribbean or the Middle East/North Africa. Males tended to have a higher proclivity for death, especially from age 14 years onwards, nearly tripling that of females by age 24 years. The proportion of deaths occurring outside of hospital also increased with age. Deaths due to injury massively rose after age 14 years, and accounted for 77% of male deaths at age 19 years and 62% of female deaths at age 20 years. Among those aged 15–24 years who died of injury, intentional causes explained one in five of male deaths, and one in six female deaths, while nearly half of injury deaths were due to accidental poisoning or overdose.

Strengths and limitations

A strength of the current study was the inclusion of a large cohort within a universal healthcare system, enabling a representative and precise estimation and comparison of mortality across an array of age groups and covariates, as well as by causes of death. As the median age at the end of follow-up was 13.7 years, most deaths among children who died before adolescence would be captured. For deaths in late adolescence or young adulthood, fewer youth contributed to the denominator of those at risk (see figure 1, below the X-axis).

This study relied on death certificate data to determine cause of death, yet, the accuracy of injury-related causes of death has been questioned.17 Not all deaths herein had a documented cause, nor was it known which deaths were confirmed by autopsy or an on-scene investigation. In adults, accidental poisonings now largely comprise prescription and illegal drug overdoses,18 and it is at the discretion of a medical examiner or coroner to determine whether they are deemed ‘accidental’, based on the circumstances and evidence surrounding the death.19 20 Some suicides continue to be misclassified as accidental poisonings or a death of unknown intent, especially in the more recent era of opioid use.21 22 In 2013–14, of all intentional self-harm hospitalisations among Canadian youth aged 10–17 years, poisoning made up 88% of events among girls and 82% of events in boys.23 More than 50% of poisonings among Canadian youth presenting to hospital are intentional.24 Given the above, and the observations herein that accidental poisoning and overdose accounted for nearly half of all injury-related deaths, and that greater than half of all deaths occurred out of hospital by age 18 years, is it conceivable that more adolescent deaths were from suicide than realised.

The relation between some variables studied herein and the risk of death could still be explained by unmeasured confounding. However, at least for the comparison of mortality between males and females, it is unlikely that residual confounding persisted, since the sibling analysis generated a similar adjusted HR (1.46, 95% CI 1.35 to 1.58) as the non-sibling analysis (1.44, 95% CI 95% CI 1.37 to 1.51). Nevertheless, this study did not detail important factors in child development, including mental illness (eg, attention deficit hyperactivity disorder (ADHD) or bipolar disorder25 and educational attainment, nor were parental marital status or parental-child engagement assessed.

Mechanisms of disease

The elevated rate of death before age 1 year is known to be largely related to prematurity, major anomalies and sudden infant death syndrome.7 8 While those influential factors may have persisted herein over the second to fourth years of life (figure 1), they were accounted for in the multivariable models. By age 12 years, the youth mortality rate started to climb, especially from fatal injuries starting at about age 15 years. It is around this age that greater risk-taking behaviour begins, especially among males and individuals with ADHD or conduct disorder,26 as well as those with psychosocial factors that limit self-regulation.27 As succinctly summarised elsewhere, self-destructive behaviours typically arise from predisposing biological, personality and cognitive vulnerabilities, worsened by further exposure to negative life events and psychiatric disorders.28 In the USA, intentional injury has been correlated with degree of community-level poverty,29 30 as was similarly shown herein. Low-income status and parental education were each associated with poorer coping strategies among Danish adolescents.31 Consistent with a previous Swedish study,15 we found that the risk of all-cause mortality increased with the number of older siblings; however, the Swedish authors did not assess adolescent mortality at age 11–19 years or cause-specific mortality in children/youth due to small numbers.

Recent mortality data from the USA has shown a concerning rise in premature adult deaths due to drug poisoning, especially among white people and those residing in rural counties.5 Herein, mortality followed a similar pattern, with a 50% overall relative higher mortality in rural areas, including a 237% relative increase within the 11-year age group. While the health of Canadian youth born to immigrants is generally better than their non-immigrant peers,32 33 the current study did not show such a protective effect among the offspring of mothers from the Caribbean, sub-Saharan Africa and the Middle East (table 3). We do not know how our findings might differ from those of children that arrived as immigrants to Canada. Thus, while there are clearly many factors that contribute to child death, prevention programmes should not lose sight of these sizeable groups, whose susceptibility to premature death may be rising over time.5

Implications for policy and clinical practice

These findings reinforce observations from the USA and elsewhere4: there is clearly a susceptibility to dying among specific subgroups of children and youth, with a predilection for death due to intentional and unintentional injury with rising age,34 especially at age 15 years. These data certainly suggest that injury prevention strategies should be initiated before age 15 years, including ongoing efforts to reduce pedestrian, passenger and driver injury,35 36 for land vehicles, and for watercraft as well.37

A concerning high proportion of injury deaths during adolescence and early adulthood are from intentional causes, including self-harm. Others have shown that many of these victims have an existing or unrealised psychiatric condition, especially a mood disorder or ADHD.38 As another highly susceptible group, a large Danish case-control study of adults aged 15–35 years showed that the risk of unnatural death was 5.3 times higher in those with a history of violence, 12.7 times greater with a history of self-harm and 29.4 times higher with a record of both behaviours.39 There may also be time-specific climbs in self-harm among youth, such as during periods of economic recession.40

In response, it has been suggested that multiple strategies are needed to reduce self-harm among youth. At the individual level, one approach includes cognitive behavioural therapy, to reduce rash reactivity to emotions and to improve self-regulation.41 At the broader public health level, examples include improving access to suicide prevention clinics, restricting the purchase of non-opioid over-the-counter analgesics and the quantity within each medication pack,40 and limiting adolescent access to firearms in high-availability regions.42

Given their high occurrence in childhood and young adulthood, more in-depth study is clearly needed on poisonings and overdoses, and the degree to which some are intentional, rather than accidental, in nature. The latter creates an important opportunity for improving our community processes in confirming the probability that suicide accounts for a larger than realised proportion of such deaths. Likewise, linkages of mortality data to other health and demographic data can facilitate more nuanced deconstruction of mortality rates, to enable more precise public health initiatives.

Conclusion

There is a heightened susceptibility of dying starting at age 15 years, especially among males, from injury, and arising outside of hospital.

Prior research on deaths in childhood largely included those occurring in the first year of life and thereafter, while omitting those occurring after age 20 years.

Most prior studies were done in jurisdictions where universal access to healthcare was not available, potentially biasing their findings on the frequency, causes, and nature of deaths in childhood, adolescence and early adulthood.

A novel birth cohort was created within a universal healthcare system, using linked population-based datasets.

This enabled the exploration of associations between deaths at age 1–24 years with maternal immigration, time-varying income status, sibling order and cause of death.

Valid measures (internally and externally) generated precise estimates of risk by age, and from 1 to 24 years.

There is clearly a susceptibility to dying among specific subgroups of children and youth, with a predilection for death due to intentional and unintentional injury with rising age.

In addition to initiation of injury prevention strategies at age 15 years, the focus should be on a disproportionate number of injury deaths during adolescence and early adulthood from self-harm and unrealised psychiatric conditions.

More in-depth study is especially needed on poisonings and overdoses, and the degree to which some are intentional, rather than accidental, in nature.